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I. Boundary Layer Growth on a Spinning Body 


By C. R. IntinawortH 
Department of Mathematics, The University, Manchester * 


[Received October 19, 1953] 


ABSTRACT 


The skin frictional force is calculated for a body of revolution in the 
early stages of a suddenly initiated (and subsequently constant) screw 
motion about the axis of symmetry. In particular the effect of spin on 
the time of separation of the boundary layer of a sphere in such a motion 
is determined. 


§ 1. INTRODUCTION 


ONE question that naturally arises in considering the boundary layer 
of a body in a screw motion (i.e. moving forward in the direction of some 
axis and at the same time spinning about that axis) is: how does the 
spin influence the position of separation of the boundary layer ? 
Of course, one can immediately give a qualitative answer to this question. 
Since spin imparts a centrifugal acceleration to fluid elements near the 
surface of the body, it produces a favourable contribution to the pressure 
gradient upstream of the maximum section of the body and an adverse 
contribution to the gradient downstream of this section, so that the effect 
of increasing the spin is always to move the position of separation nearer 
to the maximum section of the body. To give a quantitative answer to 
the question when the flow is steady would require rather laborious 
although straightforward calculations involving the integration of a 
series of pairs of simultaneous ordinary differential equations of second 
and third order. It is easier to answer the question when the flow has 
just been started up, and the boundary layer is in the initial stages of 
its growth. These conditions are just those that apply to the initial 
motion of a projectile fired from a rifled projector. For one can suppose 
that at t=0, the moment of firing, the projectile, a body of revolution, 


* Communicated by M. J. Lighthill, 
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is instantaneously given a forward velocity V and an angular velocity 2. 
However, the speeds involved in this external ballistic problem are SO 
high that compressibility effects would have to be taken into account 
in solving it. The problem that is investigated below is simpler than 
this: it is concerned with the case when V and 2 are sufficiently small 
for compressibility to be negligible. This type of investigation for a 
non-spinning body of revolution was originally carried out by Boltze 
(1908), and the present work is a simple extension of Boltze’s. 


§2. Bounpary LAYER GROWTH FOR A GENERAL Bopy oF REVOLUTION 


The independent variables to be used are t, x and y, where ¢ denotes 
time from the beginning of the motion, x distance from the nose of 
the body measured along a meridian curve whose orientation is fixed in 
space, and y distance along the normal to this curve. This means that 
the frame of reference moves forward with the body but does not rotate 
with it. The radius of the body at the station x will be denoted by R(z). 
With u, v representing the x and y components of velocity, and w the 
transverse component due to the spin, the equations for incompressible 
flow in the boundary layer are 


0 ra) 
Pe lild iy MeL heheh oo | § 3 ores 
Ou ou ou wdR dU o7u 
ae BG 01k ye dee oye ofS es (2) 
Ah i ae Ow uw aR 02w 
Ot ax | Nog WRRAeea O72c cin) oh te # Wes he ae (3) 


where U(x) is the velocity of the main stream, and yp is the kinematic 
viscosity of the fluid. 

It is supposed that the fluid motion is caused by suddenly starting 
the body to move at t=0 with constant linear velocity V in the 
direction of its axis and with constant angular velocity Q about its 
axis. The quantity U(x) is therefore the velocity the main stream 
would have at the station x if the body were fixed in a flow of constant 
velocity V. As Goldstein and Rosenhead (1936) explain in discussing 
boundary layer growth on a cylinder, the boundary layer is initially a 
vortex sheet which begins to thicken because of diffusion. Since the 
boundary layer growth is essentially a process of diffusion, the con- 
vection terms and the pressure gradient U dU/dx can bas neglected 
in (2) and (3) in the first place. Hence the first approximations w,, v 
w, to the velocity components satisfy vial 


0 

alRn)t S(Ro)=, 
Ou, Ou, 
TOE ye a ea 
Ow, wy 


ot Oye . . ° ° 7 A 5 (6) 
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The required solutions of (5) and (6) are 
Uy UE, . : . . . ° . . . e (7) 
W,= QR(1—E#), eens 42. 2 (8) 
UT 
where B=erf=22 | e* dz, .a=m 2 and n=y/2(vt)¥2. 
J0 
These satisfy u,=0, wy=QR at y=0, t>0 and u,=—U, w,=0 for 
large values of y. From (4), 
lid 
ye tet 
| 04= — 2(vt) Pas Ci) ome eee ss ays) (0) 
n 
where me H(z) dz=nH—«(1—e-"). Actually, (9) is not Bh entirely 
i d 
satisfactory expression for v, because it behaves like nies aE (RU) 
for large 7 instead of becoming very small. However, w, is practically 
equal to U and w, practically equal to 0 at moderate values of 7, for 
which v,=O[(vt)!/?]. 
The second approximate expressions for the velocity components are 
Uz +Us, Vy tVe, Wi +We, Where Us, Vo, We satisfy 


7) 0 
ae (Ru,)+ oy (Rv,)=0, Se “Oso > te PG (10) 
07s OU, Ou, Ou, .w,7 dk dU 
f "Oye Ot Ot UA ae +04 Gime deca: de? 0 6 (11) 
an 
OW, OW, dw, Ow, , uyw, dk (12) 


: sds. « .0baae Ulere <)  yaaeg aren 
The solutions of (11) and (12) may be written as 


U2dR 
Ug=t Gea Ait pp R Fg leet PRG Rk zfs) > 
dU aR 
W2=1 (ors, Jor t 2U = 22) 


in which f,,, 72; are functions of y satisfying 


foi fe ae! ae ] 


fig er he fae . ae 
2i(0 HEPA 

Y2i "+ 2ngoi’ “st ei see (14) 
iene — G2(9)=G2( 0 


Hee Rares denote differentiation with ann . n, and 
8.,;=H—E'F—1, 


S99 —EF, 
Seseary (1 —E), 
d 
me tie 


Page ee (oR (=) er. 
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Now, the two independent functions occurring in the complementary 
functions of eqns. (13) and (14) are 


j=1+ 27?, 
pene hied oe 
k=(1+2n?)E-+nE’. 

Also it is easy to verify that the particular integrals corresponding to 
the terms 1—E, (1—H)? and H'F, of which the s,,; and t,,; are linear 
combinations, are respectively 

I=E—-1, 

m= 27?H?+ 2nHE'+4H?4 2h—1, . + pe 
n=(n?+4)H2+4yHE'+ 30k’. 


and 


It follows that 


for=— (1+ $0") 9+(5 + §07)k—2l n, | 
Foo=3a7j-+(4—Ga*)k—n, 

= —(1—2a?)9+ (2—207)k— 
23 ( ay aie a2)k—m, ate 


and 
Jar= — 30°) — (33a) n, 
Jag Sa?) + (3—Sa?)k+ 21—2m-+n. 
Finally, from (10) we find 


d?U aU y 
cae ( =) ie dU 1dkR 


Vo. — Dy lize aiciee | ; 
ge [u dda rE: aa BR dg (tert 2hee) 


———— 


ce d?R 
+04, Tak Meat 2 Bea hig + 2.2? (Fe) bes | ; 
where h,(7) =|. Fei(2) dz. 
In the third approximation (the last one we shall consider) the velocity 


components are u,+u,+ 
1TUg Ug, Vy+Ve+v w ) 
satisfy , 2tV3, Wi+W2+W 3, where us and ws 


oes uz Ou du a 

Ou, Us __ : OF Ou w,wodR 

Oy? Ob: se lipgetiame ns Soy ar Acre: ee oe) 
"9 . 

O20, ws OW, Ow, OW» Ow OW, | UyWe+uw, dR 


Vv ay 
yh Oh oe Nps aye dy 1”? Sy io a ee 


The solutions of (18) and (19) are 


BAe) 
Us= tC dU | 
ie act fart 0 (F z) fat dx RaelotU fete 
x 


dR 
SAY pi ees dU dR 
(ie =“ Fost Org, aia Fafa} PURSE fay + u(F =) fsa 
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aU dU dR d*k 
Wa=?? | eu 5 hgs1 +2 (= i RG3.-+ QU — Tn dg 7 QU? da? 934 


dR ad?R dR 
se Ua alae y Jss+ QR? Te G36 $2 R(j =) tor |; 


where 
Sst’ +2 fai —8f3,= 483), (20) 
fs(0)=fs,( 0)=0, 
9J3i. + 2793; — 893 ;=4ts;, (21) 
93(0)=93;( 0) =0, 


in which 
833=LHfe,—E'hgy, 
S39 2H fy,— os ae 
S33= 3H fzo—F (for foe) —E (hart 2hoo), 
S834=Hfao—E'hgp, 
S835= —Hfe2—F' foo’, 
S3g= Efo3 —Ffo3’— 2(1—E)ga1, 
837 Efy3 —E'Ngs, 
S3g=Efo3 —P’f3’ — 2E"ho3 —2(1—E) a0, 

and 
t33 =H + Bho, 
ts9= —Pgoy'+ E'hai, 
ts3—=2(1—E) fo + E2901 +922) —F (21 +920) +H (hoi t+2hoo), 
t34= HJ 22+ Hho», 
t35—=2(1—E) foot Ho 22—LG 20’, 
t3g=L'ho;, 
tyy—= 2(1—E) fog + 2H ‘hos. 

The complementary functions of (20) and (21) contain the functions 


p=34+127n?+4n* and g=(3+12n?+4n*)H-+ (5n+273)H’, and so the 
solution of (20) is 


for= qa [14 |, Caesse de— | e*pledsale) deh —p | alesse) a |. 


The corresponding formula for gs; is obtained by substituting ¢,; for s3,. 
It would be rather laborious to evaluate the integrals involved in these 
expressions for f,; and g3,. However, we can determine the skin friction 
at the surface of the body with much less effort. For the longitudinal 
and transverse components of skin friction are given by 


="(Fe2(7) (3) 
cee oy senor t on non 
dw 1 fup\? (ow 
and t= (=) .=-3 a an ioe 
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where j4, p are the viscosity and density of the fluid respectively. Clearly 
one does not need f3,(7), 93;(y) but bi fs: (0), 93; (0) in order to obtain 
the skin friction. Now, fail G(ss,), 

37 3i/? 


93; i ae 
; 4 co ‘ 
where | G)= 3) tae) —P@)}6e\de 
Since ss; and ¢,; are linear combinations of terms of the types f.;, go; 


Efe; E9o: Ffoi, FGo;, E’he;, it is convenient to evaluate G(f,;) ete. 
From these quantities we easily obtain :— 


7 11 i) 


‘ 11 89 36 256 
fs2'(0)= eau (5 —— v3) o3— ae 05, 


o 


7 193 9 
fs3'(0)= gat (<4 ela Jes 5 V3) 08 1792 


1 eid ifsy: 3} 
fsa (0)=— =a > ie v3) o3— ea! 


fs6 (0)= ote 56+ 184/3)a3+ ee 


pel : 2 128 
37 ()=— gat (7 F v3) ee ae 


» 2 370 126 +6 
Ss (0)= (- ad pal v3) 3 — adie 


9 5 135 
and 
23 253 39 ‘ 
Ys1 = Tat (F -F v3) ae aes 
23 451 84 x 
A: (= Tat ( 5 V3) 08 — PE, 
1 1 491 i R1< 
Yoo (0O)= — nig 12 
w= get (— 4H V3) at ys, 


220 ‘1024 
: ae 
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35 DOD 64 
Y —————— ae ae Oe Sie 


13 58 
937 (0)= 9 %t (-3 ae 5 V3) ae Th Oa 
Also, H’(0)=2a and, from (15), (16) and (17), 


; 8 ’ 16 a 
for’ (0)= 20+ 3%, fog (0)=2a— zs fos (0) =4a—8a3, 
and 


99 


16 29 
Joi (0)= — 20+ Pe Joo (0) =2a— eS Os 


We may now summarize our calculations in the following formulae for 
the skin friction components :— 


U*dk 


7 ,=0-56419 (“)" “| u+ {1 4244 = 4-0-1512 =< 


dU 
+-0-7267602R 7} 1 — {0 05975U2 +0: a1osTU (F y 


dic 
a, 0287602 5 ie ea +-0-02005U35, eae —0-03154U8 5, (ae) 
=O: 1028692 RE +0: 0075922UR 
+0: 22073180 (F :) bat. Ae ie 1522) 


1g (ue\t” ue 
T,—=0-56419 (“?) Ez {0 15120R—— +0-6977 QU fe: t 


U dU dR 
+ {o. 0401 QU R+0-07162 (= R+1: ‘9568.2U = Fie cs 


x 


aR 1 (dR 
_ —0-0879.2U? 5 +1-6404.0U alae) 


2h dR\2 
+0-018260R? +0-1669823R (x) beg...) er (25) 


§3. Skin FRACTION ON A SPINNING SPHERE 


We are now in a position to calculate from (22) and (23) the skin 
friction on any particular body of revolution. We choose the sphere as 
a convenient example. In the early stages of boundary layer growth, | 
the velocity U(x) must be the velocity of inviscid flow round the body 
concerned. Thus, for a sphere, U(x)=3V sin w/a, where a is the radius 
of the sphere, this being the velocity at the surface in potential flow. 
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Also, R(z)=a sin za. It follows that 
T p= 0°84628 (“?) = V sin x/a[1+ {(2:3634-+-0-4845 18) cos x/a}& 


— {—(0-17955-+0-007592B) sin? x/a 
+ (0-48844-+ 0-323601) cos? x/a}é2+ . . .], (24) 


and 
1/2 
7,=0:56419 (‘““) Qa sin x/a{ 1+ {1-2732 cos x/a}& 


+ {(0-1077—0-018268) sin? x/a 
+. (8-2548-+0:16698) cos? a/a}e@t+...J, - - - + + + (28) 
where B= 22a?/V? and = Vi/a. 

The flow will separate from the surface along the circle of latitude at 
which 7,=0. If, in the first instance, we neglect the £ term in (24), we 
find that separation will first occur at w=az, the rear stagnation point 
on the sphere, when 

é=£,=(2:36344 0-484518)-1. 
Also, if we retain the é term, we find that separation will still occur 
first at e=a7, when €=€,, where €, is the positive root of 
(0-48844-+ 0-323601f) + (2:3634-+ 0-484518)é—1=0. 


Values of €, and &, for different 8 are as follows :— 


B g, ae 
0 0-423 0-391 
0-1 0-415 0-383 
1-0 0-351 0-322 

10 0-139 0-130 


The quantity ¢ is of course the number of radii through which the sphere 
travels before separation, and, as separation occurs behind the maximum 
section of the sphere in this problem for a non-rotating sphere, we expect 
that it will occur sooner when spin is introduced. The table gives an 
indication of how much sooner for different values of the spin. It 
suggests too that the time of separation can be found to the second 
approximation, whatever the spin, by calculating it to the first approxi- 
mation (€,) and reducing the answer by about 7%. This does not mean 
that one can dispense with the quadratic terms in the series for +. and 7 
Reference to (25) shows that the quadratic term in 7, is of the ae 
order of magnitude as the linear term when eae ; 
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II. On the Validity of Mott’s Theory of the Beta-Flow in Polycrystals 


By Pav Fre_tHam 
Standard Telecommunication Laboratories, Enfield, Middlesex * 


[Received November 3, 1953] 


SUMMARY 


A survey is given of the experimental evidence about the activation 
energy for transient creep. It is concluded that, if the Andrade equation 
is written in the form <?= Aft, the dependence of A on stress and tempera: 
ture is governed by the same constants as for steady state creep. 


Iv a recent paper Mott (1953) proposed a statistical theory of the Andrade 
beta-flow. The theory postulates a hardening of the material during 
flow, in agreement with the ‘ lamellar break-down ’ observed by Andrade 
(1914) in lead and subsequently studied by means of x-rays by Gibbs 
and Ramlal (1934) and other workers, and more recently by McLean 
(1953) in aluminium by optical means ; and, given this hardening, it is 
of a sufficiently general nature to be applicable to a wide range of crystalline 
and non-crystalline materials (e.g. celluloid, Jessop and Filon 1928) in 
which Andrade’s relation between strain and the cube root of the time 
has been observed. The theory leads to the equation 


de \1/3 
r= 0(F') : Wee eee it 


where C is a constant, independent of stress and temperature, ¢,,. is the 
(transient) strain and de/dt the strain-rate of the steady, second-stage, 
flow. 

The object of the present note is to present results of an investigation 
of the validity of this relation, based on data either published as such or 
derived from the flow-curves of a number of metals and alloys given in 
the literature. 

From a detailed study of the flow under constant stress of high-purity 
copper (99-98%), aluminium (99-98%) and cadmium (99-99%), and of 
three grades of aluminium of commercial purity, all polycrystalline, in 
the temperature range —196° to 140°c, Wyatt (1953) concludes that at 
low temperatures the flow curves fitted 


€,,=« log t+ ¢,, Me eels os sw eh) 
and at higher temperatures 
én =pbAe,, Ae eer reer 
while in the intermediate range the curves fitted 
€ p= Old te BUUe eal yea rs en eee) 
The «-flow obeyed an equation of state 
F(a, €,, de,/dt)=0, Set eo EO) 


* Communicated by the Author. 
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and was interpreted by him by an exhaustion mechanism; the f-flow 
was related by him to a recovery process. 

Although Wyatt did not specify the nature of this process, it appears 
to be associated with a capacity for inter-granular deformation, and 
B-flow practically free from the z-component would probably occur in 
metals only at temperatures (generally above about 4T neat °K), and 
stresses outside the domain of ‘equi-cohesion ’. The results for aluminium 
and magnesium at 20°C probably do not conform to these requirements ; 
they are given in the table alongside the remaining data for the sake of 
completeness. . 

The analysis upon which the determination of the results tabulated 
here was based depends on the assumption, verified experimentally in 
many cases (Novick and Machlin 1947, Fastov 1950, Johnson and Frost 
1952, Feltham 1953, etc.) that the dependence on stress and temperature 
of the steady flow-rate can be represented by the formula 


de : qa 
Ti =2A exp (—A/kT) sinh (iz) ; iv. eee 
or, for stresses high enough to permit the replacement of 2 sinh (qo/k7’) 
by exp (go/kT’), by 
de —H-+ qo 
Thame exp (=r). eee a (7) 


Here o is the applied tensile stress and A, H and q are constants. We 
write 


—H+qo=Q() ; 
in general, if the linear dependence of log € on o is not valid, we have 
de/dt=A exp {—Q(c)/kT}. 
The validity of Mott’s relation (eqn. (1)) then implies 


aes €7,=B exp {—Q,,(c)/kT }t™=Bt”, en ee} 


Qi, /Q=H ,,/H=4q,/q=m=}. S Q . a : (9) 


With reference to the table, in which the above ratios and measured 
values of the exponent of time, m, are given for a number of pure metals 
and for some alloys, the following comments must be made : 

(2) The assumption of the occurrence of stress fluctuations of random 
sign, of the type assumed by Mott in the derivation of eqn. (1), is probably 
pene to the truth in the case of well annealed, pure polyorystale and 
1omogeneous one-phase alloys. The results, given in the table, which 


relate to impure metals (i 
Pe oe p S (i.e. steels) must therefore be accepted with 


(6) According to Mott’s theor 


the p- | 
smoothly into the steady flo y the B-flow would be expected to go-over 


w after a certain ar r 
ae mount. of , hardening, | 
; Ree bare magnitude of the stress and the temperature ha 
a P aoe ne transient part of the flow should therefore be the sole 
w up to the point at which its rate becomes equal to that of the 
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11 
Metal #: 
eta Melting Flow So SSS Sar a ate era ears ee Reference 
temp. “c}| temp. °c | Experimental On|Q Hig \ geld 
Aluminium 20 035-045 
.995 0/ 660 150 0:50-0:55 | 0-33-0-41 | — ae I 
99°995% 
200 0:18-0:65 
| Al allo 
Risse | 1080 | 50-250 0:33 0:37 | oe I 
popper |" 1083 410 0-33 
pure. 0-36 sa te TH 
Lead ¢ 
Sure’ 327 17 0-33 0:37 — = III 
20 0:3 
Magnesium) 651 150 0:42-0:45 0:33-0:46 — — 1 
- 99-83% 200 0:75 
Ze 270 0-39-0°85 
Tron fa 
(alpha) 1535 444 0:33 0-36 — — Til 
pure ” 
Tron 
(gamma) 1535 | 950-1270 0-33 — 30 Or ieee Oh IV 
B9-99%, 
|Steel, plain 
| 0:05%C 1530 | 950-1300 | 0-18-0-43 — 1:00 | 0-94 IV 
| 0-23%%C 1520 | 950-1350 0:33 — 0-42 | 0:38 V 
0:44%C 1515 | 950-1300 | 0-26-0-44 — 0-45 | 0-53 IV 
0:56%C 1510 | 950-1150 | 0-25-0-43 — 0-40 | 0-61 IV 
79°C 1490 | 950-1220 | 0-39-0-65 — 0-48 | 0-53 IV 
1-15%C 1475 | 950-1170 | 0-37-0-59 — 0-97 1:11 RV 
Zine 
99-99% 419 20 — 0:42-0:63 I 


I. Crussard 1946. 


II. Johnson and Frost 1952. RR 59 contains from 1 to 2% each of Cu, Fe, 


Duration of tests up to 2000 hours. 


_ Andrade 1914. Duration of tests: about half an hour. 
_ Feltham 1953 a. Duration of tests: one minute to one hour approximately. 


-| - 


Ni, Mg. 


V. Feltham 1953 b. Duration of tests: one minute to one hour approximately. 


The 0:05 %C 
was a deoxidized acid open-hearth steel, 


open-hearth 


type. 


steel was of the non-deoxidized basic open-hearth type ; the 0:79%C steel 
the other steels were of the deoxidized basic 
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non-linear part of the flow should therefore 
and not by the commonly adopted form 


=e. -b BHesL yt. Olas |e ee 


In fact, if eqn. (10) is used, B sometimes appears to attain a constant 
‘saturation’ value above a given stress, while y continues to increase 


with stress. Now from eqns. (7), (8), (9) : 


0 log P/0g: 2 Ser in, 
dlogy/dc 4g 
and m could thus not remain constant if an increase in y were not accom- 

nied by a corresponding increase in f. 

a B- af ee Aerieen by Andrade (1914) for lead at 17°c, by the use 
of eqn. (10), are employed for the evaluation of m by means of eqn. (11), 
it is found to fall to zero at high stresses, while the measured exponent of 
the time maintains its value of approximately 1/3. However, if the total 
strain of the initial flow is plotted against ¢1/? (eqn. (3)), B-values are 
obtained which, substituted into (11), give a constant m-value of about 
1/3 at all stresses used. 

(c) Flow curves were not available for copper and alpha-iron, and 
consequently 8 and y values determined only by means of (10) were avail- 
able. The tabulated results for these metals were determined for stresses 
up to only 750 and 2150 kg/cm? respectively, for which the ‘ saturation ’ 
of 8 did not appear to have been attained. 

(d) The measured m-values for aluminium and magnesium, given by 
Crussard (1946), are also based on the assumed validity of (10); the 
remaining results were obtained by the use of (3). 

(ec) The values of the parameters A, H and q (eqn. (6)) in the case of the 
aluminium alloy, determined by the method described by the author 
(Feltham 1953 a) from curves given in a paper by Johnson and Frost 
(1952), were found to be A=0-15 per hour, H=13-9 keals/g atom, 
g/k=2-44°K/kg/em? which, except for the last, are different from those 
estimated by the authors of that paper. 

If the qualifications relating to the results given in the table are taken 
into account, the data seem to point to the essential correctness of Mott’s 
interpretation of the B-flow. 


steady flow. The initial, 
be represented by eqn. (3) 
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ABSTRACT 


From the observed momenta and angles of decays in flight, variables 
are constructed which are better suited for the analysis than the measured 
quantities themselves. These variables are 1/P, the reciprocal of the 
primary momentum, «, a combination of the longitudinal momenta of the 
secondaries, and either p,, the transverse momentum, or the quantity 
e=2p,/P. The regularities in the distribution of the new variables, to be 
expected from a two-body decay, are discussed. The decay scheme may 
be derived by comparison with the experimental distributions. 


INTRODUCTION 


A METHOD of analysing the dynamics of two-particle decays will be des- 
cribed here. The need for such a method was felt after the inhomogeneity 
of neutral V-particles had been discovered by the Pic du Midi group. 
In order to treat the experimental data it was essential to devise a method 
of separating the events into different decay schemes. 

An examination of the energy released in the decay, the Q-value, is in 
general not sufficient for this purpose. This test can show up inconsis- 
tencies between the experimental data and the assumed decay schemes, 
but it cannot reveal the reasons for these inconsistencies. The Q-values, 
therefore, provide no clear indication as to how such a discrepancy can be 
removed. 

A more detailed (graphical) representation of decays has been given by 
Blaton (1950), but it does not lend itself to a quantitative treatment as we 
shall see in §3. The transformation of variables discussed in this paper 
provides a suitable approach. The first analysis based on the new 
variables was made by Armenteros et al. (1951). Some of the properties 
of the quantity « have been employed by Leighton e¢ al. (1953) and 
Fretter et al. (1953) in the examination of their results, while a more 
extensive use of the new variables has been made to great advantage by 
Thompson e¢ al. (1952, 1953). 

It is essential to realize that the method is merely a convenient way of 
displaying experimental results. The analysis in terms of « and « cannot 
reveal more than has been put into it: conservation of momentum and 
energy. Its superiority over the Q-value method derives from the fact 
that the latter checks energy conservation only. 


* Now at the Laboratoire de l’Ecole Polytechnique, Paris. 
+ Communicated by Dr. G. D. Rochester. 
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I. Tut Decay SURFACE 


$1. THe Decay SURFACE IN TrRMs oF Two MoMENTA AND THE 
IncLUDED ANGLE 
(i) The Decay Relations 
We consider a decay M > m +m. If the decay scheme is paris, 
momenta P, p,, P2 in the laboratory system (L.S.) and the aie ie . fe ; 
b1, 62, P=bit $e, are completely determined when f, ie velocity o hei 
primary, and 6, the angle of ejection in the centre of mass sys S 
(C.M.S.), are given. We assume that three quantities have been nea a 
the momenta p,, p, of two particles and the angle da» included. ere 
p, and p, stand for any pair of the three momenta P, 1, Py. AS a decay 


Fig. 1 Fig. 2 


Decay in the laboratory Decay in the centre of mass system. 
system. 


is fixed by two parameters (8, 0), the measured quantities have to satisfy 
a relation of the formt 


EE —PaPo cos bab = (m,?m,?+ M?p*?)t, ee Pe no ef (1) 
with H=(m?+p*)!. The constant on the r.h.s. is fixed by the decay 
scheme : M is the mass of the primary and p* the momentum of the decay 
products in the C.M.S. : 

] 
(eS ay WEP (my +m)? LMP — (m,—m,)?}. 2 sat a eae 


Equation (1) stands for three ‘ decay relations’ : 


EE .—p py. Cos ¢= 3(M?—m,?—m,°), a" cet pee Ghee) 
EE,—Pp, cos ¢, = 3(M?+-m,2—m,"), . . . (1 bd) 
EE.—Pp, cos ¢, = 4(M*@—m2+m,2). .» . . (1 c) 


The first relation applies to a neutral, the other two to a charged, decay. 
In principle, three events—three sets of measurements Par Po» Pap—Should 
be sufficient to fix the parameters M, M4, My determining the form of the 
decay relation. Because of the errors of measurement, the decay scheme 


is, in practice, adjusted so that the decay relation is satisfied as well as 
possible. 


} In §7 (v) we deal with the case when only two angles are measured. 


{ Here, as elsewhere, masses and momenta are expressed in energy units. 
Velocities are referred to the velocity of light. 
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The decay relation may be visualized as a surface in the space of the 
three measured quantities p,, py, ¢,». The decay surface is fixed by the 
decay scheme. The parameters M, m,, m, have to be chosen so that the 
points representing the decay events lie as close to the surface as possible. 
The decay surfaces in the space p,, ~,, ¢4, ate, however, intricate 
structures.t It is impossible to guess, from a series of observations, the 
parameters M,m,, mg fixing the surface. The essence of the method to 
be proposed is the construction, from p,, p,, ¢,, (and from these only), of 
another set of three variables, in terms of which the decay surface becomes 
manageable. With the new variables, it will be possible to construct the 
decay surface directly from the observations, whereas using P,, Po Pad 
one can only check the decay relations (1) against an assumed decay scheme 
(by examining the constancy of the Q-values). 


(ii) Statistical Distribution 
Momenta and angles in the L.S. can be expressed as functions 


PalB, 6), Pol B; 6), Pavl(B; oD) ELSE hi Sn eed yn? (3) 


of £, the velocity of the decaying particle, and of 0, the angle of ejection in 
the C.M.S. The eqns. (3) are a parametric representation of the decay 
surface ; elimination of 8 and 6 results in the decay relation (1). 

From (3), the statistical distribution of events over the decay surface— 
corresponding to the (supposed) isotropy of decays—can be derived. 
The isotropic 6-distribution, 


w(6)d0=1sin0d0,  (0<@<180°),.... (4) 


gives, with (3), the distribution of decays along a line B=const on the 
surface. 

Again, the distribution over the p,, p,, $,,-surface is much too compli- 
cated for a detailed comparison with experiment. In contradistinction, 
the theoretical distribution in terms of the new variables is simple enough 
to allow of a sensitive test for isotropy in the C.M.S. 


§2. THE VARIABLES ¢, «, 1/P 


The real trouble with the quantities p,, P»), dy, is that their values are 
determined primarily by B, the velocity of the decaying particle. Our 
aim is to eliminate the dependence on B, as far as it is feasible to do so 
without any presumed knowledge of the decay scheme. ‘The dependence on 
M, m,, m, and 6 will then come to the fore, and this will make it easier 
both to design a decay scheme and to test the statistical distribution. 

(i) A quantity well fitted for the purpose is the transverse momentum 


Dp ee? a Oe eat nt) 


This variable is not affected by the transformation to the L.S. and is 
therefore independent of 8. It is obvious what information we expect to 


NS 2 eee ee ee ae ee 
+ A graphical representation of the surface (1 6) has been given by Stanton 
(1944). See also Butler (1951). 
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gain from p,: its maximum value is the momentum p* with which the 
secondaries are thrown out in the C.M.S. For the reasons presented in 
§3 we shall use, besides p,, the dimensionless variable 


2p, 2sin J, sin do (6) 
==) ———) — + oa ae . . . . . . 
Jie sin ¢ 


(ii) The longitudinal components, 
Pru = P1008 $1, Por = P2 COS po, 


contrariwise, take up the momentum of the primary completely 
(P=p1,+p.,) and thus depend strongly on f. The longitudinal com- 
ponents by themselves are, therefore, unsuited for the analysis, but they 
convey information about the masses of the secondaries by the way they 
share the incoming momentum. The heavier particle takes up, on the 
average, the greater part of P: 


pilP = 3J+m—m,?)/M*, PalP = 3(M?—m?+m,?)/M?. 
For large momenta (P> p*), the ratios p,,/P and p.,/P depend only slightly 
on B. We suggest as second variable 


a, LOFTON ew Di — De ae sin (62—¢) 

Co — P ye = 7 ‘exiniah ae Sy . . . (7) 
Its mean value is a measure of the mass difference of the decay products 
and should be equal to 


of = ts 


Depending on the angle of ejection, « will fluctuate about «* with an 
amplitude which, for high energies, approaches 
eS = 2p"/M. es 
(ili) As third variable, 1/P, the reciprocal of the primary momentum 
has proved to be a convenient choice. 
Thus we propose to analyse decay data in terms of the variables 
UP; Pp om eres RR 
or LP; ee Cs we Main. a ie 
defined by (5), (6) and (7), instead of the directly measured ones 
Disa: ¢ (neutral decay) or P, p,, ¢, (charged decay). 
The introduction of the set of variables (10) unifies the treatment of 


neutral and charged d. : 
one, namely, ; ecays. The three decay relations (1) reduce to 


2 1 
Pr A(1/ MP1 P% ° (a— oF )P = 8). ee 
or 2 
ee Ge eet 
(Pp tifa... . (118) 


with p*, «*, e* given by (2), (8), (9). 


The second 
through these quantities only. For t ndary masses m,, mz, enter 


he purpose of the proposed method 
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it is therefore more natural to define the decay scheme by p*, «*, ¢* or 
M, «*, «* rather than by M, m,, my. We shall devote § 4 to a discussion 
of the first two sets of parameters, before we analyse, in §§5, 6 and 7, 
the distribution of events over the three projections of the decay 
surface (11). 

But first of all, we wish to give another justification for the introduction 
of the new variables. We shall, in particular, see that the e, «-projection 
of the decay surface (11 5) allows of a simple geometrical interpretation. 


§3. THe Decay ELLIPse 


In the C.M.S., the end-points of the secondary momenta move on a 
sphere of radius p* and centre at the source S (fig. 2). The L.S. momenta 
Pi, P, are obtained, in a classical fashion, by vector addition of m,P/M 


Fig. 3 Fig. 4 


$1 


Decay in the L.S. Decay ellipse P=const. The plane e, «. 


and m,P/M, respectively. For a given primary momentum P, the 
end-points of p,, pz move, therefore, in the plane of decay, on circles 
of radius p*, but whose centres are now shifted from the source in the 
direction P by different amounts. To picture the decay in the L.S., 
one circle suffices, if the source is represented by two images S, and So, 
separated by P, In figs. 3 and 4 the picture is corrected for relativity : 


CZ 
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the loci P= const are ellipses | whose minor semi-axes are constant 
and equal to p*. The centre of the ellipse divides S,S, in the ratio 


x mz—mM~s" 
Tees Swi) oa orem 


This ratio is at the same time the average value of P1,/Po. for events on 
the upper half of the ellipse are as frequent as those on the lower half. 
The sources S, and Sy will be outside the ellipse as long as f, the velocity 
of the primary, is greater than both f,* and f,*, the velocities of the 
secondaries in the C.M.S. 

The advantages of such a method of visualizing the decays need not 
be stressed (Blaton 1950). However, for the display of decay data and 
their analysis the scheme is quite unsuited: if we fix, say, Sy, then S, 
and with it the ellipse recede from S, with increasing momentum P. 
Ellipses belonging to different momenta intersect. For a wide range 
of momenta, occurring in practice, the ellipses are elongated in one 
direction so much that, for the analysis, it is of no avail to know that 
they are ellipses. 

The remedy is to fix the distance S,S, by reducing the momentum 
seale of each event proportionally to 1/P. We choose as scale factor 2/P, 
fixing thus S;S,=2. From the centre of S,S,, which we make the origin, 
we have therefore to set out the vector 

9 ah 

plP1—4P) = Pete 
The transverse and longitudinal components of this vector are just 
the variables « and «, introduced in §2. In the ec, «-plane (fig. 5) each 
decay is represented in its original shape, all angles being preserved. 
The momentum scale has been chosen so that all events have the same 
‘size’. The rather abstract concept of the decay surface acquires, in 
the variables «, «, 1/P, a concrete meaning: the planes 1/P=const are 
faithful pictures of the decays, compatible with the P in question. 

After this mapping, the ellipses P=const (see §7 (iii) and fig. 13) are 
centred in a fixed point on S,S, and are confocal with eccentricity 
«*—=2p*/M. The axes decrease with increasing momentum, the vertical 
(major) axis as 1/P, the horizontal axis as 1/f. 

These are the facts which make the analysis in the «,«-plane so 
attractive, compared with the p,, a-diagram.t 


§4. Tot Decay PARAMETERS p*, €*6a™ AND M.e*.a™ 


The mass M of the unstable particle is, of course, greater than the 
masses of the secondaries together: M =M1+Mp. 


For most purposes 
we may, furthermore, restrict ourselves to M1>Mz. 


} As a consequence of the convention to t 
P=const cover only the right half of the ellipse. 

{In this projection, used by Thompson et al. (1952, 1953 
re-scaled only in the longitudinal direction, 


ake p, positive, the events 


), the decays are 
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(i) We denote quantities referring to the C.M.S. with an asterisk. 


The momenta and the energies of the secondaries are related to the 
masses by 


l 

ae ay UP (m+ me)?! MP — (m,—m2)?}! - + . (12a) 
1 

ie — IM (M?+-m,?2—m,?) ch Ja arp ech Bo: ae ee ars (12 b) 
1 

Tits a ay Am? +m,?), ‘ ‘ 2 . ° ° . . . (12 Cc) 


and the velocities follow then from 
ue meee bate pte Ao. <=.” (13) 
(ii) The essential features of a decay are shaped by the two dimension- 


less constants ¢«*=2p*/M and «*=(H,*—E,*)/M, which can be 
expressed by the mass ratios : 


a= (1- (5) | 1-A5) | ae dela") 


Fig. 6 


O 1 


e— 


Classification of decays according to «*, a*. 


The parameters <«*, «* determine completely the relations connecting 
dimensionless decay quantities, e.g. the velocities and the angles in the 
L.S. The number «* is a measure of the energy release. No energy is 
freed if «*—0. When «* reaches its maximum value 1, the whole rest 
mass of the primary is converted into kinetic energy. The parameter «*, 
on the other hand, characterizes the asymmetry of the decay. It runs, 
for a fixed «*, from «*=0 (symmetric decay, m,=m,) to a*¥—1—e%, 
when m,=0. All decays with m,>m, can thus be located in the triangle 
of fig. 6, in which the positions of some two-body decays are indicated. 
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(iii) The knowledge of e* and «* suffices for the calculation of the 

decay velocities in the C.M.S. : 
e* € 


Bx — Tae Bs = TO A ae CO (15) 


dition, either M or p* are given. 


The energy scale is determined when, in ad 
can then be calculated from 


Momenta and energies in the C.M.S. 


Meteo te OL ee hc 
pra ye Oa 
BY M tas) = pt. a - ae 
E* = = (1—a*) = p*, ae . » 3 ot Ee 


These formulae are to be compared with (12), (13) which express the same 
quantities in terms of M, m,, mg. 

(iv) The method proposed aims in the first instance at the determination 
of p*, «*, «* or M, «*, «*. With either of these sets the masses involved 
in the decay can be calculated from 


2 
Np Dan x Me RE 
M , 1 4 
ea “5 ker ok Sem lt =the (; <3 -1) oe 
1 
= M feet * see th i ‘ 
Min = = L a)’ ec. pe (z=—1) ae) 


with £,*, B,* from (15). 


Il. THe THREE PROJECTIONS oF THE Decay SURFACE 
§5. Tue Proszctions 1/P, « aND 1/P, p, 
eee Determination of p* 
ae nea ee, € is of interest only in combination with the 
es or visualizing the complete decay surface (fig. 11). 
Pilate a oak 
ip -p sind, (0<@<180°) 


the decay events will fill 
the re <i i 
hie aaa ne region between the axis «=0 and a straight 
aN si fe Cae une projection 1/P, p,, which leads to the determin- 
p", the C.M.S. momentum of the secondaries (fig. 7). Since 
A PP, = p* sin 6, (0<@<180°) 
nay ser i 
8 ne ea events should lie within the strip between p,=0 and p,=p* 
stribution in the horizontal direction is given ae he eee. 
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spectrum of the decaying particles. In the vertical direction, the 
distribution, assuming isotropy in the C.M.S., ist 
Dy Apr 


w(p;) dp, = Pp —pPl (18 a) 


which leads to an increasing number of events as p, approaches p*. 


For the determination of p* it may be better to use the integral 
probability distribution 


O\MEL py eel 
WO@>0/) =| 1—|-5 ee ee at (18:0) 
Fig. 7 Fig. 8 


§=90° 
9 = 60° 120° 


N¢ pe pt) 
0 =0° 180° 
ne 
A/P —>- Pe = 
The projection 1/P, p;. Number of decays with transverse 


momentum greater than p,°. 


The number of particles with a transverse momentum greater than p/°, 
plotted against p,°, will lie on an ellipse (fig. 8) whose horizontal semi-axis 
is p*. The broken line in fig. 7 corresponds to W=4 and divides the 
inhabited part of the plane into regions of equal population. 


§ 6. THe Progection 1/P, « 
(i) The Limiting Hyperbola 
Remembering the definition 
Pu-Pa sin (¢2—¢1) 


C= SSS 


iP sin ¢ 
one finds, by a Lorentz transformation, 
a=a*te*[1+(M/P)*}icos6 . . . . . (194) 
with a* == (m2—m,2)/M2,  «* = 2p*/M. 
The decays are distributed between the two branches (@=0, 180°) 
omg eee Nie tee me, Ghee ce (LOL) 


of a hyperbola, drawn in fig. 12. The position of the limiting hyperbola 
is fixed by the difference in mass of the secondaries, since its axis of 


+ For the influence of experimental errors on the p, distribution, see 
Barker (1953). 
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symmetry (0=90°) is given by x—a*. About this axis the decays are 
uniformly distributed (see (ii)) with an amplitude increasing as 1/f. 
The shape of the limiting hyperbola is determined by the energy release : 
the real semi-axis is equal to «*, and the slopes of the asymptotes are 
equal to +2p*. ‘ 

This slope is the same as that of the limiting line in the 1/P, e-diagram 
(fig. 11). The reason for this is clear : the space between the asymptotes 
is the region allowed by classical mechanics. The broadening to a 
hyperbola is a relativistic effect. According to classical mechanics the 
events with P= are concentrated in the point «=«* of the vertical 
axis, whilst, in fact, they extend over the segment FF. The transverse 
component, «, however, does not undergo a Lorentz dilatation. In 
classical mechanics, the decay ellipses of §3 are circles and the decay 
surface is (half) a circular cone, whose contours are given by the limiting 
line in fig. 11 and the asymptotes in fig. 12. This cone is, by relativity, 
deformed into the surface described in §7 (v). 


(ii) Statistics. Determination of «* and «* 


The assumption of isotropy in the C.M.S. leads to two consequences 
regarding the distribution of «-values for a given primary momentum P, 
ie. in the 1/P,«-diagram, the distribution of decays along a vertical 
strip: (a) For any P, the distribution is symmetric about a fixed value «*. 
(6) The distribution is rectangular with a width «*/8. Both these 
inferences can be checked independently and enable «* and «* to be 
calculated : 


g 
| 


: Siz (20) 
Fe*2+ g¥*2 — 92 _1(2p*/P)2, ey sae st ote ae (21) 


Comparison of these averages for different ranges of the primary 
momentum provides a check for isotropy. In (21) p* is assumed 
to be known from inspection of the 1/P, p-diagram.t The corresponding 
eqn. (25) of the next section is probably more useful. 

The over-all distribution of « is flat between «*—e* and a*+te* and 
tails off in a way determined only by the momentum spectra of the 
decaying particles. 

In principle, thus, the three decay parameters p*, a*, «* have been 
derived from the distribution of events in the planes 1 IP. p, and 1/P, « 
In §8, we shall indicate how the limitations of this necnedate due to 
experimental errors, can be overcome. 

Pea decay scheme may be checked in the 1/P,«-plane by 
sate i ying ee prvalues with those to be expected from the 
tae a ie or this purpose one draws the curves 6=const which, 
Rig nee é ), are hyperbolae. Since p,=p* sin 6, they are also the 
Gem teas a e corresponding to the assumed decay scheme. 

ty, the check may be done by means of the curves Pr Pz Or 


} Another method of finding <*, when p* 


o 


and «* are known, is described in §9. 
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p=const. Methods of constructing these loci are not too difficult to 
design. Other ways of examining a decay scheme will be described in § 9. 

(ui) The Events S,, S, 

These are the ‘ sources ’ of § 3, projected on to the decay surface. They 
are the points «=F1 of the limiting hyperbola (fig. 12) and represent 
decays in which one of the secondaries is at rest in the L.S. (p,=0 or 
P2=0). This can happen only when the primary has a velocity equal 
to the velocity of the secondary in the C.M.S. (B=f,* or B=8,*) and 


when the secondary is ejected in a direction opposite to the primary 
momentum (06=180° or 6=0): 


Sy: B = B,*, Pi = 9, Det i p* lm, o 
22a 
So: B = 8,*, Po = 9, Duet = pa. ) 


The secondary which does not come to rest emerges in the L.S. with 
the velocity 


1 mym.\2 |]  M*%—m,?—m,? 
z= [1+ (ai) | Se (22 b) 


The points S,, S, are important for many constructions. Some of their 
properties are summarized here (and in §7 (iii)) : 

The loci ¢,=const start from S,, the loci ¢,=const from S,. The 
horizontal lines through S, and S, enclose the region where both 
secondaries are, in the L.S., ejected in the forward direction. Above S3, 
particle 2, and below S,, particle 1, is thrown backwards. The probability 
of finding particle 1 or 2 ejected in the backward direction is, for a fixed 
velocity B of the primary, given by 4(1—f/8,*) and 4(1—8/f,*), 
respectively. For primary velocities+ greater than f,*, i.e. to the left 
of S,, both ¢d, and d=¢,+¢, go through a maximum. To the right 
of S,, the angle ¢ has a minimum. The loci é=const connect S, and S,. 

Furthermore, the points S,, S, play a réle in the construction of the 
curves p,=—const and p,=const. 


§7. THE PROJECTION e€, « 


(i) The third projection, already discussed in § 3, has as coordinates the 
dimensionless quantities { 


__ 2p, _ 2sin ¢, sin $2 
niger ye = sin (23) 
Vive tie sin ($2—$4) i 


— 


Le sin d 
The events fill the whole of the half-plane «>0. Energetic events will 
be found near the segment F',F', which represents, as in the 1/P, «-diagram, 


+ Sy is to the left of S,, ie. By* >B,*, because of m, > mz. 
+ An alternative is the plot p, «, used by Thompson et al. (1952, 1953). 
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decays with P=oo. For low energies of the primary the events recede 
to a distance ~M/P~=1/B from the centre of f,/'5. 

The position of an event in the e«, a-plane is, according to (23), 
determined by the angles only. This representation is therefore the 
appropriate choice when the momenta are not known, but two angles 
are measured. (In (iv) more attention will be given to this case.) 

Also, if only one angle is measured (together with the momenta), 
this projection possesses one striking virtue when the errors of measure- 
ment are taken into account: the position of an event is much better 
defined than in the other planes. The errors in the angle are usually 
negligible compared with the momentum errors. The event is, therefore, 
bound to lie on a line é=const (or ¢;=const, ¢.=const), i.e. on a circle 
through 8, and S, (or on a straight line through S, or S,) as is evident 
from fig. 9. Because of the momentum errors, the event will be 
represented by a segment of this line. In the other projections a decay 
event will be spread, by the momentum errors, over a quadrangle. 


Fig. 9 Fig. 10 


The projection e, a. Calculation of the momenta from 
the angles 41, do. 


In the «, a-projection—and only in this projection—the loci of constant 
Day Pz or p are independent of the decay scheme. The same is true for 
any dimensionless combination of measured quantities, e Poll 
P1/P2 which are constant along semi-circles Oe 


(11) Determination of «* and e* 
The distributi 
distribution of events over the « a-plane yield d 
1/P, «-projection, the dimensionless decay para tx ‘ “ mee 
s Tameters «*, e* ; 


oe 
a = y Be een ee 


1.43 ce es 
+E + % = g2— <2 
Cy x Cais . . 5 ° . . 5 (25) 
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Compared with (21), formula (25) has the advantage of referring to 
the plots in the ¢,«-plane only. Application of (21), employing the 
1/P, «-plane, depends on the previous scrutiny of the 1/P, p, diagram 
and is sensitive to errors in the determination of p*. 

In order to check an assumed decay scheme, the measured values of 
both P and p, may be compared with those expected. The curves 
P=const and p,=const, corresponding to the assumed decay scheme, 
serve this purpose. These loci are constructed below. They can, of 
course, also be used to evaluate, from the decay scheme, the momenta 
if these are not known. This will be done, however, in a more efficient 
way in (iv). 

(ili) P=const and 6=const (fig. 13) 

The lines P=const, 

(a—a*)? € 


[4CHPe Gps 8) 


2 


are confocal ellipses with foci F/',, #, and eccentricity «*=2p*/M. The 
major (vertical) and minor (horizontal) semi-axes are given by 


e*/B = e* . [1+(M/P)?}* and wp) Pie" ea Pe. 


the ratio of their lengths is y=(1—f?)-+. 

The curves P=const of fig. 13, together with fig. 9, make it possible 
to follow the behaviour of the angles ¢,, ¢,, ¢ for a given P. This 
behaviour will be different depending on whether the sources S,, S, are 
outside or inside the ellipse P=const. This happens when B>f,* and 
when 8<f,*, respectively, where 8,*, 8,* are the C.M.S. velocities of the 
secondaries (see footnote, page 23). The maxima of ¢, and ¢, are 
found by drawing the tangents from S, and S, on to the ellipse. 
Similarly, dnax and ¢yi, are found by constructing the circle(s) through. 
S, and S, tangential to the ellipse P=const. The points where, in the 
«,a-diagram, ¢, reaches a maximum lie on a semi-circle through S, 
and the point .=«*—e*B,* on the «a-axis (S, and «*+e*/6,* for ¢5).f 

Since p,—p* sin 6, the curves p,=const coincide with 6=const. They 
constitute a set of confocal hyperbolae, 

(x—a*)2 2 

Se0s2 Om = sin? 6 = eae ° . . . . . (27) 
with the same foci as the ellipses P=const, which they intersect at 
right angles. The asymptote of the hyperbola branch 6=const makes 
an angle @ with the positive direction of the a-axis. The loci 6=const 
can be used for testing the isotropy in the C.M.S. For instance, half 
the decays should fall between the two branches 6=60° and 6=120° of 
a hyperbola (fig. 13). 


+ The corresponding curves in the 1/P, a-diagram are parabolae through the 
same points, with vertex on the «-axis. 
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(iv) Momenta from Angles and Decay Scheme : 
| Dia e 
When the angles ¢,, ¢, are measured, the momenta Py Pa ae on 
derived if a decay scheme is assumed, Le. if p*, «* and e«* are ; 


Figs. 11, 12, 13 


ith li =c =const. 
The three projections of the decay surface with lines P Beer ber §=cons 
The statistical distribution of decays for constant P is indic f 


Projection 1/P,«. The decays are contained 
between the two branches (@=0, 180°) of a 
hyperbola. Population I of fig. 11 is to be 
found between the broken lines 6—60° and 
§6=120°, while population II is equally divided 
between the two marginal regions, 


We give below both a graphical and ar 


Projection 1/P, €. The 
decays are contained 
within the sector between 
the 1/P-axis and_ the 
straight line 6=90°. Half 
of the events (population 
I) lie above, the other half 
(population II) below, the 
broken line. 


Projection ¢«, «. Population 
I is again contained between 
the lines 0=60° and @=120°; 
population IT is distributed 
over the rest of the half-plane. 


1 analytical method of evaluating 


the primary velocity 6 and the ejection angle 6. Once, 6 and @ have 


been found from (28) or (29) 


. . SI] SI : 
P, Sin d; = pssind, =P. sin $ sin $2 


, the momenta follow from 


= p* sin 6, 
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For the graphical construction (fig. 10) one measures the distances 

ff, of the event from the points F’, and F, and finds 6 and 6 by means of 
1 aoe 

_ tithe jy ate 


aes” Der 


To follow this construction analytically,t first the angles 1, wo. are 
calculated : 


(28) 


eesca3 2 1 *)__ 
cs pee at nee Rew EE) a aod 
Borer ae Een a Sea 68 
Then f£ and @ are obtained from 
1 COs }(.—Y1) sin }(%2—7 1) 
5 ee ee 6 = ——= (296 
B ~ cos 2Fy,) OE Car Une 
or By=P/M and @ from 
M Re (sin %, sin #,) ees (sin %, sin fy) (29) 


Us COS 3 (y+) ” Sin 3(fo-+y1) 
The problem just treated arises when the momenta have not been 
measured, but when the origin of a neutral V-event has been located. 
In this case, however, it may be impossible to distinguish between the 
particles 1 and 2; the sign of « is then undetermined. Evaluation of 
f and @, according to (28) or (29), gives two different answers, but often 
one of them can be discarded. 


(v) The Decay Surface 


The three projections of the decay surface in the spacet 1/P,«, « are 
shown in figs. 11, 12, 13. The surface constitutes half an elliptical cone 
with axis parallel to the 1/P-axis. Towards 1/P=0 the semi-ellipses 
contract into a straight line. Hence the surface develops into a wedge 
whose thin end is formed by the segment F',F’, of the a-axis. Towards 
low momenta the cross section changes into a semi-circle, and the decay 
surface approaches the ‘ classical’ cone of § 6 (i). 


Ill. DETERMINATION OF THE DercAay SCHEME 


§8. THE Decay SCHEME 


In principle, the parameters p*, a*, —«* can be derived from the 
distribution of the decay events in the three projections, and thus the 
decay scheme determined from the observations without any pre- 
conceptions: the 1/P,p,;plane gives p*, the 1/P,« or the «,a yields 
a* and e*. This method, however, is only feasible with very good 


experimental material. 


+ Equivalent expressions have been obtained by M. M. May (Fretter 1953) ; 


formulae (29) are however simpler. 
+ Thompson et al. (1953) have also pointed out the usefulness of plotting 


the decay events three-dimensionally, in the space I/P, «, py. 
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In practice, with not very accurate, or a limited amount of, data, 
only the parameters p* and «* may be expected to be determined with 
reasonable accuracy. To fix the decay scheme, further information 
must be invoked. This is provided by assuming the masses of two of 
the particles to be known. 

Together with p*, «*, this additional information over-determines the 
decay scheme. Fixing two masses leads thus to a relation between 
p* and «*. This relation may be used to check the evaluation of p* and «*, 
and enables the values of all decay parameters to be adjusted within 
narrow limits. For decays into a proton and a 7-meson, the relation 
between p* and «* is illustrated in fig. 14. Besides «*, the energy 
release Q and the quantity «*+e«* are plotted against p*. 


Fig. 14 


—— 


Relation between the decay parameters for decays into a proton and a 7-meson 


The following procedure has been found to be very convenient for 
calculating the curves exemplified by fig. 14. The decay quantities are 
expressed in terms of two angles p and o: 

a* — COS p Cos a | 
«* = sin psing 


a) ; 
p* = $(m,+m,) sin p tana = 3M sin p sino at 
Q = (m,+m,) tan o tan do = 2M sin? to. 


With two of the masses oi 
Iven, p and 
following equalities : & p a are connected by one of the 


M1+Ms 27 
ot ee mee LT 4m, 2m 
p M1—™m, - COS a. = Tr 08 ¢ = cos oF. . (31) 


For a range of o-values, p is calculated from (31) and then the decay 
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parameters from (30). For each value of oc, one obtains a set of 
decay parameters compatible with the assumed mass values. 


§9. CHECK OF THE Decay SCHEME 


(i) A method of checking the decay events individually against an 
assumed decay scheme has been described in §§6 and 7. The measured 
values of p, and P are compared with those expected from the assumed 
decay scheme. For this comparison one draws the curves p,=const 
in the 1/P, «-projection or the curves p,—const and P=const in the 
e, x-diagram. 

(ii) Another test consists in taking two of the decay parameters as 
known and evaluating the third one for each event separately. This 
method can be applied in different ways : 

(a) Starting from the assumption that two masses are known, the 
third mass is calculated by means of the decay relation (1). This is the 
usual check for the constancy of the Q-values. 

(6) Alternatively one may start from p*, «*, calculate «*=2p*/M, 
and check against the value obtained, e.g. by the method of §8. The 
calculation of «* from the observed values of p,, « and P of a single decay 
is conveniently done in two steps: 


sin 0 = p,/p* ne aes (ol GC) 
#2 __ %—a* a 2p* 9 
6c (=) og) ee ears (32 b) 


This test is, however, better performed graphically in the 1/P, «-projection. 
With 6 from (32 a), one plots 


ay * 1 
—< against pee he (33) 


The event is thereby shifted on to the upper branch of the limiting 
hyperbola (cf. (326) and (196)). To the points (33), therefore, a 
hyperbola has to be fitted with an asymptote whose slope (2p*) is fixed 
from the outset. The value of e* is given by the length of the semi-axis. 

(c) Finally one may start by assuming M and «* to be known and 
again check the value of <* (Thompson, see Barker (1953)). By plotting 


a—a* ne € (34) 
(1+ W2/P2)3 agains MP ' 0 O 5 
the events are shifted, in the ¢, «-plane, on to a semi-circle of radius e* 


ict, (11 )). 
We have, however, a preference for the graph (33). It has the 
advantage of giving to the high energy events the greater weight they 


+ Low momenta reveal, indeed, the value of p* only—at least from momenta 
measurements. In the 1/P, «-graph this shows up in the fact that, for small P, 
the branches of the limiting hyperbola are widely apart: the ordinate of an 
event is of little help in determining the length «* of the semi-axis. The better 
position, in this respect, of the high energy events may be counteracted by 


the greater errors in the momenta, 
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carry in the determination of e*. But, more important, method (5) is 
still applicable if no assumptions regarding the masses are made: only 
the values of p* and «*, available directly from the analysis, are needed 


for the construction. 
(d) The method (b) can be exploited in another way by plotting 


—ay*\2 2n*\ 2 
(=) against (=) . «eee 


and fitting to the points a straight line of slope 1. This intersects the 
ordinate axis at a distance «*? from the origin. 
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ABSTRACT 


The half-lives, mass assignments and radiation characteristics of the 
neutron excess iridium isotopes have been determined as follows : 


194m, 195],y 196Typ 197Ty 1987p 
Half-life 7-3 SC 140 min 9-7 days 7min 50 sec 
f-energy (Mev) — 2:1, 1:25 0-08 1-6 3°6 
y-energy (Mev) 5-6, 0-125 0-42, 0-66, 0-58, 0-76, 1-8 0-78 


OS ee ke Omer = () 


The gamma-ray energies of 19-hour !%Ir have also been examined. 

The 7-second 1°4"Ir isomer is either a 125 kev E3 transition followed by 
the 5-6 Mev gamma-ray or possibly a 5-6 Mev M6, K6 or E7 transition 
followed by a 125 key M1, E2 or E3 gamma-ray. The latter alternative 
would be the highest spin-change gamma-decaying isomer yet discovered 
~ and both are unusual in the very high energy transition involved. 


§ 1. INTRODUCTION 


Wirs the partial exception of 19-hour '™Ir the neutron excess isotopes of 
iridium have not previously been accurately characterized. Butement 
(1951) found evidence for 1°°Ir or 197] with a half-life of somewhat above 
one hour, produced by a (y, p) reaction on platinum, using 23 Mev x-rays ; 
but the very low source intensity prevented any accurate study. While 
the work described below was in progress Christian, Mitchell and Martin 
(1952) found a 140-minute and a 7-minute iridium activity, the maximum 
beta-energy of the former being given as about one Mev and the activities 
being assigned to 1%°Ir and to !°°Ir or 1°7Ir respectively. This paper 
describes the production and characterization of the iridium isotopes of 
mass 194 to 198 inclusive. 


§ 2. EXPERIMENTAL TECHNIQUES 


The techniques for preparing the isotopes and studying their radiations 
were similar to those previously described (Butement and Shillito 1952). 


The isotopes were prepared by (n, p), (n, pn) and (y, p) reactions on 

platinum using neutrons of energies up to 48 Mev and 28 Mev x-rays from 

a synchrotron. For half-lives so short that no chemical separation of the 
* Communicate by the Authors. 
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iridium was possible, discs of platinum metal were irradiated. ei 
intermediate half-lives platinum chloride, and for long half-lives nee 
divided platinum metal, were irradiated, chemical separation oy a e 
iridium using an inactive carrier then being possible by standard methods. 


§ 3. RESULTS 
(i) Half-lives and Mass Assignments 


(n, p) and (n, pn) reactions in the cyclotron together produced activities 
with half-lives of 50 seconds, 7 minutes, 140 minutes, 19 hours, 9-7 days 
and 70 days in both beta- and gamma-decay, while a 7-second activity was 
observed in gamma-decay only. (n, p) reactions alone produced all the 
above activities except the 7-minute. (y, p) reactions produced the 
7-minute, 140-minute and 19-hour activities but the 50-second, 9-7-day 
and 70-day activities were not detected. All activities except the 7-second 
and 50-second were radiochemically identified as iridium. The 19-hour 
and 70-day activities correspond to the well known !Ir and }Ir. Of 
the remaining activities the 140-minute and 7-minute can be unambig- 
uously assigned to }%*Ir and 1°7Ir respectively while the 9-7-day and 50- 
second activities are most probably associated with 1°*Ir and 1°%Ir 
respectively. The approximate yields of the 19-hour, 140-minute and 
7-minute isotopes from (y, p) reactions on platinum were in the ratio 
8:4: 1 as compared with the relative abundances of the !°°Pt, 19*Pt and 
198Pt isotopes of 4-7 : 3-5 : 1 and to the yield ratio 5-3 : 3-7 : 1 obtained by 
Christian, Mitchell and Martin (1952). 

The 7-second gamma-decay was also observed when ‘ Specpure ’ 
iridium metal was irradiated in the cyclotron by neutrons of energy 
=—15 Mev and it must therefore be associated with 1%Ir, 192Ir or 1%[r, 
The yields of the 7-second and 50-second isotopes after (n, p) reactions on 
platinum were about the same, so that the 7-second decay cannot be 
assigned to °°Ir, which would be formed from the very scarce Pt. It is 
therefore either 1°Ir or Tr, Theoretical consideration of the radiation 
characteristics supports an assignment to 1% Tr (see below). 


(ii) Radiations 

(a) 50-second }°%Ir. The maximum beta-energy was found from the 
aluminium absorber curves to be 3:6--0-3 Mey, using Behrens’ nomogram 
(1950). Each point on the curve (fig. 1) was found from a separate 
irradiation. The gamma-ray energy was measured by the photographic 
scintillation spectrometer which showed one 0-78 Mev gamma-ray 
decaying with about the right half-life. Knowing the efficiencies of the 
counters, the ratio of beta- to gamma-transitions was found to be 1: 1 
and the gamma-rays were found to be in coincidence with the beta 
particles. 
The decay scheme is therefore simply a 3-6 Mev beta-transition followed 
y # 0-78 Mev gamma-ray. The total decay energy of 4:4 Mev supports 


b 


b 
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the assignment of this activity to 1°*Ir, for which the decay energy 
calculated from the semi-empirical mass tables of Metropolis and 
Reitwiesner (1950) is 3-3 Mev. 

(b) 7-minute }87Ir. The aluminium absorber curve found as the result 
of several runs indicated a maximum beta-particle energy of 1-6--0-2 Mev, 
the shape of the curve indicating one or more lower energy beta-transitions. 
Aluminium absorption of beta-particles in coincidence with gamma-rays 
showed that the higher energy beta-decay does not go to the ground state 
of the daughter nucleus. The gamma-ray spectrum was examined 
photographically and a gamma-ray of 1-8 Mev was found. 


Fig. 1 


%% Activity transmitted 


0 0-2 0-4 0-6 0-8 xe) \-2 |-4 16 
Aluminium g/cm? 


Absorption of radiation from 50sec 8Ir 


Any other gamma-ray energies were not resolvable from the large 
amount of gamma-rays from the 140-minute and 19-hour activities 
present. j 

(c) 9-7-day 1°°Ir, Aluminium absorption of the beta-particles from the 
mixture of this and the 70-day isotope indicated two beta-transitions of 
energy 80+20kev and 0-67 Mev. Beta-activity _ passing through 
12 mg/cm? of aluminium decayed with a 70-day half-life. The principal 
beta-energy of the 9-7-day isotope is therefore 80 kev, though there may be 


D2 
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up to about 5% of a higher energy beta-transition which ONG Se ae 
shown up in the decay through 12 mg/em? of aluminium. ‘ if 
absorption curve for the 9-7-day isotope was not obtained because the 
corrections for the 70-day background were large. 
Examination of the gamma-ray spectrum was effected te a scin- 
tillation spectrometer and single channel pulse analyser. Gamma-ray 
energies of 0-58 and 0-76 Mev were observed, the latter probably being 
two unresolved gamma-rays of similar energy. There were also some 
unresolved gamma-rays of about one Mev energy. The 0:58 Mev gamma- 
ray was resolved from a '**Ir gamma-ray of similar energy by following 
its decay relative to the 0-31 and 0-47 Mey '**Ir gamma-rays. 0-80 and 
0-94 Mev gamma-rays from 1%Ir were also observed in fair agreement 
with observations by Roulston and Pringle (1952) of 0-775 and 0-870 Mev. 

The ratio of beta- to gamma-counts and the beta- and gamma-counting 
efficiencies, as measured by coincidence counting techniques, indicated 
that there are relatively more beta-particles than gamma-rays and the 
decay scheme is therefore complex, some at least of the apparent beta- 
particles probably being conversion electrons. 

The assignment of the 9-7-day activity to 1°%Ir is plausible but not 
certain. 

(d) 140-minute 1°>°Ir. Absorption of beta-particles indicated a beta- 
transition of 2-1-.0-1 Mev energy together with a lower energy component. 
Absorption of beta-particles coincident with gamma-rays, however, 
indicated a beta-transition of energy 1-25--0-10 Mev, a softer component 
again being present (fig. 2). The 2-1 Mev beta-transition therefore goes to 
the ground state of the daughter nucleus. 

The activity of this isotope obtained for the investigation of the gamma- 
ray spectrum was rather low for examination with the pulse analyser, but 
gamma-rays of energy 0:42, 0-66 and 0-81 Mev were indicated together 
with a 70 kev component, presumably an x-ray, and other unresolved 
gamma-rays of low intensity above one Mev. The 0-66 Mev gamma-ray 
was resolved from a '4Ir gamma-ray of similar energy by following the 
decay. The 126 and 97 kev gamma-rays known to be associated with the 
Pt nucleus, following the decay of 195Au, were not observed. 
However, they would not follow directly after a beta-transition, as no high 
energy beta-disintegration followed by a low energy gamma-ray occurs 
and gamma-transitions from higher energy levels would not necessarily go 
via these low lying levels. 

The total decay energy of 2-1 Mev is considerably higher than the 
theoretical prediction of 1-0 mev based on the tables of Metropolis and 
Reitwiesner (1950). 


(e) 19-hour 1°4Ir, The gamma-rays from this isotope were examined 
using the pulse analyser scintillation spectrometer 
energy 0:32, 0:61, 1-18, 1-45 and ~1:8 mey were dete 
intensities being in the ratio 100: 22-10 ROL 
0-61 and 1:18 Mev gamma-rays 


and gamma-rays of 
cted, the approximate 
‘2 respectively. The 
have not previously been reported. 
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Sources for these experiments were obtained both by pile irradiation of 
‘Specpure ’ iridium and by cyclotron irradiations of platinum. 

(f) 7-second "Ir, This was prepared by the irradiation of platinum or 
iridium by fast neutrons in the cyclotron, the samples being transferred to 
the counter within ten seconds of the end of the irradiation. The gamma- 
rays of this isomer were examined with the photographic scintillation 
spectrometer and found to have energies of 125 kev and 5-6-+-0-2 ev, the 
latter being compared with the 4:45 Mev gamma-ray from a polonium- 
beryllium source, and measurements being made both on the peak due to 


Fig. 2 


A—Absorption of all 6-particles 


B—Absorption of 6-particles in 
coincidence with y-rays 


% Activity transmitted 


Rance ~=1000 mg Jem? 


) Qa) /0-4p aC GaemO-8 Ome awe G 
Aluminium g/cm? 


Absorption of radiation from 140 min 1°°Ir 


pair production in the Nal crystal and on that due to pair production with 
one annihilation quantum absorbed. The iridium K x-ray was also 
present in slightly greater intensity than the 125 kev gamma-ray. 

The half-life of the isomer was found to be 7:3--0-1 seconds by using a 
fast tape recorder and applying a discriminator bias to the counter so that 
only the 5-6 Mev gamma-rays were counted. The 125 kev gamma-ray 
and the x-ray were also found to decay with a 7-second half-life, using a 
differential bias discriminator. 
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‘Applying Weisskopf’s equations (1952) relating transition probability 
for gamma-ray emission with nuclear mass, spin change eo energy ea 
half-life of a 5:6 Mev isomeric transition in 1%Jr or 1%Ir (interna 


conversion being negligible at such high energy) will be as follows, using 
Goldhaber and Sunyar’s classification (1951) : 


Transition Type E M E M E M 
Spin change 5 5 6 6 7 7 
Parity change Yes No Now) #Yes Yes No ; 
Half-life (seconds) 610-8 10%) 2x 10-2 3°35) 92 1-6 x 10 


Similarly the half-life of a 125 kev transition (making no correction for 
internal conversion) will be : 


Transition Type M E M E 
Spin change 2 3 3 4 
Parity change Yes Yes No No 
Half-life (seconds) 2°3%10-* 0-32 56-2 2-7 x 108 


The K conversion ratio of the 125 key transition, as measured by the 
ratio of the intensities of the x- and gamma-rays in the spectrometer 
oscilloscope photographs and in pulse-analyser experiments, is about 1-5. 
From the calculations of Rose et al. (1951) this defines the 125 kev 
gamma-ray as an E2, K3 or M1 transition. The K/L conversion ratio for 
the E3 transition obtained from the graphs of Goldhaber and Sunyar 
(1951) is ~0-2. The calculated half-life of the 125 kev E3 transition, 
correcting the value obtained above for internal conversion, is therefore 
0-32 (1+1-5++5 x 1:5)=3-2 seconds. 

The decay from the metastable state may therefore occur in one of two 
ways : 

(a) A 7-second half-life, 5-6 Mev transition of type E6, M6 or E7 
followed by a 125 kev, partly converted, gamma-ray of type M1, E2 or 
E3. 

(b) A 7-second half-life, 125 kev transition of type E3 followed by a 
5-6 Mev gamma-ray of unknown type. 

The isomer cannot be assigned to !°2Ir which is known to have a meta- 
ee state of spin 3 above the ground state and it must therefore be 

94mT yp 

It is not found possible to detect this isomer after irradiation of iridium 
with thermal neutrons in the pile or with neutrons of energy less than 
about 15 Mev in the cyclotron. Since !®[r has a spin of 3/2 (Bene et al. 
1949) this indicates that the metastable "Tr has a fairly high spin, being 
formed only via a high spin compound nucleus resulting from capture of a 
neutron carrying high angular momentum. This is only possible with 
relatively high energy neutrons and this is why the cross section for the 
reaction increases with energy in spite of the increasing probability of a 
neutron being re-emitted instead of being captured. 
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4Tr is an odd-odd nucleus with 77 protons and 117 neutrons. Several 
shell model interpretations of the above decay schemes can be formulated. 
For example, shell theory suggests that the odd proton will be in the dy, 
level and the odd neutron in the pj). level. A ground level state of spin 1 
and odd parity can be obtained by the subtraction of these spins, and a 
first excited level of spin 2 and odd parity by their addition. A higher 
level of spin 8 and even parity from the addition of the ds, proton spin 
and an ij; neutron spin would then give a consistent interpretation of 
decay scheme (a), the 5-6 Mev rate-determining transition being M6 and 
the 125 kev transition M1. At such high energies, however, the shell 
model is probably not strictly applicable and multiple excitation of both a 
neutron and a proton may well occur. 

The occurrence of a high energy and high spin change transition as in 
decay scheme (a) would be very surprising in a heavy odd-odd nucleus 
where the many combinations of neutron and proton spins would be 
expected to give several levels of intermediate energy and spin, so that 
cascade emission could occur. In decay scheme (b) the energy, spin and 
parity of levels up to 5-6 Mev must be limited in such a way that the 
probability of transitions to them from the meta-stable level is less than the 
probability of the 125 kev transition, as well as in such a way that the 
5-6 Mev transition is preferred to a cascade. 
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ABSTRACT 


The mean lifetime of V,°-particles has been determined from measure- 
ments on 26 selected V-decays observed in a large cloud chamber. The 
chief criterion for selection was that the momentum of the positive 
secondary particle should be at least three times as great as that of the 
negative particle. It is shown that this selects events which are pre- 
dominantly V,°-decays. The decay scheme of the particles is assumed to 
be 

V,° > Pt+7-+~37 Mev. 

The result of the analysis gives a value of (3-7 17:3) x 10-10 sec for the 
mean lifetime. Combining this with other published values, it is concluded 
that the lifetime is probably between 3-0 and 4-5 x 10~” sec. 


§ 1. INTRODUCTION 


AxouT 200 neutral V-particle decays have been observed in an experiment 
carried out at the Jungfraujoch (3580 m). Half the decays were observed 
in a cloud chamber that has been described previously (Astbury e¢ al. 
1952). The remainder occurred in a new, deeper chamber with an 
illuminated volume measuring 55cm x 55cm x 16 cm situated in a magnetic 
field of 5200 gauss. The V-particles came from nuclear interactions in lead 
placed above the cloud chamber. 

An outstanding feature of the V-particle decays is their lack of sym- 
metry. In the majority of the events the momentum of the positively 
charged secondary particle is much greater than that of the negative 
particle ; events for which the converse is true are rare. This asymmetry 
is readily explained if a large proportion of the V-particles decay into a 
thai and a negative 7-meson. This mode of decay was suggested by 
pee cae i al. (1951) and further evidence for its correctness has come 
eee wee ees et al. (1953), Bridge et al. (1953) and Millar and 

; e decay scheme that is deduced for these particles is 


e VPs Ptta--~ 87 Mev, =. eo 
. oe asymmetry of their mode of decay has been used to separate 
1 decays from other neutral V-decays. From these events the mean 


lifet; i 

ae of the V 1"-particles has been estimated. The details of the 
paration and the analysis are discussed in § 2 and §3 below 

eater athena pea be EC seis : 
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§2. SELECTION oF V,°-PARTICLES 


The separation of V,°-particles is most easily described by reference to a 
study by Podolanski and Armenteros (1954) of the dynamics of two-body 
decay processes. They introduce a dimensionless parameter («) which is 
related to the masses (IM, m,, m_) and momenta (P, p,, p_) of the primary 
and secondary particles by the equations 

2 2 mm 2m 2 1/2 
a= Pt Pe + 2p" cos 6* (in B) sean (2) 
where p* is the momentum of either secondary particle in the centre of 
mass frame of reference and 6* is, in the same frame of reference, the angle 
of emission of the positive secondary particle measured from the direction 
of flight of the primary particle. 

Since the direction of decay in the centre of mass system is expected to 
be randomly orientated the average value of cos 6* for a number of decays 
will be zero. The average value of « for particles with the same mode of 
decay is thus 

m,2—m_? 
he 

For the V,°-particles «* has the value 0-69 while for any decay scheme in 
which m,—m_ the value of «* is zero. 

The range of the allowed values of « increases as P decreases. This is 
displayed in the figure where « is plotted against 1/P. The region within 
which « must lie for a particular mode of decay is bounded by the hyper- 
bolae (2) with cos 0*=-+1. In the figure these boundaries are shown for 
the V,°-decay scheme and for the suggested scheme V ,°->7+-+-7- + ~210 Mev 
(Armenteros et al. 1951, Thompson et al. 1953). 

Equation (2) reveals a property of « which is important in a statistical 
treatment of V-decays. If the decay process is isotropic in the centre of 
mass frame of reference the distribution of values of cos @* is rectangular 
between +1. This means that the «-values obtained from V-decays 
should be uniformly distributed between the limits shown in the figure. 
(The distribution is only uniform along lines of constant momentum ; 
along the 1/P axis it is determined by the spectrum of the decaying 
particles.) 

Most of the V-particles observed to decay in our experiment had 
momenta greater than 10° ev/c. The figure shows that « is greater than 
0-5 for nearly all V,°-particles with momenta in this range. Moreover, 
only about 25% of V,°-particles have «>0-5. 

Rather than use « itself, the ratio p,/p_ can be used to select decays 
with large values of «. If ¢ is the angle between the trajectories of the 
secondary particles in the cloud chamber 


P?=p,2+p_?+ 2p, p_ cos 6<(p,+p-_)*. 


ao 


From (2) 
1 Dk saa OTE sty 
oe ee eae eS —_—_$___— 
. Tee Isis 
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This shows that events for which p,>3p- recessarily satisfy the 
condition « >0°5. 

The three conditions satisfied by the events selected for ane wen 

(a) the measurable track of each secondary particle was at least cm 
long, 

(i the ratio p,/p— obtained by measuring 10 cm lengths of the tracks 
was greater than 3 and a 

the values of p,, p-, ¢ and the ionization densities of the tracks 
were consistent with the event being a V,°-decay. 


or] O-2 Z| O73 
/P [10° evic] 


Diagram showing the allowed values of « as a function of 1/P. The solid lines 
correspond to the decay scheme of the V,°-particles with a Q-value of 


37 Mev and the interrupted lines are for V,°-particles decaying into two 
m-mesons with a Q-value of 210 Mev. 


The points on the diagram represent the 26 events referred to in the 


text and the values of w have been calculated on the assumption that the 
events are V,°-decays. 


The upper limit of momentum measurement on tracks 10 em long is 
15x 10% ev/e¢ so that condition (b) selected, in general, events with 
p-<5xX 108 ev/c. However, for three events there was an evident origin 
for the V-particle and condition (b) was satisfied by angle measurements 
using the fact that p,/p sin 6_/sin 6, where 6 , and @_ are the angles 
the secondary tracks make with the line of flight of the primary. 


§3. MEASUREMENTS AND ANALYSIS 
For each of the 26 selected events the momenta p, and p_— were derived 
from the track curvatures. Frequently only a lower limit could be put 
to p,. However, by assuming the decay scheme (1) the values of Ps 
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and P were found from the more accurate measurements of p. and ¢. 
The possible error introduced by this procedure is discussed in § 4. 
The events were reprojected stereoscopically and the angle (4) between 
the secondary tracks was measured. The line of flight of the V-particle 
was determined and the distance (1) between the point of decay and the 
point where the V-particle entered the illuminated volume of the chamber 
was measured. In addition, the potential path (L) inside the chamber 


Table 1. Time of Flight Observations on V,°-Particles 


Event P l L t 7 
(108 ev/c) (cm) (cm) (10-19 sec) | (10-1 sec) 

JB 85 18-2 0-6 35:9 O-1 8-2 
LG 96 9-5 1-8 2:5 0-7 1:0 
NC 79 37-0 28:2 36-6 2:5 3°3 
NJ 97 9-5 2-6 6-1 1-0 2-4 
NR 14 30-4 15-5 32-0 1:8 3°8 
NU 178 31-0 1-2 32-2 0-1 3-9 
NW 35 15-8 11-0 18-5 2°5 4:2 
NY 21 58-5 26-0 28-0 1-6 ley 
NZ 49 27°5 215 32-0 2°8 4-2 
OC 64 28-1 16-7 37-5 2-2 4-9 
OD 555 19-5 10-0 33-2 1-9 6:3 
OE 119 18-2 2-0 8-4 0-4 13 
PB 150 9-7 13-6 37-0 5:2 14:1 
PB 206 47-5* 15-5 18-5 1-2 15 
PC) 151 25-2 10-5 34-0 1-6 51 
PE 346 16-0 16-4 23:8 3°8 5:5 
PEF 27 30-5* 30-0 35-5 3°6 4:3 
PH 449 33:8 2°2 35-6 0-2 3°9 
PK 28 20-0 2-6 35:6 0-5 6:8 
QI 575 12-0 4-5 33-0 1-4 10-4 
QJ 118 29-2 25:6 39-7 3-3 5-2 
QJ 176 28-8 12-0 37-0 1-6 4:8 
QL 46 6-7 2-2 6-4 Meg 3°5 
QL 272 16-5 10-0 15:3 2:2 3-4 
QL 329 8-9 0:2 3-0 0-1 1-3 
QL 759 39-0 1-8 35-0 0-1 3°2 


Details of the 26 events analysed in the experiment. P is the momentum 
of the neutral V-particle, / and L are the length of the trajectory in the chamber 
before decay and the potential length respectively. ¢ and T' are the corres- 
ponding proper times of flight. 


was found. This potential path is the distance between the extreme 
positions on the V-particle’s trajectory where the decay would have 
satisfied the conditions given in §2. This meant that every point where 
the decay would have produced secondary tracks 10 em long was 


included in the potential path. 
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The times of flight (¢, 7’) corresponding to the distances / and L were 
then obtained from the relation t=lx M/P. These times are given 1n 
table 1 together with the values of /, L and P for the 26 events. 

Bartlett (1953 a) has described in detail how a maximum likelihood 
estimate of the mean lifetime (+) can be obtained from a set of values 
(t;, 7’;). He shows that the statistical error in the estimate of the decay 
constant (1/7) is approximately symmetrical and gives an interpolation 
procedure for finding the best value and standard error of 1 [T. 

This procedure, when applied to the times of flight listed in table 1, 
leads to the result that 


1/7=(0-27-0-14) x 101° sec? 
corresponding to 7= (3-743) x 10-10 sec 


where only the statistical error is quoted. Other sources of error are 
discussed in § 4. 

From the measured values of p— and the derived values of p, and P 
the values of « for the decays were calculated. The results are displayed 
in the figure. The fact that all the values of « are greater than 0-6 
means only that the selection of events with p,/p_>3 was over-cautious. 
If the events are V,°-decays the ratio p,/p_ is greater than 4 in every 
case. It must be emphasized that the points in the figure do not represent 
measured values of «; they are calculated on the assumption that the 
events were V’,°-decays. 


§4. Non-SratisticaAL ERRORS 


There are two important sources of error in the above result that 
have to be considered. The first is the presence among the selected 
events of decays not due to V,°-particles and the second is the uncertainty 
in determining the momenta P from the measured values of p_ and ¢. 
By comparison, the errors involved in measuring / and L are trivial. 
The energy release in the decay of V,°-particles may well be greater than 
the value of 37 Mev that we have assumed (Armenteros et al. 1953 give 
a value of 40 to 45 Mev) but this uncertainty has only a small effect on 
the derived lifetime. 

Tke number of particles erroneously classified as V ,°-particles can be 
estimated if the «-values for all other V-decays are symmetrically 
distributed about «=0. This is true of the decay into two 7-mesons 
suggested for V,°-particles and would also hold for a decay into a heavy 
and a light meson (Leighton et al. 1953) unless the heavy meson were 
preferentially positively charged. 

Only three of the 200 V-decays examined have «<—0-5 and would 
pute ee selection conditions of §2 if the signs of the decay products 
were reversed. From the symmetry 


ee pes assumed for the «-values we expect 
events which are FiO hector : . 
ch are not V,°-decays to have been included among the 
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26 events analysed.t| The bias due to this dilution of the V,°-decays is 
certainly small—especially since the lifetime estimated for the V ,°-particles 
is of the same order of size as that of the V,°-particles (Fretter et al. 1953, 
Bridge et al. 1953). 

An error arises in the determination of P from p_ and ¢ since, for 
given values of p and 4, P is double-valued. In many cases one of 
the possible values can be excluded by measurement of p, but for some 
events values of P differing by a factor of two would fit the measurements 
equally well. In these cases (marked by an asterisk in table 1) the 
higher of the two values of P was used in the analysis. Fortunately, 
the events for which this uncertainty arises contribute very little to the 
final value of 7. 

Apart from this ambiguity there is a random error of ~25°% in the 
values of P arising directly from the errors in p_ and ¢. The effect of 
this error is much smaller than the present statistical error in rv. 


§ 5. Discussion 


The result quoted in §3 can be compared with several other reported 
values for the mean lifetime. These are set out in table 2. 


Table 2. Estimates of the Mean Lifetime of V,°-Particles 


At, Number Mean Lifetime 
ee of events (10-48 26c) 
Fretter e¢ al. (1953) 22 10+7 
Alford and Leighton (1953) 74 2-5+0:7 
a : (‘ High Q’ events) 20 1340-5 
x af (‘ Low Q ’ events) 37 2:9+0°8 
Bridge et al. (1953) 21 3541-2 
Deutschmann (1953, unpublished) 22 4-87 26 
Gayther (1953, unpublished) 21 4-073 
Present paper 26 Ouielee 


The values of the mean lifetime of V,°-particles obtained in six different 
experiments. The results of Alford and Leighton and that reported in this 
paper were obtained from V-decays observed in magnet cloud chambers. 
The other four results have been obtained from multiplate cloud chamber 

experiments. 
The various results are discussed briefly in § 5. 


+ Newth and James (1953) concluded that only about 50% of the fast 
V-decays observed in these experiments were )’,°-decays. ie the remainder 
are V,0-decays inspection of figure shows that a quarter of them (ees 129, OF 
all V-decays) will have a>0°5. This leads to the conclusion that 80 % of the 
decays with a>0-5 are V,°-decays. The proportion of V,°-decays is thought 
to be higher among the 26 events reported here partly because the value of a 
was not the only criterion of selection and partly because Newth and James 
considered V-particles with momenta above 10° ev/c only—for lower momenta 
the proportion of V,°-particles may well be higher than 50%. 
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The result of Fretter and his collaborators may be biased since ee) 
disregarded the presence of the lead plates in their cloud chamber w * 
analysing their measurements. They state that if the plates et he J 
spaced ignoring their presence introduces no bias in the resu . his 
is certainly true if the plates are sufficiently thin. However, if the aa 
spent by a V-particle in such an unobservable region of the cloud cham er 
is comparable with its mean life this fact must be taken into account in 
the analysis. If it is ignored a bias tending to increase the derived 
lifetime is introduced. How important this bias is in any particular 
experiment is difficult to assess but Bridge et al. (1953) have shown how 
the bias may be avoided and Bartlett (1953 b) has extended his estimation 
equations to deal with this kind of experiment. The published data of 
Fretter et al. are not sufficiently detailed for an unbiased estimate to be 
made from them. In addition, since the statistical error in the result 
depends upon the potential times of flight in a rather complicated way, 
ignoring the lead plates which obscure about 20% of the trajectories 
of their particles has led Fretter et al. to underestimate the error in their 
result. 

Alford and Leighton have analysed by far the largest number of 
V-decays and have estimated lifetimes from the maximum likelihood 
equation. They use an approximate formula for the statistical error 
which gives a result symmetrical in +. The classification of their 
V-particles is based on the work of Leighton et al. (1953) who state that 
only 20% or less of their V-decays could be of the V,° type where the 
secondaries are of equal mass. Nearly all the remainder are consistent 
with decay into a proton and a 7~-meson with a measured energy release 
varying from 10 to 100 Mev. This wide spread in Q-values leads Alford 
and Leighton to make tentative lifetime estimates for the ‘low Q’ 
events and the ‘ high Q’ events, the division being made at a Q-value 
of 50 Mev. The results of these separate estimates are given in table 2. 

There is no difficulty in identifying the ‘low Q’ decays with the 
V,°-decay scheme (1) but Leighton et al. consider it ‘quite unlikely ’ 
that the whole spread in @-values is due solely to errors in measurement. 
It seems possible to us that some of the ‘ high Q’ events are, in fact, 
V,°-decays. This interpretation would reduce the discrepancy between 
the small number of V,°-decays found by Leighton et al. and the larger 
proportion found by Newth and James (1953) and by Thompson e¢ al. 
(1953) in similar experiments. It would also account for the shorter 
lifetime of the ‘ high Q’ events since the lifetime of the V ,°-particles is 
shorter than that of the V,°-particles (Fretter et al. 1953, Bridge et al. 
1953, Astbury 1953). However, the ‘high Q’ events with identified 
proton secondaries could not be explained in this way and the scarcity 
of decays with negative «-values found by Leighton et al. would have 
to be regarded as an unlikely statistical fluctuation. It is clear that 


further experimental evidence is necessary before the 


se difficulties can 
be resolved. 
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Bridge et al., Deutschmann, and Gayther observed their V-decays in 
multiplate cloud chambers. In each experiment the measurements have 
been analysed to avoid the bias referred to in connection with the result 


of Fretter et al. Bridge et al. consider that their result can best be 
expressed in the form 


1/7=(0-29+-0-11) x 10° sec-1 
or T=(3°5 75) X10 °° sec 
to take account of the asymmetry of the error in r. 

These workers also state that they have made no allowance for the 
nuclear interaction of the V-particles in the lead plates of their cloud 
chamber. If the cross section for nuclear absorption were appreciable 
their estimate of the lifetime would be biased to a low value. A very 
rough calculation of this effect using their data and assuming an absorption 
length of 150 g/cm? for V-particles in lead gives a bias of less than 5%. 
It is clear that the determination of the nuclear interaction of V-particles 


from lifetime measurements made in chambers with and without lead 
plates is a very remote possibility. 


§ 6. CONCLUSIONS 


The values for the mean lifetime given in table 2 can be combined 
together to give a weighted mean value. This has been done by weighting 
the separate values of 1/7 inversely as their variance and then finding 
the weighted mean value and standard error in 1/7. This procedure is 
justifiable since the errors in the separate results are largely statistical 
and the values of 1/7 are expected to be approximately normally 
distributed. 

The result of Fretter et al. has been excluded since it may contain a 
large bias. The result obtained by Alford and Leighton for their ‘ low Q ’ 
events has been used in view of the admitted uncertainty about the 
homogeneity of the events included in their main analysis. The fractional 
error in + quoted for their result by Alford and Leighton has been 
transferred directly to the 1/7 scale. 

With these restrictions the weighted mean lifetime is found to be 
given by 

1/r=(0-27+0-05) x 102° see“? 
or T==(3' Te) 6 10 BEC. 


The agreement between the results obtained by workers with different 
types of apparatus observing V-particles of different energy is extremely 
satisfactory. The error in the ‘ best value’ given above is still mainly 
the inherent statistical error. Since this value is derived from measure- 
ments on 127 events it is clear that many hundreds of events will be 
required before the statistical uncertainty in the lifetime can be made as 
small as the error due to our ignorance of the exact Q-value of the 
V,°-particles. 


46 On the Mean Lifetime of the V,°-Particles 


ACKNOWLEDGMENTS 


The events discussed in this paper were recorded at the J ungfraujoch 
Research Station in Switzerland. We are grateful to the Administration 
of the Station and to the manager, Mr. Hans Wiederkehr, for making 


our work there possible. 

Several colleagues in Manchester have cooperated in the maintenance 
of the cloud chamber and in interpreting the results. In particular, 
Dr. J. P. Astbury made many calculations and measurements that we 


have used. 
Dr. M. Deutschmann and Mr. D. B. Gayther have very kindly allowed 


us to quote their conclusions before publication. Dr. R. Safford of the 
Massachusetts Institute of Technology and Professor M. 8. Bartlett in 
Manchester have been most helpful in discussing the results of experiments 


with multiplate cloud chambers. 
We have profited greatly from the advice of Professor P. M. 8. Blackett 


and are grateful to him for supporting our work. One of us (D. I. P.) 
has received a grant from the Department of Scientific and Industrial 


Research during the course of this work. 


REFERENCES 


AuLForD, W. L., and Letcuton, R. B., 1953, Phys. Rev., 90, 622. 

ARMENTEROS, R., Barker, K. H., Burzer, C. C., and Cacuon, A., 1951, Phil. 
Mag., 42, 1113. 

ARMENTEROS, R., BarKer, K. H., Burumr, C. C., Coates, M. S8., and SowERBY, 
M. G., 1953, Phil. Mag., 44, 861. 

Astpury, J. P., 1953, Report of the Bagnéres Conference on Cosmic Radiation. 

Astsury, J. P., CaIpprInpALE, P., Mizuar, D. D., Newru, J. A., Page, D. I, 
Rytz, A., and Saurar, A. B., 1952, Phil. Mag., 43, 1283. 

BarteTt, M.8., 1953 a, Phil. Mag., 44, 249 ; 1953 b, Ibid., 44, 1407. 

Dare Ht S., Peyrou, C., Rossr, B., and Sarrorp, R., 1953, Phys. Rev., 91, 

62. 
Fretter, W. B., May, M. M., and Naxapa, M. P., 1953, Phys. Rev., 89, 168. 
ey R. B., Wanuass, 8. D., and AnpERson, O. D., 1953, Phys. Rev., 89, 
48, 

Miciar, D. D., and Paas, D. I., 1953, Phil. Mag., 44, 1049. 

Newta, J. A., and Jams, G. D., 1953, Report of the Bagnéres Conference on 
Cosmic Radiation. 

PopoLanskt, J., and ARMENTEROS, R., 1954, Phil. Mag., 45, 13. 

Tsompson, R. W., Buskirk, A. V., Errer, L. R., Karzmarx, C. J., and 
Repiker, R. H., 1953, Phys. Rev., 90, 329. 


a 


bat] 


VI. Angular and Energy Distributions in Energetic Pion Interactions 
with Nuclei 


By A. H. Morrisu 
Radiation Laboratory, McGill University, Montreal, Canada* 


[Received September 7, 1953] 


ABSTRACT 


A total of 405 events due to the interaction of 210-+-20 mev 7--mesons 
with the nuclei of photographic emulsion are analysed. Angular and 
energy distributions of the charged reaction products are interpreted as 
giving evidence of a cascade process for those events in which the meson 
survives. The characteristics of the remaining events, in which the meson 
is either absorbed or undergoes charge exchange, are primarily determined 
by the evaporation process. 


§ 1. INTRODUCTION 


PHOTOGRAPHIC emulsions, exposed to a beam of 210+-20 Mev z~-mesons, 
have been searched for events which have either at least one grey track 
produced by a proton (200>H>30 Mev), or at least one black prong 
(H#<30 Mev). The large number of events found now permits reliable 
studies of angular and energy distributions of the charged reaction 
products to be made. Such studies can, in principle, give evidence for 
or against the existence of a cascade process inside the target nucleus. 
There have been several papers on evidence for cascades produced by high 
energy protons in nuclei. (Bernardini et al. 1952 a, b, Randle et al. 1952, 
Hofman and Strauch 1953, Nelson et al. 1952, Fishman and Perry 1952.) 
However, this is the first detailed experiment using energetic mesons. 
Goldberger (1948) has developed a model for the behaviour of a fast 
nucleon inside a nucleus. He represents a heavy nucleus as a Fermi 
gas of nucleons, and considers the incoming energetic nucleon to interact 
with this gas by a series of single nucleon-nucleon scatterings. The 
recoil nucleons in such a process, are assumed to interact in a like manner 
with other nucleons in the gas, and in this way the cascade develops. 
Some of the fast nucleons escape from the nucleus almost immediately. 
Finally, the remaining excited nucleus ejects some slow particles in 
accordance with the evaporation theory. (Weisskopf 1937.) Goldberger 
retained the statistical fluctuations which occur for each interaction in 
the cascade by using the Monte Carlo method (Ulam and Von Neumann 


* Communicated by the Author, 
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1947), when carrying out calculations for 90 Mev neutrons interacting 
with Pb nuclei. Bernardini et al. (1952 b) subsequently made calculations 
for 300-400 Mev protons interacting with nuclei with A~100. This 
value of A was chosen since it had been found that 80% of the stars 
produced by 300-400 Mev protons in emulsions were due to collisions 
with Ag and Br nuclei. Experimental results with emulsions supported 
the theory. ; 

In the above calculations, the experimentally determined free nucleon— 
nucleon cross sections were employed. A similar calculation for meson 
interactions with nuclei could be made if the appropriate meson cross 
sections are known. Differential cross sections for elastic 7~-meson—proton 
collisions and charge exchange interactions in hydrogen have been 
determined at 120 and 144 mev by Anderson et al. (1953), and further 
work at 227 Mev is in progress now (Anderson 1953) at Chicago. There is 
almost no data available regarding the absorption of a fast 7~-meson in a 
nucleus. While Goldberger type calculations of meson-nucleonic cascades 
in nuclei cannot be made at the time of writing, they may be possible in 
the near future. In this work the conclusions drawn from the data are 
essentially qualitative. 


§ 2. ExPERIMENTAL METHOD 


Electron sensitive Ilford G5 emulsions, 400 microns thick, were exposed 
to the 227 Mev 7~-meson beam from the University of Chicago synchro- 
cyclotron, and suitably processed. Considerations of energy loss by 
ionization in the stack of plates and measurements using the multiple 
scattering technique led to the conclusion that the average pion energy 
was about 210+20 Mev. Further, a scan for tracks of about minimum 
ionization perpendicular to the pion beam direction showed that less than 
2%, of the tracks were due to energetic cosmic ray particles. About 5% 
of the tracks were produced by .-mesons which resulted from 7—u decays 
in flight. 

Previously, values of cross sections for various types of interactions 
of 210 Mev z~-mesons with the nuclei of photographic emulsions were 
determined in this laboratory (Morrish 1953 a), by scanning along indi- 
vidual tracks. Because this method is arduous and tedious, only 2823 em 
of pion track were followed. -Since that time, the number of events has 
been greatly increased by ‘ volume’ or ‘ area’ scanning using a high power 
(x53) oil objective. However, even with this great magnification; 
certain types of events are missed. These are ‘stop’ events, in which 
the meson disappears in flight with no visible disintegration products 
and ‘ elastic scattering ’ events, in which the meson is scattered through 
an angle greater than 5° with a grain density change of less than 4%, 
habe aes Henin by any visible disintegration products. There 

idence that a few events of other types are missed ; but in 
general this is not serious. This loss will be discussed in the next section. 
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The following classification of the remaining events was adopted : 

Type 1. Stars. Events in which the 7~-meson disappears and produces 
a nuclear disintegration with one or more black prongs. (E<30 Mev.) 

Type 2. Inelastic scatterings. Events in which the meson survives,* and 
is also accompanied by one or more black prongs. (£<30 Mey.) Such 
events with no black prongs were missed in the ‘ volume ’ scan. 

Type 3. Stars with one or two fast protons. Stars accompanied, in 
addition, by one or two grey tracks produced by protons. (H>30 Mev.) 

Type 4. Uncertain events. Hither inelastic scatterings or stars accom- 
panied by a fast proton. 

Type 5. Inelastic scatterings with a fast proton. Inelastic scatterings, 
accompanied by a grey track produced by a proton. (H>30 Mev.) 

Black prongs were usually identified visually. If of sufficient length, 
tracks of fast protons or of inelastically scattered mesons were identified 
by multiple scattering and grain density measurements. 

Identification of short grey tracks as pions was made if (a) the grain 
density indicated the particle, if a proton, had an energy higher than the 
total energy of the incoming meson, or (b) the grain density indicated 
the particle, if a proton, had an energy greater than the kinetic energy, 
but less than the total energy of the incoming meson. The latter is 
believed to be an improbable process in negative meson absorption. 
Unambiguous identification was not possible for other short grey tracks ; 
such events are listed as uncertain. 

Also, it should be noted that determinations of the energy of the pions, 
on the basis of grain density measurements alone, are subject to large 
errors at energies above 100 Mev. This is because of the small change 
in grain density in this region (between 117 and 222 Mev the change is 
only about 6-5°% ; Morrish 1953 b) and because of the relatively large 
statistical fluctuations which exist for short tracks. 

It is very difficult to discriminate between slow protons, deuterons, 
tritons and alpha particles in electron sensitive emulsions. It appears 
that most of the black prongs were due to protons, but no detailed 
classification of origin of the black prongs has been attempted. 


§ 3. Man FREE Paros FOR INTERACTIONS 


A total of 405 events of the type listed in the previous section were 
found. Of these, 74 had been observed in the previous ‘ along the track ’ 
scan, the remainder were found in the ‘ volume’ scan. Table | records 
the frequency of events found as a function of the number of black 
prongs. 

The relative frequency distributions for stars and inelastic scatterings 
for the two methods of scanning were compared, and, within the statistical 
limits, were found to be the same. In addition, the relative numbers 
of events of types 1, 2, 3 and 4 were the same in both scans. However, 


*This assumes that the number of ‘double charge exchange’ events 
(1n->7r°->7+) is small. 
EZ 


50 A. H. Morrish on the Angular and Energy 


there was a definite loss of type 5 events (about 75-80%) mae 
by no black prongs in the volume scan. This is nt Boe ; 
since usually all tracks in such events have fairly low grain densi am a 

The total track length in the volume scan was determined by oy e fe 
the track length in numerous small sample volumes. — ee % i he 
mean free path for the type 1—5 interactions was determine : to be ‘ ee 
The corresponding ‘along the track’ scan value was 38-2 cm. 


Table 1. Frequency Distribution of Events 


Number of black prongs 0 ||. 2.137) 4 Sie Ge ee ee oat 


Stars — |.28 | 57 | 63 |. 55°} 30°) TT] vogweese 
Inelastic scatterings shi) 40 \) 23°) TEST 2 ee 89 
Stars with one fast proton Fo That. Gil O01 Leah eee io 
Stars with two fast A eae OF Oe] 2599 05" Or) SOs a ater Ree 
Uncertain events : stars with 
ee para ge tik LU tele 62 eee 28 
Inelastic scatterings with a f 
fast proton TA Bt 0:1) CE 0h a ae Il 
Total ‘ an Maar eo 405 


type 5 events with zero black prongs are excluded, the m.f.p.’s are 47-5 
and 40:4 cm respectively. In other words, there is a 17° difference. 
Probable errors are about 7° and 8°% respectively, where errors in the 
determination of the total track length for the volume scan are included. 
There is, therefore, some indication that a small percentage of other 
events have been missed, probably primarily those with only one black 
prong, and then most often when the prong is roughly perpendicular to 
the emulsion plane. The considerations of the previous paragraph 
however, seem to indicate that the total loss of events of type 1—4 is less 
than 10%. One has considerable confidence that such a small loss will 
not appreciably affect the angular and energy distribution data. 
Previously, comparison with results obtained with 50-100 mev 
m-mesons at Columbia University showed that the number of black 
prongs per event increased at our energy, as expected (Morrish 1953 a). 
From table 1, it is seen that the stars have more black prongs per event 
than the inelastic scatterings. The stars represent interactions where 
either the incoming meson has been absorbed or has undergone charge 
exchange, and emerged as a neutral meson. It is clear therefore, that 
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the increase in the number of black prongs in the stars is primarily due 
to the additional meson rest mass energy available in an absorption. 
Indeed, the distributions suggest that the absorption process may 
predominate over the charge exchange process. If it is assumed that, 
for the 7 prong stars, all the energy available is transferred as excitation 
energy, it follows that the average energy per black prong is about 
50 Mev. Bernardini et al. (1952 a) obtained a similar figure in their work 
with protons. 

The fact that there were 89 inelastic scatterings found as compared 
to 28 stars accompanied by a fast proton, suggests that the majority 
of the uncertain events were actually inelastic scatterings. Arguments 
concerning the frequency distributions are not conclusive, but comparison 
of the angular distributions strongly supports the suggestion. We will 
not discuss the uncertain events further in this work. 


Table 2. Geometrical Mean Free Paths for Emulsion 


Element Ag Br I C H O S N 


PAUS\ oye 1 
Nah? | — ; 2 

pec 71-5 | .86°9 | 5820 | 236 498 261 | 8850 | 1010 
(cm) 


m.f.p. heavy nuclei : 39-2 cm m.f.p. light nuclei: 89-4 cm 


m.f.p. light nuclei ex- 
cluding hydrogen : 109 cm 


m.f.p. all nuclei : 27-25 cm 


m.f.p. all nuclei, excluding hydrogen : 28:8 cm 


Table 2 records the ‘ geometrical ’ mean free paths for the constituents 
of G5 emulsions assuming the nuclear radius to be equal to hA'?/uc, 
where p is the mass of the 7~-meson. If the cross sections for 210 Mev 
a--meson interactions with emulsion nuclei are roughly proportional to 
the geometrical one, then it follows that about 70% of our events occur 
in Ag, Br, or I nuclei. Further, elastic pion—proton events can be 
distinguished. If we neglect the hydrogen of the emulsion, this figure 
is raised to about 74%. It is reasonable to assume that the great majority 
of the interactions are with the heavier nuclei. 

Tt would be of value, however, to determine the precise contribution 
that the light elements play on the results reported in this paper. There- 
fore, further experiments along the lines of this investigation using either 
sandwich or diluted emulsions would be worthwhile. 
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§ 4, ANGULAR AND ENERGY DISTRIBUTIONS 


The angle of the prongs, relative to the incident pion direction, was 
calculated from measurements of the prongs’ projected and ‘ dip ° angles, 
taking into account the shrinkage factor of the emulsion. The results 
for the black prongs occurring in the stars, and in the stars accompanied 
by a fast proton, are shown by the histogram of fig. 1. The results are 
considered per equal solid angles by plotting them in equal cosine intervals. 
Separate histograms were also plotted for stars with the same black prong 
size, and also for those accompanied by a fast proton. Each, within 
the statistics available, showed the same general shape as the one of fig. 1, 


Fig. 1 


FREQUENCY 


'O 86 e§ 4 2 QO -2 --§@ --6 -—-g -r'o 


COSINE 6 
Histogram of the angles of the black prongs (#<30 Mev) in stars 


Figure 2 illustrates Ke 
Beis fe te: ates the black prong angular distribution for the inelastic 
§ a Ai ‘ a a . T 
ee ngs. Also, separate histograms for events with equal numbers of 
ack prongs were very similar in shape to fig. 2 
iy 


; oe : 2. It is immediatel 
obvious that there is a considerable difference between the distributions 


in figs. 1 and 2 icts 
ee ae s ea predicts that the particles emitted during the 
es Bees are isotropic in direction. The work with fast 
oe Si brata 
» previously referred to, showed a preferential emission of slow 
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Angles of black prongs, relative to the incident pion direction, for inelastic 
scattering events. 


Fig. 3 
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Fast proton angular distribution from (a) stars, and () inelastic scatterings. 


protons in the forward direction. This is predicted by the Goldberger 
calculations, and is due to a ‘knock on’ process resulting from the 
nucleon cascade. It seems likely, therefore, that the distribution of 
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fig. 2 is mainly due to a cascade process, generated by the incoming 
meson. The distribution of fig. 1 is much more isotropic in nature, but 
even here there is a tendency for an increase in the number of particles 
ejected forwards. It is tempting to consider this distribution to be made 
up of two parts ; one part, a cascade type of distribution which may be 
due to either the charge exchange events or to those absorption events in 
which one or two fast nucleons are produced which initiate a cascade. 
The other part, which appears to predominate, may be due to those 


Fig. 4 


FREQUENCY 


COSINE 6 


Angular distribution of pions emerging after a nuclear interaction. 


interactions in which the meson is absorbed by a group of nucleons, with 
the result that the nucleus is left with a certain amount of excitation 
energy. The amount of this energy depends upon the number and energy 
of the fast neutrons which escape almost immediately, apparently without 
initiating a cascade. The presence of one or more fast neutrons coming 
from the majority of these stars is inferred by the conservation of energy 
principle. 

It is unfortunate that in this type of experiment these neutrons cannot 
be detected and analysed. Instead, we will now consider those few 
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events which are accompanied by a fast proton. Histograms of the 
angular distribution of fast protons associated with (a) stars, and 
(b) inelastic scatterings, are shown in fig. 3. It is interesting to note that 
none of the fast protons of the inelastic scattering events emerge in the 
backward direction in the laboratory system, whereas, for the stars, more 
than 25° emerge in the backward direction. This appears to give 
evidence that the former events may be due to a quasi-single pion—proton 
elastic scattering inside the nucleus. The argument is supported also by 
the observation, as shown in table 3, that these pions usually emerge 
in the backward direction. The fast protons ejected from stars in he 
backward direction could again be explained as meson absorption by a 


Table 4. Energy Distribution of Emergent Pions and Fast Protons 


Hnergy Mev 0-30|30-60|60—90]90—1 20] 120-150]150—180|180—210|over 210 


31 26 7 7 3 8 2 


Or 


Pions from 
inelastic 
scattering 
events 


Pions from 
inelastic 
scattering plus 
a fast proton 
event 


—_ 
— 
e 
_ 
—_ 
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> 


Fast protons — ayy eal 10 3 0 1 0 
from stars 


Fast protons — 1 4 6 0 ' 
; 0 

from inelastic : Y 
scatterings 


as or ist of nucleons in a nucleus. The limited statistics 
ilable here show the desirabilit <peri 
Sag en oe bility of extending this type of experiment 
Figure 4 is a histogram of the angular distribution of the ejected pions 
of the inelastic scattering events. There is some indication of a minimum 
at about 60°, and a slight increase at about 160°. A little more than half 
(49 of the 89 pions) are emitted at angles greater than 90°. However 
when ie angular distribution is considered for different energy eons 
ecease wanes apparent. As shown by the results in table 3, 
nihen ieee ae Betas to be ejected backwards, while those of 
aide : By end to be in the forward direction in the laboratory system. 
recor ds the energy distribution of the emerging 7~-mesons ; th 
great majority have an energy less than 90 Mev . pai 


aft) . 
energy a 215 Mev 7- Since the maximum 


-mes 18 1N ¢« 1 j 
eson can lose in an elastic single nucleon collision 
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is about 135 Mev, it is clear that these results indicate these mesons 
usually undergo more than one nucleon collision inside the Ag or Br 
nucleus. 

Finally, the energy distribution of the fast protons from stars and 
inelastic scatterings is tabulated in table 4. 


§ 5. CONCLUSIONS 


The observations have shown a striking difference between the angular 
distributions of the black prongs from stars and inelastic scatterings. 
The latter distribution appears to imply that cascade generation in the 
nucleus determines the main characteristics of the star, rather than the 
evaporation process. On the other hand, the black prong distribution 
for the stars appears to be primarily due to evaporation. 

On the basis of earlier cosmic ray work, it has been the fashion to 
conclude that such black prongs are almost entirely due to the evaporation 
process. The present results, together with those of Bernardini e¢ al. 
(1952 b), now provide ample evidence that this conclusion is not true in 
general. 

The emulsions were exposed to the University of Chicago synchrocyclo- 

tron through the kind cooperation of H. L. Anderson. 
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VII. Gamma Radiation from the Decay of *5>Th 
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ABSTRACT 


The energies and relative intensities of the y-rays following the «-particle 
decay of 228Th have been investigated with proportional counters. Lines 
were observed at 84:4+0-2, 132:3-1-3, 167+2 and 214+3 kev. The 
multipolarity of some of the y-ray transitions is discussed, and two possible 
decay schemes are presented. 


THE gamma radiation which follows the «-particle decay of radiothorium 
has been investigated in the range 30-220 kev using krypton and xenon 
filled proportional counters. When this work was begun, the only gamma 
radiation reported from this isotope was a line at 84-3 kev (Rosenblum 
et al. 1952) which had previously been interpreted as a doublet (Surugue 
and Tsien 1941). The resolution obtainable with the proportional 
counters would allow confirmation of this result. 

The #8Th was obtained by extraction from a pure (<10~¢ MsTh) 
solution of radiothorium chloride in equilibrium with its decay products. 
Two different samples of radiothorium solution were obtained from the 
Radio Chemical Centre, Amersham, and two different separation 
procedures were adopted during the course of the work; these gave 
similar results. 

The first. separation involved the precipitation of the daughters as 
sulphides and sulphates (Beling 1951). In the second separation the 
228'Th as nitrate was extracted by washing with ethyl acetate (Bock 1950). 

The energy measurements were made using the 59-78 kev y-ray from 
the decay of #41Am (Browne 1952) as a standard, other work having shown 
that the pulse height output of the detection system was linear with y-ray 
energy to +0:5% up to 200kev. The spectrum was measured at 
intervals up to 24 hours after the separation of each source in order to 
decide whether an observed line was due to ?8Th or to one of its decay 
products. The principal y-activity from the decay products is due to 
ThB and subsequent members of the chain. A gamma ray from an 
early member of the chain will build up with the half-life of ThX, namely 
3 days, whilst activities due to ThB or later members will remain low for 
the first 10 hours (the half-life of ThB) and then increase with the half-life 
of ThX. A gamma ray from 28Th (T4=1-9 years) will of course remain 


* Communicated by the Authors. 
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constant during this period of measurement. Typical curves taken 
with a xenon filled counter at 1 hour and 7 hours after the source separa- 
tion are shown in fig. 1 (A and B). 

In this way it was possible to ascribe to the decay of 228Th lines at 
84:4+-0-2 kev, 132:3-41-3 kev, 167+2 kev and 214+3 kev. The energy 
of 84-4 kev is to be compared with that of 84-3 kev obtained by Rosenblum, 
Valladares and Guillot (1952). The line width we obtain supports their 
conclusions that this is a single y-ray and not a double line as had been 
proposed by earlier work (Surugue and Tsien 1941). The gamma 
radiation from the decay of ??°Th has also been investigated by Boussiéres 
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Curves taken with a xenon counter (A) | hour and (B) 7 hours after separation. 
The peak at 102 kev is the escape peak of the 132 kev line. The line 
at 75 kev in (B) is attributed to bismuth K radiation from ThB decay. 


et al. (1953) using a sodium iodide crystal. Their interpretation of their 
curves is in fair agreement with our work. 

The counters used had an internal diameter of 5 cm and were filled to 
a gas pressure of ~60 cm. The proportion of photoelectrons produced 
in the counter gas which lose all their energy in it becomes small as the 
range of the photoelectron approaches the radius of the counter. It is 
not therefore possible to infer the relative intensities of gamma rays 
directly from a knowledge of the intensities of their photo peaks and 
corresponding mass absorption coefficients in the counter gas when this 
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situation occurs, i.e. for gamma rays of energy greater than ~60 kev 
in our case. Correction curves for this effect published by West (1952) 
do not cover the energy range in which we were concerned. We were 
able to show (Newton and Rose 1953) that it was possible to obtain the 
necessary correction factor by observation of the change of counter 
efficiency with applied longitudinal magnetic field. The electron tracks 
are curled up by the field, thus increasing the number which lose all their 
energy in the counter gas. In the case of the krypton filled counter, a 
longitudinal field of ~2700 gauss was used in order to get peaks for the 
167 kev and 214 kev radiation. 

The relative intensities of the lines which we attribute to the decay of 
228Th are shown in the table. 


Energies and Relative Intensities of the Gamma Rays from 
the Decay of Radiothorium 


Agta Energy in | 94.440-2 | 13234153 167+2 21443 
ev 
8 


Relative Intensity 1 0:105+0-01 | 0:058+0-012 | 0-18+0-05 


We were also able to set an upper limit of 5 x 10~? to the relative intensity 
of ThX K-radiation, which is sufficient to rule out the possibility that any 
of the three higher energy transitions is predominantly magnetic. 

The «-fine structure of ?28Th as measured by Rosenblum, Valadares 
and Perey (1949) shows groups «) and «, corresponding to transitions to 
the ground and first excited states of ThX respectively, and also a group 
which he interpreted as being due to the «, group from the decay of ThX 
(Rosenblum et al. 1949) which was present as a contaminant. However, 
the intensity of this group seems too large for this interpretation and it 
could also be due to the group «, from radiothorium corresponding to an 
excited state in ThX at about 218 kev. 

In view of this evidence and also of the fact that there is no indication 
from Rosenblum’s data (Rosenblum, Valadares and Perey 1949) of an 
alpha particle group corresponding to an excited state in ThX at 167 kev, 
the two most likely decay schemes are those shown in fig. 2 (a) and (b). 
The Bohr—Mottelson collective excitation theory (1953) predicts levels of 
a rotational character in even—even nuclei far removed from closed shells 
with energies #, spins J and even parity given by 
E= e J(J+1) where J=0, 2,4... and J is an effective moment of inertia. 


These states decay by cascade E2 transitions with energy ratios of 
3:7:11.... The first excited state of ThX is known to have spin 2 
and even parity (Rosenblum et al. 1952) so that if the simple rule above 
is followed exactly we should expect a level of spin 4 and even parity at 
282 kev. In fact, deviations from this simple rule are expected as a 
closed shell is approached, in such a way as to depress the higher levels 
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proportionately more than the lower ones. In view of these considerations 
and also of the fact that rotational levels are found to exist in other 
nuclei it is tempting to assume that scheme (a) is the correct one and 
that the level at 252 kev is the second rotational state. Further work 
is however required to establish this point. 

The 217 kev level in both schemes can be either 1— or 2+ ; spin values 
23— are excluded because the ground state gamma ray transition is 
observed to have the same order of intensity as that to the first excited 
state, and other values are excluded by the rules for angular momentum 
and parity. 


Fig. 2 


(4+?) 252 


I-, (2+) 


+ Ra 224 Ra22e4 


(Q) (d) 


Possible schemes for the decay of ?°8Th. The intensity of the a-particle group 
leading to the 217 kev state is estimated from the curves of Rosenblum, 
Valadares and Perey(1949). 


By extrapolation from the tables of Rose et al. (1951) it is possible to 
estimate the amount of K radiation from the various transitions. This 
gives a relative intensity of (2-1--0-4)10~% for the combined 132 and 
214 kev transitions if both are El, and (7-2--0-8)10-? if both are E2. 
This is to be compared with the upper limit of 5 10~? obtained experi- 
mentally for the total K radiation. This evidence suggests that the 
transition is E1 and therefore that the state at 217 kev is 1— ; but it 
cannot be considered conclusive since the K conversion coefficients are 

“not known accurately. The intensity ratio between the ‘214’ and 
132 kev gamma rays calculated from Weisskopf’s formula is 4:5 if both 
radiations are El, and 11 if both are E2 whilst the observed ratio is 
1:8+0:5. This again agrees rather better when both transitions are E1. 
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After this study had been completed we received a report of some work 
done by Asaro (1953) on this problem; he reached similar conclusions 
to those reported above. 

We should like to express our gratitude to K. L. Wilkinson and 
J. Waldron of the Chemistry Division at Harwell who carried out the 
separations. We also thank the Director, A.E.R.E. Harwell, for 
permission to publish. 

REFERENCES 


Asaro, F., 1953, Thesis. U.C.R.L. Report No. 2180. 

IBELING J. Ke lOb le ihesis) Mole 

Bock, R., 1950, Zeitschrift fiir Anorganische und Allgemeine Chemie, 263, 146. 

Bour, A., and Morrrtson, B. R., 1953, Kongelige Danske Viderslabernes 
Selskab. Mat. Fys. Medd., 27, No. 16. 

BoussizERRES, G., FatK-VARIART, P., Riou, M., Termuac, J., and Victor, C., 
1953, Comptes Rendus (Paris), 256,, 1874. 

Browne, C. I., 1952, U.C.R.L. Report No. 1764. 

Newron, J. O., and Rossz, B., 1953, to be published. 

Rose, M. E., Gorrtze., G. H., and SprnraD, B. I., 1951, Phys. Rev., 88, 79. 

RosENBLUM, S., VALADARES, M., and Perey, M., 1949, Comptes Rendus (Paris), 
228, 385. 

RosENBLUM, S., VALADARES, M., Perry, M., and Viau, J., 1949, Comptes 
Rendus (Paris), 229, 1009. 

RosENBLUM, S., VALADARES, M., and GuiLLot, M., 1952, Comptes Rendus 
(Paris), 235, 238. 

SuruGue, J., and San-Tstanea Tsien, 1941, Comptes Rendus (Paris), 213, 172. 

West, D., Dawson, J. K., and MANDLEBURG, C. J., 1952, Phil. Mag., 48, 875. 


eS er 


VUI. The Brownian Movement of Linear and Non-Linear Systems 
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ABSTRACT 


It has not so far proved possible to extend the conventional Brownian 
movement analysis to a system with an intrinsically non-linear relaxation 
mechanism (e.g. viscosity or electrical conductivity). In this paper a new 
statistical hypothesis is proposed which for the ideally linear circuit yields 
results in agreement with those already accepted and enables a calculation 
to be made of the Brownian movement (and corresponding frequency- 
spectrum) of a non-linear system. 

The results may be of interest in a number of practical problems, 
particularly in relation to proposals for the use of electrical ‘ noise ’ as an 
_ absolute thermometer at high temperatures. 


§ 1. InTRODUCTION 


SrtncE the first classic treatment of the Brownian motion by Einstein 
(1905, 1906), the problem has been treated by many writers on the basis of 
various assumptions. The analysis of Langevin (1908) is in its essentials 
frequently followed today (cf. e.g. Uhlenbeck and Ornstein 1930) ; therein 
the equation of motion of classical mechanics for a free particle of mass m, 
in a fluid of viscosity coefficient, f, say 


du 
m 7 =—fu, Eee ae «at C1) 

du 
is replaced by m a ae fai (bie Gs) (2) 
where w is the velocity of the particle (=dz/dt). To quote Uhlenbeck and 
Ornstein: “‘...Characteristically of this equation, the influence of the 


surrounding medium is split into two parts: 
(1) a systematic part —fu, which causes the friction 


b) 


(2) a fluctuating part F(t)...”. 

More or less detailed (but very reasonable) statistical assumptions about 

the character of this semi-molecular force are then made to enable (2) to 

be ‘ solved ’, yielding statistical parameters susceptible to experimental 
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observations, such as @, 22, etc. The average may either be one over an 
ensemble of similar individual systems, or be a time-average for a single 
(typical) system. 

This type of analysis has since been applied to many kinds of system— 
mechanical, electrical, and acoustical, for example—but it has generally 
been assumed that it is adequate to specify the frictional * force ’ as ideally 
linear in the velocity—i.e. that f in eqn. (2) is a constant. It is true that 
Kramers (1940) has pointed out that “... other laws of friction should 
also be possible... ”’, saying “‘... I do not know if the more complicated 
cases might have some physical application’. However, in his actual 
analysis (using the model of a particle under the influence of Brownian 
movement to discuss the rate of chemical reactions) he says ““... we will 
restrict ourselves to the simple Einstein case ”’ (i.e. of a linear friction). 

It appears to us that in electrical circuits, for example, which may 
involve essentially non-linear elements such as crystal rectifiers, an analysis 
embracing a non-linear resistance (i.e. friction) is now very desirable on 
many grounds; the more so in view of proposals to use electrical © noise ’ 
of thermal origin as an absolute thermometer for high temperatures. Also 
in problems involving radiative equilibrium of a body, it is evident that the 
radiative component of thermal conductance is intrinsically non-linear as a 
consequence of the quartic dependence on temperature (Stefan’s Law). 

It appears, however, that an extension of the analysis starting from 
eqn. (2) with —fu replaced by a more general function of velocity is rather 
intractable. We therefore propose a new statistical hypothesis of a 
non-molecular character ; if applied to an ideally linear system this yields 
immediately the usually accepted results, while in the non-linear problem 
it enables us to derive closed expressions for the Brownian movement 
(and the corresponding frequency-spectrum). The degree of validity 
of the hypothesis must naturally depend on future experiments, but it 
appears already clear how we may expect a non-linear viscous force to 
influence the character of the Brownian movement. 


§ 2. ANALYSIS 


Let us deal, for the sake of definiteness, with a simple electrical circuit 
(see figure) of capacity C, and resistance R, whose thermodynamic state 


“(== R 


is defined by the parameter q (the charge on the condenser at any time ¢). 
Now let us consider the fluctuating quantity {q(é) —q(0)? = Aq?. where 
the bar denotes an average over similar time-intervals, ¢, for a single 


Movement of Linear and Non-Linear Systems 65 


circuit or over an ensemble of similar circuits during the same time- 
interval, ¢.* Then 


Agq’=4(t)? +9(0)?—2q(t)q(0). ee eevee, (3) 
In statistical equilibrium, the first two terms on the right-hand side are 
equal and independent of t. Such terms can readily be evaluated quite 
generally if we accept, following Einstein (1906), that for a physical system 
placed in a ‘ heat-bath’ at temperature 7’, the probability that certain 
observable parameters—such as the displacement of a particle, or the 
electric charge in a circuit—lie between, say, « and «+d is given by 


pla) da=Ae~ PONT Gy, cine hoy tet) 


where (x) is the potential energy necessary to displace the coordinate « 
from its thermodynamic equilibrium value.t+ 

In order to deal with the time-dependent term (involving q(t)q(0) which 
remains in (3), we now propose the following general hypothesis : 

The average charge in an assembly of similar circuits subject to thermal 
fluctuations, all having a given initial charge (say q(0)), is governed by the 
same differential equation of ‘ relaxation ’ as holds for a single circuit in 

the absence of spontaneous fluctuations. 

Hence for an assembly of conventional linear circuits we have 


at Ro = © 
therefore om 
q(t)=4(0) exp (—t/RC), 
and 
q(t)q(0)=4q(0)? exp (—t/RC) 
Thus 


Aq 2=2¢?[1—exp (—t/RC)]. ear a s22 (6) 


But now from (4), since the thermodynamic equilibrium value of q is 
zero, we have ®(q)=q?/2C, and therefore 


= he exp (—q2/2CkT)q? ay] { exp (—@2/2CkT) dq=kTC. (7) 


Hence 


Ag 2=2kT Ci sexp (=.t) RO) Io 2. ey (8) 
For t< RC, this gives G 
2hT' 

Ad ae aha t, et Hie ioe LP ge PSEA sc (8 a) 


* The equality of such averages is essentially the assumption of the 
‘ ergodic hypothesis ’. ; 

This equation, although apparently closely related to the basic equation of 
statistical mechanics, presents certain difficulties in its derivation and 


application which will not be discussed here. 
F2 
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which is simply the well-known ‘ Brownian movement > equation of 
Einstein when applied to electricity. 
Further now consider 


2kT ; 
— Aq?= Aq? = cope exp (—t/RC). : en Sa ee (9) 


A transform theorem, derived earlier by the writer (MacDonald 1949) 


enables us then immediately to derive the power spectrum W(f), (or 2,) 
of the current fluctuation ; in fact 


Wif=4nf | Agesin 2nfdt, . . . . . . (0a) 
0 


which, on inserting (9), gives 


4kT ( (2nfCR)? 
W= =e oom} Gln Mabie Ki 


The spectrum of the charge-fluctuations is then given by 
qP=(1/2nf)P Wf) 
=4RkTC2/{1-+(2nfER)}y > «ee Oe) 


and we may check immediately that the tofal fluctuation of charge 
rc e | a if) —hTC. 
0 


Corresponding to (106), the voltage-fluctuations are given by 


ie 4RkT 


Vi=g7/C= T+QnfORP °° (11 a) 


which for f<1/CR gives simply 
Vi=4RkT os , Se 


usually known as Nyquist’s theorem. 


§3. THE Non-Linear PROBLEM 


Let us now generalize our circuit so that the electrical conductance, G, is 
no longer a constant but a function of the voltage across it and, therefore 
of the charge q(t) on the condenser. Let us assume Giq)=ayas 
Returning to eqn. (3) we may now write 


Ag P= 2gt2g(g(0), a. 
where the symbol ‘ signifies an average first over a sub-assembly 


of similar circuits which all have initially charge q(0)=q;._ The upper bar 
then denotes the final average over all sub-assemblies. 
Now in general, following our hypothesis, 


dq , gt). Gq) 
Apis ie CS ches od, ia 
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(in place of eqn. (5) for a linear circuit); whence, solving with G(q)=« + yy" 
and q(0)=q; we have 

_ exp (=at/C) 

1) = “OBE 


where B=(y/«)[1—exp (—2at/C)]. Hénce 


9 (upo_( EP (24/0) 
Aq?= 2 (ero (Gree )) 


; Pop aexp (97 KT Cg? d 
=2 {k7C—exp (—aii0) [" Se eee isaac) 


The upper integral may be evaluated in terms of the Bessel functions K, 
and we have 


Aj? =2kPCY 1—exp ( eee E — exp | ay K,—K,)—(4,+K, |}. 
(14) 
where a?={2kTC(y/«)[1—exp (—2at/C)]}-! and the Bessel function 


- argument is a?/2 throughout. 


We may now consider the behaviour of 4q¢/? for small and large values 
of a” either by use of appropriate expansions of the Bessel functions or by 
direct approximation to the integral in eqn. (14a) above. If we recall 
that G(q)=a+yq? and therefore G(q)=a+yq?=«[1+(y/a)kT'C], we see 
that the non-linearity will be small or large according as (y/«)kT'C is small 
or large eeabore’ with unity. 

If now 2(y/a)k7C <1, then a? is large for all values of t. Hence 


Aq? (a= BIC H1— < — exp (—at/C) ne exp (—2?)x (1 =) ae} 
AA MGs {1-exp (—at/C) | 1-5 5 ATC E [l—exp (—2a1icy1 | (15 
Evidently, as we expect, in the limit y/x=0 (i.e. zero non-linearity), 


Aq?=2kTC[1—exp (—at/C)], in agreement with (8) above, setting 
R=1/«. More generally, 


: 9 
Aq? =20kT { (1- z kTC ") exp (—at/C)+ 5 kTC E exp (—3ati0)| 


therefore using the transform eqn. (10 a) we find the ‘ noise ’ spectrum of 
the current 


(w/a)? 3 a { d(wC/3a)? (wO/a)? 
WES et {ort Se feemanne 33 | 
(16) 


where w=2nf. These expressions, with their non-linear temperature 
dependence, then appear of particular interest in connection with the 
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proposals, mentioned earlier, of using electrical ‘ noise’ as an absolute 
thermometer at high temperatures. 

The spectrum of the charge fluctuations is given, as before, by 
W(f)/(27f)2 and we may immediately verify that again the ‘ total ’ 
fluctuation 


re (= I, G if) =kTC. 


The case of strong non-linearity [(y/«)kTC>1] requires somewhat 
more detailed examination since it is evident that a? cannot -> 0 for all 
values of ¢t. We find that for t<1/ykT 

‘Aq 2=2kTC[1—exp (—at/C)] 
and for t>1/ykT’ 

eS 2 exp (—at/C) 

2__ 9], y es |] FT 

lA {1 (nkTO(y]a)[1—exp (—2ad/C)] = 
These two expressions are, in fact, very nearly equal for t=1/7ykT, and 
we may therefore expect that an excellent approximation will result by 
extending both to this common limit. However, the subsequent deter- 
mination of the noise spectrum will not be attempted. 

Finally, we may analyse readily the behaviour of an idealized ‘ rectifier ’ 
—a problem discussed qualitatively in the past. If we let 


(17) 


G(q)=«. ; q>0,* 
ie q<0. 
Then 
ss x rp . __ BHO 
AqP=2kTC} 1— ph esp ta 
therefore Agf=2kT 5 exp (—at/C)+ z exp (—ay0)} bey en 


Whence the spectrum of the noise current is 


Nair oe Oia aR agtOC Bs 
WN=AkP 15 aap THC aR: Si 
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SUMMARY 


A hypothesis recently advanced by MacDonald is shown to lead to a 
contradiction. It is indicated how a formalism developed by Kramers 
may be used in this connection. 


$1 

In the preceding paper MacDonald has drawn attention to the fact 
that relatively little is known about the theory of the Brownian motion 
in non-linear systems. As an example MacDonald considered a simple 
#-C circuit in which the resistance R is a function of the charge g on 
the condenser C. In order to be able to work out the noise spectrum 
of such a circuit MacDonald advanced the following hypothesis: ¢, the 
average value of q of a great number of identical systems, all starting 
at ¢=0 with the same value of g=q satisfies the equation 

(Gage GiGaOnw ee i (1) 
In other words, and slightly less precisely, it was assumed that the time 
dependence of the average charge of the ensemble is the same as that of 
the charge of a single system not subject to Brownian motion. 

It is obvious that for the case R is a constant hypothesis (1) is in 
agreement with existing theoretical ideas about the Brownian motion 
in linear circuits. At first sight the extension of eqn. (1) to the non-linear 
case seems to be a hypothesis of which the validity is no more questionable 
than that of any of the assumptions usually made in the theory of 
Brownian motion, i.e. the test of its validity depends on a detailed 
molecular theory of the relevant random process. However, this is not 
the case. It can be shown that MacDonald’s hypothesis is in fact not 
self-consistent. This will be done in §2. In §3 we shall indicate how a 
paper by Kramers (1940) may help to illuminate the problem under 
discussion. We shall not discuss questions concerning the validity of 
the usual assumptions in the theory or the problem of the minimum 
number that suffices to set up the theory, though both points might be 
relevant to a satisfactory discussion of the problem of the noise spectrum. 


§2 
In order to test the self-consistency of MacDonald’s hypothesis (1) 
we consider a special case in which for g>0 RF is equal to a constant R,, 
and for q<0 R is a constant R,; k,>R, say. If at t=0 all systems of 
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the ensemble start at g—=qo, eqn. (1) predicts that g=qo exp (—t/R,C) 
since g always remains in the positive region of g. After a finite time 7’, 
however, we must expect to find in the ensemble a certain distribution 
of g values around g. More in particular a certain fraction f of the 
systems will have drifted into the region <0, this fraction being the 
larger, the larger 7’. Now let us split up the ensemble at =T' into 
sub-ensembles of systems, each to be found in different intervals 4g, 
together covering the whole range of q-values. According to the 
hypothesis (1), all sub-ensembles in the positive region will again show 
the exponential decay exp (—t/R,C). The sub-ensembles in the region 
q<0, however, will show a decay with the characteristic time R,C. This 
means that for (>7' we now have 


gq=—a exp (—t/R,C)—b exp (—#/R,C), ~ 2. ee 


where a and 6 are certain constants >0. This obviously contradicts 
the previous deduction from the theory. 

As to the origin of this failure, we wish to emphasize that (1) has been 
introduced as a purely ad hoc hypothesis, based on similarity rather 
than on existing theory. Nevertheless it may be of interest to discuss 
one point that should not be used as an argument against (1). It may 
seem surprising that the initial condition on the members of the ensemble 
specifies only g=q) and not the value of the current 7=dgq/dt. The 
answer to this is that the circuit should really be considered as one with 
an extremely small self-inductance L in series with R and C, with L—0. 
Now, whatever is the initial current 7) of the members of the ensemble, 
a distribution of currents around an average —4q/RC will have established 
itself within an extremely short time of the order L/R. It is this 
distribution which alternatively may be taken as the initial condition 
for the currents. The width of this current distribution is, of course, 


very wide, (47)?~kT/L but the width of the q-distribution after that 


short time is still very small (4¢)?~2kTt/R~2LkT/R? as required in the 
hypothesis. 


$3 

It may be of interest to point out that a non-linear problem has been 
discussed before in a paper by Kramers (1940). In fact his treatment 
concerns an even more general case : a mass point subject to a frictional 
force —B dq/dt, and a systematic force K, where both 8 and K are functions 
of the coordinate g. This corresponds to an electric circuit with 
constant L, but Rk and C depending on g. Kramers soon specializes to 
the case of constant 8 and variable K, but it is not difficult to see what 
happens if we keep B=B(q) and let m0, thus returning to MacDonald’s 
problem with also the capacity depending on the charge. 

For details we must refer to Kramers’ paper, but we mention that 
the fundamental assumptions in this paper are, as usual, certain state- 
ments about the integral over a small time 7 of the fluctuating force 
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acting on the particle. This time is long compared with the time scale 
of a molecular collision, but so small that the velocity of the particle 
has not substantially changed (i.e. smaller than the small time mentioned 
at the end of §2). This, incidentally, means that m—->0 should be 
considered as m is macroscopically small, not as a limiting process into 
the molecular region. The average of the integral, of its square, etc., 
over an ensemble with identical initial g and dgq/dt, are supposed to be 
developable into a power series of +. The terms proportional to 7 involve 
the value of f at the initial value of gq and the q-dependence of 8 
affects only the higher order terms. Only the terms proportional to 7 are 
of interest for setting up a differential equation for the distribution 
function p of any arbitrary ensemble in phase space. Kramers’ equation 
for an arbitrary ensemble runs as follows : 


a 


The terms on the right-hand side describe ee be effect of the 
velocity p/m, of the systematic force K, of the average of the fluctuating 
force or frictional force and of the Brownian motion. 

Following Kramers with a minor variation, we introduce as variables 
@ and p instead of ¢ and p, where 


@= |" BU’) aa’ +p. . Ten Aig Oe OP, 5 


Introducing partial derivatives with respect to the new coordinates we 
find from (3) 


lap o/ p x(x 5) — 3a ( s) (= x) 5 
3 Hap (Pn) ~ ap (KB) aq XB) +4? (Se ta@) & ©) 


Again following Kramers we integrate over p at constant Q’. If m is 
very small, constant Q’ means approximately the same as constant gq, 


since the spread in p is at most of the order VmkT. Thus we have for 
any, not too quickly varying functions f of q, if m0 


[olw, Qf (@) dp>of atq=qQ.p=0) . . . . (6) 


where o is the ensemble density in configuration space. The integration 
cancels the terms with 0/dp since o+0 for large p and we have 


1 do 0 (Ko 0?a 
5a SS kT —,. a 
parole) Hae o 
Re-writing with the aid of ¢, we finally have 
da 


= =— al a(Ke -«r) |. ih ek peta oor) 


This is a well-known diffusion equation in q-space. It is easily seen that 
the Boltzmann distribution is a stationary solution if AK is a conservative 
force. 
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Equation (8) gives, in principle, the answer to all questions regarding 
the time-dependence of any average quantity of interest. If the initial 
d=; we have to find the solution of (8) for c=6(q—q) at t=0. It is 
not unlikely, however, that it will be very difficult to find an explicit 
solution in all cases but in the well known simple ones. It is also 
evident, however, that a simple suggestion such as given by MacDonald 
does not provide the correct answer to the problem. 


The author wishes to acknowledge discussions with Dr. D. K. C. 
MacDonald and to express his thanks for having had the opportunity to 
see his paper before publication. 
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ABSTRACT 


The electrical and thermal resistivities of the same specimen of 
magnesium have been measured together between 2 and 35°x. The 
electrical resistance had a minimum at 6°K and an effect of a similar order 
of magnitude occurred for the thermal resistance W. This is best demon- 
strated by plotting a graph of WT against 7? which shows a marked 
departure from linearity in the region of the electrical resistance anomaly. 


§1. INTRODUCTION 


In some samples of certain metals a minimum is known to exist in the 
electrical resistance at low temperatures. The experiments to be des- 
cribed were designed to see whether there was a corresponding effect in the 
thermal resistance. An anomaly of a similar order of magnitude has 
indeed now been found in a batch of magnesium, which was known to 
exhibit a minimum in the electrical resistance (MacDonald and 
Mendelssohn 1950). 

The only previous work on these lines is an observation by Berman and 
MacDonald (1952) on a copper specimen which showed an electrical 
resistance minimum. In this case, however, the authors stated that any 
deviation of the thermal resistance from its normal curve must have been 
considerably less than the anomalous increase in electrical resistance and 
that a careful analysis of the results was necessary to indicate its existence. 


§ 2. EXPERIMENTAL TECHNIQUE 


Earlier work on the correspondence between electrical and thermal 
conductivities has on some occasions been rendered difficult of interpre- 
tation because the thermal and electrical measurements were made at 
different times and in some cases separate specimens have been used for 
each measurement. Even when the same specimen had been used for 
both the thermal and electrical experiments it has usually been handled 
between the two determinations thereby rendering it liable to a small 
amount of distortion which often will affect the conductivities considerably. 
For this reason it was decided to measure both the thermal and electrical 
conductivities together on the same specimen and during the same experi- 
mental run. The thermal conductivity was measured in a way similar to 
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that described in our previous work (Mendelssohn and Rosenberg 1952). 
The accuracy of the thermal conductivity measurements over most of the 
experimental range was about 1%. 

The electrical resistance of the specimen was measured by using a 
galvanometer amplifier (MacDonald 1947) to determine the potential 
difference between the gas thermometer contacts when a current of about 
0-5 amp was passed through the specimen. The constantan potential 
leads were actually connected to the gas thermometers. In this way the 
form factors for both the electrical and thermal measurements were the 
same and no errors were introduced in comparing the two sets of values and, 
in particular, when calculating the Lorenz number. It might be thought 
at first that the German silver capillaries attached to the gas thermometers 
which are anchored together to the liquefier framework higher up in the 
apparatus might act as a short circuit to the potential measurement. In 
actual fact, however, the resistance of this ‘ short circuit ’ between the gas 
thermometers is about 2 ohms whereas that of the specimen is about 
10-® ohms. Hence no appreciable error is introduced. Since the cold 
end of the specimen is attached to the metal apparatus and to earth, the 
latter was used as one current lead whilst a 38 gauge copper wire attached 
to the warm end of the specimen was used as the other connection. 

The specimen and the gas thermometers were mounted in a high vacuum 
enclosure and were surrounded by a copper radiation shield which was 
thermally attached to the helium container. Such an arrangement 
whilst being the most favourable for the measurement of the thermal 
resistance is not so satisfactory for the electrical measurements because the 
measuring current inevitably introduces a small heating both in the 
specimen and in the current lead which can be appreciable as the specimen 
is thermally isolated. Thus thermal e.m.f.’s might be produced which 
would lower the accuracy of the measurements. This will be appreciated 
when it is realized that the potentials to be measured at liquid helium 
temperatures were only of the order of 3x 10-° volts. For this reason the 
measuring current had to be made as small as possible and hence the 
relative accuracy of the electrical measurements was about 3°% at liquid 
helium temperatures. Since the electrical measurements were undertaken 
mainly in order to check that this specimen did in fact exhibit an electrical 
resistance minimum such an accuracy was considered satisfactory. 

The specimen was made from Johnson Matthey spectrographically 
standardized rod, JM 1703 which had a purity greater than 99-95%. 
The main impurities were Mn, 0:03; Fe, 0:0075; Al, 0-004. It was 
turned down to a diameter of 1-5 mm and was annealed in vacuo at 500°C 
for six hours. This rod was made of the same batch of metal as were 
two of the specimens which MacDonald and Mendelssohn (1950) used to 
determine the electrical resistance and since the present electrical measure- 
ments do correspond with theirs, a more detailed comparison of the 
thermal and electrical measurements can be made by referring to these 
previous electrical results. 
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§ 3. RESULTS 


The thermal conductivity K, and the electrical resistance R, as measured 
in these experiments between 2 and 35°K are shown in fig. 1. In spite of 
the scatter a minimum in R is clearly seen at about 6°K. At first sight the 


Fig. 1 


Thermal conductivity, K watt/cm—deg. 
Electrical resistivity, R ohm cm x 108 


Temperature °K 


The thermal conductivity A, and the electrical resistivity R, of the magnesium 
specimen showing the departure of K from linearity in the region of the 
electrical resistance minimum. 


thermal conductivity curve has the appearance that one would expect 
from a normal fairly pure metal. There is a linear section from 2 to 10°K 
passing through the origin, showing electron scattering by the lattice 
impurities, followed by a maximum and then a decrease in K where the 
electrons are scattered by the lattice vibrations. A theoretical study of 
these mechanisms has been given by Wilson (1936), Makinson (1938) and 
Sondheimer (1950) and a general review of the theory by Olsen and 
Rosenberg (1953). 

A more careful study of the curve of K against 7’, however, shows that 
the lower temperature section is not quite linear, but that between about 
4 and 9-5°x it is slightly higher than it ought to be. This departure from 
linearity is much better demonstrated in the following way. 
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The thermal resistance, W, can be written as the sum of two terms at 
low temperatures, 

W=AT?+B/T oi 2 
where, according to theory, A and B are constants of the metal and its 
impurity content respectively. 

Thus 

WT=AT®LB:. i oe 
and hence a graph of WT against 7? should be a straight line with a slope 
of A and an intercept on the W7' axis of B. This has been found to be 
generally true for most metals which have been investigated and is the 
standard method for determining the constants A and B. 

When such a plot is made for our magnesium specimen, however, as is 
shown in fig. 2, it is evident that there is a marked departure from a 


Fig. 2 
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The graph of W7' against 7 showing the minimum in the value of W7. 
Inset, the lower half of the curve on a larger scale. 


straight line below about 11°K. With falling temperature the value of 
WT drops below the straight line until it reaches a minimum at about 6°K 
and then rises again. At the lowest temperatures W7' again agrees with 
eqn. (2). Above 11°k A and B in eqns. (1) and (2) have the values 
8-5 10~° and 1-05 respectively. 

A graph of the difference between the observed WT and AT8+B is 
shown in fig. 3. The scatter in the points is understandable since we are 


Conductivity of Magnesium at Low Temperatures ae 


dealing with the difference of two quantities which are nearly the same 
and this difference is very near to the experimental accuracy. From the 
smoothed curve of fig. 3 a graph of 4W, the decrease in resistance from the 
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The graph of 4(W7') against 7°, where 4(W7’) is the difference between the 
theoretical equation A7%+B and the experimental results. 
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theoretical value of eqn. (1) against temperature is shown in fig. 4. 
Whilst this curve is approximately parabolic in shape a detailed analysis 


shows that it is not strictly so. 
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It can be seen that the minimum in the value of 4W occurs at 6°K 
which corresponds exactly to the position of the minimum in the electrical 
resistance both for our, specimen and also for those of MacDonald and 
Mendelssohn. Furthermore the value of 4W at this temperature is about 
4°/, of the value of the absolute thermal resistance. For the electrical 
resistance, R, the minimum is about 3% below the residual resistance as 
deduced from the measured points above 10°K by plotting & against 7”. 
Thus whatever mechanism is producing the minimum, it seems to affect 
the thermal and the electrical resistances equally. 


EL (watt ohm/deg.?) x 108 


O 10 20 30 


Temperature °K 


The temperature variati ian ; 
P ariation of the Lorenz number L, for the magnesium 
specimen. 


It has sometimes been suggested that the rise of R below the minimum 
may lead to infinite resistivity at the absolute zero (cf. Gorter 1938) 
An indication of the low temperature behaviour of R is provided by the 
WT against 7° curve of fig. 2. It will be noted that at the lowest temper- 
atures measured the values of W7' again fall on to the theoretical Bis ht 
line. Since the present experiments show that W varies in a similar ane 
with temperature as does R, they suggest that R, too, will as absolute oie 
is approached, return to its theoretical temperature-independent valet 
However, it should be noted that recent experimental work (Mendoza and 
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Thomas 1951, Croft, Faulkner, Hatton and Seymour 1953) still shows R 
increasing at he lowest ee atures reached. 

As a consequence of the fact that comparable effects were found in both 
the electrical and thermal resistances no anomaly occurs in the curve of 
the Lorenz number, L, against temperature, where 


L=R|WT. 


In experimental range covered, the value of LZ, shown in fig. 5, rises 
with decreasing temperature, dees off to a constant value at about 
5°K, below which temperature the Wiedemann—Franz law is valid. 
Whilst this constant value of L is slightly higher than the theoretical value 
of 2:45 x 10-8 watt ohm/deg.? this increase is not thought to be significant 
to this problem. 

The minimum in the electrical resistance is usually very shallow and is 
difficult to detect. The deviation from linearity of the W7' curves gives 
another method for showing the existence of such a minimum. It was 
therefore thought advisable to examine the results of thermal conductivity 
measurements which we have made previously on pure metals (Mendelssohn 
and Rosenberg 1952) to see whether any of them exhibited a slight dip in 
the WT against T° curve at low temperatures as does our magnesium 
specimen. In four cases such a deviation could be detected : in copper 
and silver (samples of which are known to exhibit a minimum) and in 
nickel and palladium. Further careful measurements however, would be 
necessary to confirm the existence of a minimum in the last two cases. 

It is apparent that more data on the nature of these minima are 
necessary before the reason for their existence is made clear. At the 
moment, whilst it seems likely that their presence is connected with certain 
types of impurity atoms (e.g. cf. Mendelssohn 1951, MacDonald 1952) a 
full explanation of the phenomenon is not yet possible. The present 
result, that the effect is not limited to the electrical resistance, but that 
the thermal resistance is altered in a similar way, may widen the field of 
investigation. 
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ABSTRACT 


The theory described in this paper was developed with a view to its 
application to the molecular theory of helium II. This application will be 
discussed in a subsequent paper. 

The system treated here consists of a very large number of identical 
interacting particles. Sound waves can be described using the Fourier 
components of the density operator and their time-derivatives. Approxi- 
mate commutation relations for these operators are obtained, and 
equations of motion are deduced from the assumption that the system 
can propagate sound. The commutation relations are analogous to those 
for a set of harmonic oscillators or for a radiation field; using this 
analogy, creation operators, which can be said to create phonons, are 
defined. It is shown how stationary states may be constructed by 
adding phonons to the ground state, and that the energy or momentum 
of such a stationary state equals the energy or momentum of the ground 
state plus that of the ‘ phonon gas ’ added to it. 

The theory is illustrated by its application to three particular cases— 
the imperfect Bose gas, linear imperfect Fermi gas, and linear crystal. 
In each case the present general method leads to approximately the 
same creation operators as the special treatments due to other authors. 


$1. INTRODUCTION 


In his well-known treatment of the specific heats of solids, Debye (1912) 
introduced the approximation of treating a many-particle system as an 
elastic continuum. Thermal motion is represented by density and shear 
strain fluctuations in the continuum. This motion of the continuum can 
be analysed into transverse and longitudinal sound vibrations, and the 
system therefore has the same energy levels as an assembly of harmonic 
oscillators, one for each sound mode. At low temperatures, the part of 
the specific heat due to longitudinal modes is 

(2777/15) («T'/he)>« ge ee 2, 
per unit volume, where c is the velocity of longitudinal sound waves, « is 
Boltzmann’s constant, and 7’ is the temperature. 

Born and von Karman (1913) justified Debye’s approximation, by 
showing that the motion of a crystal can be analysed into normal modes of 
vibration, similar to sound waves, and that Debye’s approximation gives 
Se ee a eee 
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the correct frequencies for long wavelengths. It follows that (1) does 
apply to a real crystal at sufficiently low temperatures. 

In his work on liquid helium IT, Landau (1941) followed Debye’s example 
and treated this many-particle system as a continuum—an ideal fluid 
rather than a solid. Just as in Debye’s work, the irrotational part of the 
motion is analysed into longitudinal sound waves and gives the 
contribution (1) to the specific heat. The incompressible part of the 
motion, however, cannot be treated by Debye’s method because an ideal 
fluid cannot support transverse vibrations. Landau claimed that the 
degrees of freedom associated with incompressible motion do not contribute 
to the specific heat at very low temperatures, so that (1) gives the entire 
specific heat at such temperatures. This conclusion was verified experi- 
mentally by the specific heat measurements of Kramers, Wasscher, and 
Gorter (1952). 

A partial justification of Landau’s assumptions was given by 
Bogolyubov (1947). Bogolyubov showed that the Hamiltonian for a 
degenerate imperfect Bose gas can be approximately expressed as a sum of 
commuting harmonic oscillator Hamiltonians. As one would expect from 
the application of Landau’s theory at very low temperatures, the lowest 
frequencies are distributed as if they were the frequencies of longitudinal 
sound modes, and there is nothing corresponding to transverse sound 
modes. 

Results confirming the value of Debye’s point of view have been obtained 
for two other types of many-particle system. Tomonaga (1950) showed 
that the Hamiltonian for a one-dimensional imperfect Fermi gas could be 
expressed as a sum of commuting harmonic oscillator Hamiltonians, with 
the low frequencies distributed in the way suggested by Debye’s theory. 
And Pines and Bohm (1952 a, b) showed that harmonic oscillator terms are 
present in the Hamiltonian for an electron gas, though other degrees of 
freedom are also present. 

The success of Debye’s method for such a variety of systems suggests 
the possibility that similar results would follow for any system in which 
sound can be propagated. This conjecture is confirmed by the general 
theory described in the present paper, which shows that the presence of 
harmonic oscillator terms in the Hamiltonian can be inferred from the 
periodic motions which constitute sound. The theory has been illustrated 
by showing that the essential results obtained for the above-mentioned 
special models by Born and von Karman, Tomonaga, Bogolyubov, and 
Pines and Bohm can also be obtained by applying this general theory. 


§ 2. Microscopic DESCRIPTION OF THE SOUND FIELD 


As a model of a physical system, consider a very large number, NV, of 


interacting particles, each of mass m, with position vectors qj, ... qy, 
momentum vectors p;,..- Py, and Hamiltonian 
N 
H= Spine 2 UG; aae = as get?) 
j=l i<j 


G2 
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Dirac’s notation (1947) will be used. To simulate the effect of the walls 
of the container, periodic boundary conditions will be imposed : that is, 
the coordinate representative of any permissible state or operator must be 
periodic in each of the q,’s. For example, in order for the Hamiltonian (2) 
to be periodic, U must be a periodic function of its argument. With these 
boundary conditions, each of the momentum operators p; has a discrete 
spectrum of eigenvalues. The system is effectively confined to a volume 
V, where V is the repeating volume of the periodicity. 

Landau (1941) defined the density operator p(x) and the current density 
operator J(x) : 


p(x 25(q,—%), 
J(x)= 23 {p; 8(q;—*) +8(4;—*) P;}; | 


where 5(x) is a three-dimensional Dirac delta-function. Their expectation 
values are respectively the density and the current density at the point x. 
They satisfy the equation of continuity 


(3) 


p(x)=—div ](x)7 29. 2) =) eee 
with Heisenberg’s definition of a time-derivative 
A=|A\ HI=(AH_ HAA. 


The Fourier components of these operators, 


A= | exp (7k . x)p(x) d’x=m 2’ exp (ik . q;), 
Vv j 
(6) 
Je= |_ exp (ik. x)J(x) x= T4{p, exp (ik «4,)-+ exp (ik. .)P,} 


therefore satisfy the equations 


Py=tk . la . e . . ° . . . e (7) 
If k40, the operators py, p_~: Pe: and p_,. describe the pair of longitu- 
dinal sound modes whose wave-vectors are +k (the wavelength is 
2n/k). Moreover, pp=Nm is the mass of the system, and J, is its 
momentum. 
From (6), (7), and the fundamental commutation relations between 
the position and momentum variables, the following commutation 
relations can be deduced : 


JoPu— PrJo=hk py, | 
JoPu— Prdo=Nk py. J 
These will be used in discussing the momentum of sound quanta. 


(8) 


§ 3. THe Sounp PropaGation Equation 
If sound vibrations are present in the system, py will vary approxi- 
mately sinusoidally with time, with frequency kc/27, where k is the length 
of the vector k. For a system at rest, this fact can be expressed by the 
oscillation equation 


Puc kee Oke ene te 
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This equation can be derived from the equations of hydrodynamics, or 
from the experimental fact that sound waves are observed. For the 
theoretical models treated in §§ 6 and 7 it will be proved directly. 

If dispersion is negligible, ¢ is independent of k and can be calculated 
from the adiabatic compressibility. For the present purpose, the 
compressibility at absolute zero will be used; consequently (9) only 
applies for temperatures near absolute zero. Equation (9) neglects 
attenuation, fluctuations, and non-linear effects such as the coupling 
of different modes ; since attenuation increases rapidly with k, the approxi- 
mation (9) is only permissible when 


k<hy. ee Pee! oe Pe (10) 


Here ky is some limiting wave-number. Since disturbances of shorter 
wavelength than the mean separation of particles will not behave like 
sound waves, ky, cannot exceed (N/V)¥/3, 

Since non-linear effects were neglected in (9), the amplitude of sound 
must be restricted to small values. That is, neither the density nor the 
current density may depart greatly from its value averaged over the 
whole volume V. Since the average velocity is 


w= J,/Nm, ee me oe, LL) 
this assumption can be written 
| p(x) —Nm|V | <NmV, 1 | 
[ J) —E{wp(x) + p(xyw}| <Nmelv. fo 7 OF 


The form of the second inequality has been chosen so as to make it 
independent of the average velocity w, which is irrelevant to the internal 
motions of the liquid. It follows from (12) that 


| pe | = || exp (ik . x){p(x)—Nm/V} d8x| <Nm,...(k40), 


| JE {Went pew} | <Nme. 


The assumption (13) need only be made for k<ky,. 
It is possible to generalize eqn. (9) to moving systems. The unitary 
operator 


(13) 


G,= expXimu.q,/h, Re go oe oe (14) 
j 


transforms the variables of the system as follows : 

Gui Fu=4s> 

G_yPj;Gu=P; +m, 

GaGa Jae ee eb) 

G_,HG,=H+u. Jo+3Nmv?. Soe te 
The boundary conditions only permit values of u such that mu is a 
possible value of the momentum of a single particle. The transformation 


corresponds to changing to a frame of reference whose velocity is —u 
relative to the old one, or alternatively to increasing the velocity of every 
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particle by u. Neither the Hamiltonian H nor the total momentum J, 
is invariant under this transformation, but the operator 
H=H-JA2Nm VA 
is so invariant. It can be regarded as the Hamiltonian referred to a 
frame of reference moving with the velocity w of the system (where w is 


given by (11)), and will be called the internal energy. 
A new kind of time-derivative can now be defined : 


° 


A= [A, 9] 

=A—d{w .[4, Joltl4, Jel: We 
e.g. p(x)=p(x)+4{w . grad p(x)+ grad p(x) . w}. 
Since, by (12), the local velocity is approximately w everywhere, this 
last equation shows that the new time-derivative is approximately 
the same as the derivative following the motion of the fluid, or convective 
derivative, used in hydrodynamics (Lamb 1924). Equation (9) may 
therefore be generalized to moving systems by replacing the ordinary 
time-derivatives by the convective derivatives defined in (18) : 


byt hect py Omen se) en 


Each term of this equation is invariant under transformation by G,, so 
that eqn. (19) is a statement of the sound propagation equation that is 
independent of the frame of reference. 
The convective derivative of p, can be found from (18), (8), (6), and 
(7): | 
pr=ik «2X 3{(p;—mw) exp (ik .q,)+ exp (ik . q;)(p;—mw)}. (20) 
J 
The sets of operators p, and p,, being invariant under transformation by 
any Gy, are suitable for a description of the sound field that is independent 
of the velocity of the system. Their commutation relations, which follow 
from (6), (20), (11), (8), and the fundamental quantum conditions, are 


PaPb— PpPa—), 
PaPb— PpPa=—th(a. b) {a+b — Papr/ Nm}, 
PaPb—PbPa=h(a. b)(a—b) . (Jasb—2{WPa+ptPa+pW}) 
+ih(a . b)(pap,—Prpa)/Nm. 
On account of (13), these simplify to 
PaPb— PoPa= 0, 
PaP-a—P_aPamihNma?, 
| paby— Pubs |<hNenab Gf vadtbs’o. | eee 
| PaPe—PoPa| <hNmab(a+b)c. | 
Thus, p_,/Nmk is an approximate canonical conjugate variable to p,. 
This is analogous to the situation for a harmonic oscillator, where the 


canonical conjugate of either the position or the momentum variable is 
proportional to the time-derivative of this variable. 
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Kronig and Thellung (1952) have given a field-theoretical treatment 
of the irrotational motion of a perfect fluid. Their treatment shows that 
the density field possesses an exact canonical conjugate—the velocity 
potential. The canonical conjugates of the p,’s are therefore the Fourier 
components of the velocity potential. These are indeed proportional to 
the approximate time-derivatives of the p,’s obtained from the approxi- 
mate Hamiltonian which contains only terms of the second order in the 
px S and their conjugates (eqn. (39) of Kronig and Thellung’s paper). 
Since the number of p,’s is infinite, the fact that canonical conjugates can 
be found for all of them simultaneously means that the ideal fluid has an 
infinite number of degrees of freedom. 

The system of particles considered in the present paper has only a 
finite number (3) of degrees of freedom. Consequently it is not 
possible to find exact canonical conjugates for all the p,’s simultaneously. 
Thus, the approximate nature of the commutation relations (21) can 
be thought of as resulting from the attempt to describe a discrete system 


of particles in a way more appropriate for the description of a continuous 
fluid. 


§ 4. THE CREATION AND DESTRUCTION OPERATORS 


Using the analogy between a sound mode and a harmonic oscillator, 
mentioned above, a quantum theory of sound waves can be built up, 
similar to Dirac’s treatment of the harmonic oscillator (1947 a). Dirac’s 
treatment is based on ‘ creation’ and ‘ destruction’ operators, which 
transfer the system from one stationary state to another. It will be 
shown here that the corresponding operators for a sound mode are 


n= (2N bike) 1p, +-ikepy.); 


Z we i i eee 22) 
= (2Nmhk3c)-Y?(p_,—tkep_,). ( ssf 


Very similar operators have already been given by Akhieser and Pomeran- 


chuk (1945). r ' 
Because of (19), the 7,’s and 7,8 have the equations of motion 


Nyeiken, ; eo —tken,, . . . . . . (23) 


and, because of (21), they satisfy the commutation relations 


Talo "Iba =O, l oy 
Tla"lb— NbNa=— Sap: | 
Thus the 7,’s and 7’s approximately satisfy the commutation relations 
characteristic of the creation and destruction operators for a set of 
harmonic oscillators or for a set of Bose-Einstein particles (Dirac 1947 b). 
These operators 7, and 7, may therefore be said respectively to create 
and destroy sound quanta or phonons in the mode k. 
The operators defined by 
Pekka ee oe oe ee) 
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satisfy, by (24), the relations 
ida = Malet Oak) | 

“Na = Na(’e—Sak)- 
Equation (26) means that if a state |v’) is simultaneously an approxi- 
mate eigenstate of all the v’s, corresponding to the set of eigenvalues 


ee Dal sss p ss... (this is possible because the v’s commute), then 

m|v’) and 7,|v’) are also eigenstates of all the v’s, corresponding to 
° , Jt 

the sets of eigenvalues ...v,’... vp’... (4% +1)... and... .-- Mm +> 


(,’—1)... respectively. The form of eqn. (25) shows that the expecta- 
tion of », cannot be negative, so that », has no negative eigenvalues. 
Consequently, the sequence |v’), |v’), m|v’)... cannot continue 
indefinitely ; it must end with some state that is annihilated by 7. 
By (25), such a state must be an eigenstate of »,, with the eigenvalue 
zero, and therefore »,’ can only take the values 0, 1, 2,3.... (This proof 
is due to Dirac (1947 a).) It is therefore natural to interpret 4, as the 
number of phonons in the mode k. The interpretation of , and 7 as 
operators creating and destroying one phonon in the mode k is confirmed 
by this discussion. 

Quanta of energy and momentum may be associated with each phonon 
in the following way: eqn. (23) is equivalent to 

On, = 1 (O--hhe). eg ee 

Thus, if | 5’) is an eigenstate of the internal energy operator 5, corres- 
ponding to the eigenvalue ’, then », |’) is also an eigenstate of 5, 
corresponding to the eigenvalue )’-+-Akce. In this way, the addition of 
a phonon in the mode k is accompanied by an increase of energy hkc, 
which may be called the energy of that phonon. Similarly, the equation 


Jom =7u(Jo+hk), 2 ee ae as 
which follows from (8) and (22), can be interpreted to mean that a 
phonon in the mode k has the momentum ik. This momentum of a 
phonon is related to the Doppler effect : by (18) and (28), 
Me= Me — 31K. (Wt 7W}s 

and therefore, by (23), 

Me= hi {(ke+k . W) n+, (ke+k . w)}. Meee he LY, 
Thus, the frequency of sound in the mode k, when the velocity of the 
system is w, is (ke+k.w)/27. This value agrees with the Doppler 
effect formula. 

The creation operators can be used to construct stationary states 
from the ground state. The ground state | 0) is defined as the eigenstate 
of H corresponding to the lowest eigenvalue, Ey. For simplicity this 
energy level is assumed non-degenerate. Since H commutes with the 
momentum Jy, | 0) is also an eigenstate of J,; the corresponding eigen- . 
value must be zero because the properties of the system are unaltered by 


Se ee : 
reflections. Consequently |0) corresponds to the eigenvalue Hy of the 
internal energy operator §, 
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One may also define ‘ moving ground states ’, as follows : 
where G, is defined by (14). Because G, commutes with §, while J, 
is transformed by it according to (15), | u) is an eigenstate of $ and of 
Jo, corresponding to the eigenvalues EZ, and Nmu respectively. 

Now, consider the state 7,|0). Equations (27) and (28) apparently 
indicate that this is an eigenstate of 5 and of J, corresponding to the 
eigenvalues H,—fike and —hk respectively, and hence of H corresponding 
to the eigenvalue Hy)—fhkce+2k?/2Nm. This quantity, however, cannot 
be an eigenvalue of H because it is less than Hy ; consequently, | 020 
must hold. Since G, commutes with the p,’s and p,’s, and therefore 
with the »,’s, it follows that 

TOKO OMe eR, eae ays ee get (30) 
This implies that | u) is an eigenstate of all the »,’s, corresponding to the 
eigenvalue zero ; that is, no phonons are present when the system is in 
a moving ground state. 
By operating on any moving ground state |u) with a succession of 
7 8, a State may be constructed which is an eigenstate of J, and an approxi- 
mate eigenstate of 5, H and all the »,’s. By virtue of (26), (27), and 
(28), the corresponding eigenvalues will satisfy 


See” hice wk «gh ale Ahn 4(81) 
k 

MeN eee logs ae (32) 
k 


These equations mean that the energy of a state constructed by adding 
phonons to a moving ground state can be regarded as the sum of the 
energy of the ground state and the energy of a ‘ phonon gas’, and the 
momentum similarly. 


§ 5. STATISTICAL THERMODYNAMICS OF THE SYSTEM 


For some systems, the method of constructing stationary states by 
adding phonons to the various moving ground states, described above, 
is sufficient to construct approximations to nearly all the low-energy 
stationary states. In other words, longitudinal sound vibrations are 
the only degrees of freedom effective at very low temperatures. If this 
is so, the specific heat at such low temperatures will be that for the 
equivalent assembly of harmonic oscillators, as given by (1). 

The mean number of phonons in the mode k will be 


Cu, = OXP WNC RL 1) ee Se en (33) 
(cf. Fowler 1936). Since the right-hand side is independent of JN, it 
follows from (25) and (22) that p, and p, have the order of magnitude 
N12, Therefore the approximation (13) is justified for this application. 
The approximation (19) will be justified if »,.,,<1. This will be so if 


Tih chime a ae ee me OS 
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For helium II, if we take k,,=(N/V)13, the right-hand side of (34) is 
about 5°. 

If there are low-energy states which cannot be constructed by adding 
phonons to moving ground states, one may say that the system has extra 
degrees of freedom, other than longitudinal sound modes. In this case, 
the specific heat must be greater than (1). For example, in a crystal, 
transverse vibrations are extra degrees of freedom, which give an extra 
contribution to the specific heat, comparable with the contribution of 
longitudinal vibrations. 

In this way, the theory of phonons gives a lower limit to the specific 
heat, and whether or not the lower limit is achieved shows whether or 
not all the low-energy stationary states can be constructed by the method 
described in the last section. I am grateful to Professor Lars Onsager 
for pointing this out. 

$6. THe LingEAR CRYSTAL 

In this and the next section it will be verified that the present theory 
gives correct results for some simple models whose low-energy stationary 
states are known. The first of these is the linear crystal. Equation (19) 
does not apply to a real crystal because it ignores the coupling between 
longitudinal and transverse vibrations. However, the theory can be 
applied to a linear crystal, where transverse vibrations are excluded by 
confining all the particles to a straight line. 

Consider the part of configuration space where all the quantities 
Gjai—Q—L/N (j=1, 2,...N—1) are small (L being the repeating 
distance of the periodic boundary conditions). This means that the 
particles are approximately equally spaced and are arranged in the order 
in which they are numbered. The deviations of the particles from the 
equally spaced positions are 

r= j7—ILIN—Q,,» -:\ Se 
where oor q/N—3(N+1)L/N. (This choice of ‘equally spaced posi- 
tions ’ serves to make 2’; equal to zero.) 


Jj 
A good approximation to the Hamiltonian for the part of configuration 
space where all the 7,’s are small is obtained by Born and von Karman’s 
method—neglecting third and higher order terms in a Taylor expansion 
of the potential energy : 
N-1 


Hi = Zp #/2m+ (metN*/212) a (rar (ry rah + const. 
d 


j=1 

The interaction between the first and N’th particles arises from the 

mathematical device of making the Hamiltonian periodic, mentioned in § 2. 
The canonical transformation | 


P,.= Lexp (tkjL/N)p,, (all k) 
e] 

rr Lexp (ikjL/N)r,/N, (k 0) | , Ste 
2 
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where k takes the values 


—(N—1)n/L, —(N—3)n/L, ... —2n/L, 0,...(N—1)n/L (N odd), 
—(N—2)n/L, —(N—4)n/L,...—2n/L,0,...Nn/L (N even), 


transforms (36) into a sum of harmonic oscillator Hamiltonians : 


H=i _ (P,,.P_,/Nm+Nmw,2Q,,Q_;,)-+const. 


=JS,2/2Nm+ 2 ho, €,€,+const., 
k40 
where 
w, = (2cN/L) sin| k| L/2N~c| k |, 


§,= (2Nmhw,k?) Vk(P,+iNmw,Q,). 


The creation operators &, are due to Peierls (1929). In terms of Ps 
and Q,,’s, the Fourier components of density become 


py=m Zexp (ikg,) 
j 

=m Lexp (tkjL/N){1+-ikr;+O(kr,)} exp (tkQy) 
i) 


=ikNmQ, exp (tkQ,).... (&40). 


Now, the equations of motion for Q;, and P,, following from the Hamil- 
tonian (36) are, for k0, 


0,=Q,=P,/Nm, 
P,=P,~Nmw,2Q,, 
and hence ppoikP,, exp (ikQ,), 
py=—tkNmw,2Q;, exp (ikQy) = —e*p,, 
(since OO), Thus the sound propagation eqn. (19) holds, and conse- 


quently approximate stationary states may be constructed using the 


creation operators 
Her CORP AUhGa)= | ae er ce -. (38) 


For small values of k, the 7,’s differ from Peierls’ creation operators 
&, only in the factor exp (¢kQ)). This factor is analogous to the operator 
G., of (14), and changes the momentum of the system but not its internal 
- energy, while é, changes the energy but not the momentum. Although 
either factor can be used separately to construct eigenfunctions of the 
approximate Hamiltonian (36), these functions will be approximations 
to eigenfunctions of the true Hamiltonian only if the two factors are 
combined as in (38). If they are not so combined, the boundary condi- 
tions cannot be satisfied. For example, exp (¢kQ,)) by itself does not 
obey the periodic boundary conditions unless ’ either vanishes or is 
outside the range of values enumerated in (37). 

This discussion shows why Peierls’ theory apparently leads to the 
incorrect result that a phonon has no momentum. 
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§ 7. IMPERFECT GASES 


To compare the present theory with the known results for imperfect 
Bose and Fermi gases, the expressions for p, and its time-derivatives 
in the formalism of second quantization (Dirac 1947 b, c, Fock 1932) 
are necessary. These are (Tomonaga 1950, Nishiyama 1951) 


PR=M Was Pa—He ° 5 S J ° 5 : a (39) 


pr=th >) (a . k)ta+jeYa—ak> A e - c . C (40) 


a 


py —(h2/m) & (a. k)2ba4 etha—ae— © (A+ Kk) pape—aUa/m?V, 


where U.=| exp (ia. x)U(x) dx, and %, and 4, are respectively the 


} 
creation and destruction operators for a particle with momentum /k. 


In the sum which forms the second term in the expression for p,, the term 
containing p, is overwhelmingly the largest because (13) applies to only | 
one of its factors instead of to two. Therefore, if U, decreases rapidly 
enough as k-+oo, only this term need be retained : 

pee — (IF) E (ak) Pay reba te —EN py UalmV. . . . (41) 


a 
This is equivalent to Pines and Bohm’s approximation of random phase 
(1952 a). 

In an imperfect Bose gas at low temperatures, nearly all the particles 
are assumed to be in the zero-momentum single-particle state. Since 
the number of particles with momentum /k is 

ye = Db. PS 


this assumption may be written 
popo> 2 bits or o> Ye 


Bogolyubov’s approximation (1947) is to neglect %, in comparison with 
#. Applied to (39), (40), and (41), this gives 
P= (Prot YPotpt); 
Py= Tink (Pyro — Port): 
Pu —(hPK4/ Am) (Pyro + Porte) — FAN py Uy/mV 
~—k*{(hk/2m)?+ U,N/mV }p,. se) Se 
Therefore the oscillation eqn. (9) is satisfied, with 
c?=(hk/2m)?+U,N/mV. 
If the velocity w of the system is negligible, so that py p,, (44) is equiva- 
lent to (19), so that the present theory is applicable. Indeed, Bogolyubov 
showed that all the low-energy stationary states could be constructed 
from the ground state by means of certain creation operators. These are 
just the 7 8 of the present theory, with the approximations (43) and 
Pr=Px- Bogolyubov also showed that the velocity of sound as calculated 


(43) 
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from the compressibility is (U)N/mV)12 ; hence (44) is a genuine sound 
propagation equation for small enough values of h. 

In an imperfect Fermi gas at low temperatures, it may be assumed 
that all the single-particle levels with momenta <hky, where ky is the 
Fermi maximum, are full, while all those with momenta > hk p are empty. 
That is, m%=Yy,.0 and 4,20 if k>ky, while 1-1. = ih, 0 and 
W.=0 if k<ky. Hence either /, or yw, is negligible unless kykp. The 
sums in (39), (40), and (41) accordingly become, for 40 and k<kp, 


prem 3 bas ya ites 


a Sky 
: as -; ae pote ee (4D 
Pein 2 (A. K)bas yha— sk, ea 
a=kp 
puo —(h2/m) = (a. k) bas aha—ye— PU oN /mV. . . (46) 
ak p 


For one dimension, (9) follows at once, with 
C= (hkp/m)?2+U,N/mV. 


Just as for the Bose gas, the theory of §§ 3 and 4 is applicable if the 
velocity of the system is small. 

Tomonaga (1950) has shown that the Hamiltonian for this system can 
be expressed as a sum of harmonic oscillator Hamiltonians. The creation 
operators for these harmonic oscillators are, with the approximations (45) 
and py&p,, just the 7,’s defined in (22). A perturbation calculation of 
the energy of the ground state can again be made; the resulting value 
of the velocity of sound agrees at long wavelengths with the expression 
given above. 

For three dimensions, (k.a)? in the expression (46) for p, may be 
replaced by its value averaged over all directions of a and k, which is 
4k?a?,_ Then (9) follows, with 


=) (hkp/m)2+U,N/mV. 


As before, the theory described in this paper is applicable if the velocitiy 
is small, and the perturbation calculation verifies the formula for the 
velocity of sound in the limit of long wavelengths. However, even without 
interactions formula (1) will give a specific heat proportional to 7’, 
while the true specific heat is known to be proportional to 7. Therefore, 
there must be extra degrees of freedom, and the phonon theory is not 
very useful in this case. . 

For Coulomb interactions, U,,=47e?/k?, where « is the charge on each 
particle, and therefore (41) is approximately 

Py t(47e2N/mV) py, =0 

for small k. This corresponds to plasma oscillations (Tomonaga 1950, 
Pines and Bohm 1952 a). In this case, it is not clear whether or not 
extra degrees of freedom are effective at low temperatures. a 

The application of the theory developed here to the problem of liquid 
helium II will be considered in part II of this paper. 
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ABSTRACT 

A variation principle giving the energy levels of electrons in polyatomic 
lattices is presented. When used in conjunction with the cellular method 
it includes the matching of the wave functions across cell boundaries 
without the need of making ad hoc assumptions. Values of V.(k) can 
be obtained from quantities used to determine E(k) without carrying out 
any numerical differentiation. The results of some calculations on the 
valency bands of lead sulphide are reported. 


§ 1. InTRODUCTION 


Iv is often found that a convenient method of dealing with problems in 
mathematical physics is to set up a variation principle which is equivalent 
to the more obvious equations. It is then possible to approximate to 
some function associated with the problem by a convenient expression 
containing parameters which are varied to give the best answer possible. 
However the Ritz variational form of the Schrédinger equation can only 
be used to determine the energy levels of electrons in crystals if a very 
restricted class of trial wave function is used (as in the tight binding 
method and in the method of orthogonalized plane waves). Such 
functions and their variations must satisfy the correct boundary condi- 
tions on the surface enclosing the volume of integration, and the condition 
appropriate to Bloch functions in crystals gives little freedom of choice. 

Kohn (1952) has removed this restriction by formulating a variation 
principle which includes the Bloch condition as well as the Schrodinger 
equation. As usual, the latter is expressed by a volume integral, but the 
Bloch condition gives a new term, a surface integral over the boundary 
of the volume (which need only be one unit cell of the crystal). However, 
-when trial wave functions are obtained by the cellular method, Kohn’s 
form of the principle applies only to crystals with one atom in the unit 
cell. The present paper extends his variation functional to polyatomic 
crystals and also shows that it is closely related to an expression giving 
the gradient of the energy in momentum space. 

The technique has been used to redetermine some of the energy levels 
of lead sulphide giving results in agreement with those previously obtained 


by other methods (Bell e¢ al. 1953). 


* Communicated by the Authors. 
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§ 2. Tork VARIATION PRINCIPLE 


The total electronic wave function of a crystal can be written as a 
determinantal product of one-electron Bloch functions ¢(k, r), satisfying 
the approximate equation (in atomic units) 

f—4V2+ V(r)—H(k)] d(k,r)=0. . . . . . (1) 
The potential V(r) has the periodicity of the lattice and will be assumed to 


be independent of the form of (k, r). 
When each unit cell of the lattice is divided into N polyhedra it becomes 
convenient to change the phase of xs in each polyhedron and define 
po(k, r)=exp(—ik-p)p(kyr) . . . . . (2) 
where p locates an origin in the pth polyhedron. Because of the Bloch 
condition, the %,, will have the same value at corresponding points in 
equivalent polyhedra in different unit cells. Clearly they will satisfy 
(1), and for ¢ (k, r) to be continuous 
Py (k, r)=exp [vk . (q—p)] be (k, cal Vi, (k, r)=exp [ik . (q—p)] Vika (k,r) 
(3) 


for all points on the boundary between the polyhedra p and q. 
When Kohn’s analysis is applied to these conditions on the y%,, it is 
found that the generalization of his formula (5.8) is that the functional 


Ijp)=Real 2 { | Pol —} V2+ V—E(k) |b, dr+4 ¥ exp [ik . (q—p)] 
p= Tp qd 


x | bo¥%, dS, % i eee 
"pq 

must be stationary for all variations of the ,. In (4), % denotes the 
complex conjugate of %, p goes over all the polyhedra in one unit cell, 
7, being the volume of the pth polyhedron, while the g summation is over 
all polyhedra having a surface o,, in common with polyhedron p. As well 
as being stationary, J vanishes for the correct wave function and energy. 
However, the #(k) found from the principle is not an upper bound for the 
correct value (Kohn, loc. cit.). 

The application of the variation principle will depend on the form of the 
approximation to the wave function, but a very convenient technique is 
based on the conventional cellular method. Each polyhedron is chosen to 
have an atom at its centre which is taken as the origin, so that over a large 
volume it is reasonable to approximate to the potential by a spherically 
symmetric function, V(r,). Equation (1) is then separable and the wave 
function in each polyhedron can be approximated to by the expansion 


pb ,(r)=2'A,, Pys(V ps A); py(Fy, A)J=L,(9,, b») P, (75,A) . . (5) 
where r,, 0, and ¢, are spherical polar coordinates, the L, are angular 


functions depending on k and completely determined by symmetry, and 
the P, are solutions of a radial equation such that 


[—£V?+V(rg)—A]bpa(p AV=O. . 2 2 2. (6) 
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The coefficients A,, (in general complex) and the parameter X are to be 
determined by the condition that the functional J shall be stationary when: 
they are varied. For this type of trial function, the functional can be 
written in the real form 


ee WE) 2A, Ata 2 2A, A Ko, (KoA). . (7) 
p,§ P,§ GU 
Keys qu(k,A)=4 exp [ike (A—P)] | [FoulA) VF (AFA) Vaal A)] - AS, 
=K yu, yal, 2) be 


when, as is always possible, the base functions ¢,,, are chosen to be ortho- 
normal in their own polyhedra : 


| PnP dt=6,,. 
™p 

Equation (7) neglects any lack of spherical symmetry in the potential 
function in each polyhedron ; if a more exact potential is available, the 
¢,; can still be used as reasonable base functions, but the volume integral 
will contribute another term to (7). 

From the condition that J shall be stationary the following equations 
are obtained : 


LO aie ee eee Ketones (Kye \ ee te moun (8) 
0A ,, q,uU ; 
for all ps ; 
oJ oa Ar Olan May 
WH ~ ZAwAnt Z 2 A, Ae| Bett |=0, poaee 25 (9) 


The eqns. (8) are of familiar latent root form and, together with (9) 
and a normalization condition, are just enough to give those values of A 
and the A,, from which the best #(k) can be obtained. 

In principle this may be done by solving (8) for trial values of A and 
then substituting the coefficients in (9) until a value of A which satisfies 
it is found ; in practice this is a very tedious procedure. However, if an 
infinite number of terms could be included in the expansion (5), the 
value obtained for 4 would be H(k) ; it is only the use of a finite number 
of terms that makes the variation of necessary. This suggests that 
choosing }==H(k) will be a good approximation even for a finite number 
of terms. Equations (8) then become 


ZA oy Kgs, quik, A)=0 Meee ee st LO) 
q,u 
and for these equations to be consistent the determinant of the K,, 9, 
must vanish: the value of A for which this happens is the required energy, 
E(k). This simplifies computation enormously but sacrifices an important 
aspect of the variation principle : the compromise between a mismatch 
at the boundaries and a poor solution of the differential equation (1), 
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However, this is a minor approximation compared with others implicit 
in the cellular method. 

The use of (10) introduces an ambiguity for it can be shown that, 
when (10) is satisfied, the cellular functions must obey one of the two 
continuity conditions 


bo=+ exp fik.(q—=p)Wa = ar 


and, in general, only one of these gives the Bloch function required. The 
other gives spurious solutions derived from a variation functional which 
differs from J of (4) by having a negative sign between the volume 
integral and the surface integrals. This fact makes it possible to recognize 
them. But suppose that all the polyhedra in the crystal can be divided 
into two sets such that no polyhedra in the same set have a common 
surface (for instance, crystals with NaCl structure when the two polyhedra 
are taken to be of equal size). Then both continuity conditions in (11) 
correspond to a Bloch function of wave vector k obtained by taking 


f(k, r)= exp (ik.p), for polyhedra in one set, 
o(k, r)=-+ exp (tk .q), for polyhedra in the other set, 


(the sign used for the second set depending on the sign in the continuity 
conditions) and the determinant has no spurious zeros. Now, however 
(10) gives no relation between the coefficients in the two different sets 
of polyhedra. The latent root eqn. (8) does not solve this difficulty 
because it can be partitioned by putting together all elements associated 
with polyhedra of the same set 


ae 2 See pa ETB) 
2 pt [A—E(k) |, 


where J, is the unit matrix of order M, the total number of base functions 
in one set; it is assumed that the total number of functions in each set 
is the same. It is easily seen that at A=H(k) the matrix has two degen- 
erate eigenvectors which may be taken in any linear combination. No 
limiting procedure is possible because the coefficients obtained from (12) 
for AAH(k) satisfy 
we ei 2 2’ ne qu (A) 
iO 2) AeA ee ee 


A _ p,s g,u l . 
pinone 8 qin other u aaa be 2(H(k)—A) ( Z) 
set set 


which implies that for any electron the total charge in all polyhedra of 
the p set is equal to that in all the q polyhedra, a physically absurd 
conclusion, Hence for this class of crystals (10), rather than (8), must be 
considered to be the correct equation, while (9) supplies the required 
relation between the coefficients in the two sets of polyhedra, 
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§ 3. THE GRADIENT oF E(k) tN k Space 


The variation principle of eqn. (4) can be used to obtain a convenient 
formula for the vector y,(k). Since J vanishes for the correct Py 
and is stationary for their variations, a first order error in #, only intro- 
duces a second order error in the energy obtained by equating J to zero. 
Hence it is permissible to assume the ; » to be independent of k when 
calculating V,H(k). This leads to 

N 


p=1" ty 
N : 
=—$Imag. 2 Zq—p)explik.(q—p)l] ha VHp-ASy-- + (14) 
“pg 

It is generally assumed that V is independent of k so that V,H(k) can be 
found in terms of surface integrals which are evaluated when E(k) is 
determined. This is much simpler and more accurate than obtaining 
Vif by numerical differentiation of H(k) and is a very convenient feature 
of this method of finding crystal energy levels. The frequency with 
which y,/(k) appears in formulae relating to the properties of electrons 
in crystals, e.g. the density of states, makes it very desirable to have a 
simple method of calculating this quantity. 


Table 1. Energy Levels of the Valency Band of PbS 


Energy (arbitrary zero) 
Angular 5 


Symmetry Functions 


Variation | Matching conditions 
procedure | (Legendre matching) 


Og) gy Og —0-048 —0-049 


My, (Ya)i> %a —(Q)-111 ety 


—0-179 —0-18 


—0-039 = 

—0-005 —0-004 
—0-005 —0-010 
+0:007 +0-016 
+.0-005 +0-007 


§ 4. APPLICATION TO LEAD SULPHIDE 


Lead sulphide is a crystal of NaCl type so that its energy levels and 
wave functions can be obtained by using eqns. (9) and (10). Values of 
the ,, (eqn. (6)) obtained in an earlier calculation (Bell et al. 1953) have 
been used to redetermine some energy levels in the valency bands giving 
the values shown in table 1 (notation as in Bouckaert, Smoluchowski 
and Wigner (1936)). The energy of the I’, level was calculated with 


H2 
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varying numbers of terms in the expansion of the wave functions in the 
two cells, the notation for the angular functions being as in Von der Lage 
and Bethe (1947). 

When the coefficients in the expansions were calculated it was found 
that (9) did not give a reasonable value for the relation between the 
coefficients in the two polyhedra, lead and sulphur, with the limited 
number of terms used, possibly because the orthonormal relations for 
the ,, were satisfied over equal volume spheres (this approximation has 
no effect on the energy calculated from (10)). Also, because of the 
symmetry of (9) in the two polyhedra, there seems to be no way of deciding 
between two sets of coefficients which give opposite total charges to the 
two cells. To overcome these difficulties a semi-empirical procedure 
was adopted, and the relation between the coefficients in the two cells 
obtained by matching the average values of the wave function and of its 
normal derivative over each face in a unit cell. In general, this gives 
twelve different values for the ratio, one each for function and derivative 
on each of six faces, and deviations among these values give some estimate 
of the accuracy. This is illustrated in table 2 which gives the coefficients 
for the level I’, with varying numbers of terms in the expansion of # in 
each cell. | 


Table 2. Coefficients of Base Functions for J’, Level 


Pb-cell S-cell 
34 33 55 8’; 3, 35 55 0’; 
0-473 | 0-233 —0-847 | —0-075 a 
Variation | 0-472 | 0-228 | —0-003 —0-848 | —0-075 0-005 
Procedure | 0:425 | 0-254 —(0-046 | —0-864 —0:077 0-027 


0-431 | 0-285 0-020 | —0-035 | —0-849* | —0-099 | —0-012 | 0-027 


0-465 | 0-222 —0-852 | —0-094 
Legendre | 0-475 | 0-225 0-007 —0-845 | —0-103 | —0-009 
Matching | 0-401 | 0-255 —0-039 | —0-874 | —0-094 0-028 
0-429 | 0-259 0-013 | —0-037 | —0-857 | —0-110 | —0-011 | 0-027 


* = . : . 
There are only two non-equivalent faces for this value of k, and the ratio 
given between the 6, coefficients in the two cells is the mean of 


equating the average of % on x=constant : —1-962, 
equating the average of 4, on x=constant : —1-986, 

equating the average of ys on y=constant or z=constant : —1-998, 
equating the average of +, on y=constant or z=constant : —1-947, 


It is clear from tables 1 and 2 that the matching conditions employed 
in previous work (in this case matching of Legendre coefficients, see 
Bell et al., loc. cit., p. 82) and the variation procedure lead to practically 
the same results ; it appears that the variation procedure converges 
more rapidly to the correct energy. It should be mentioned that if very 
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few terms are included in the wave function, the determinant of eqn. (10) 
may not vanish for any value of 2. 

As an illustration of the calculation of the slope of the H—k curves 
by eqn. (14), table 3 gives some values of dH (k)/dk, for the 4, valency 
band at various points in the («, 0, 0) direction. 


Table 3. Values of d#/dk, 
PbS—Valency band. k=(«, 0, 0). 4, symmetry 


| S—L| o 0 0:74 | 1-151 1-341 1-596 1-81 7 
#H —()-004| 0:03 | 0-06 0:06 0-03 0 —0-11 
From eqn. (14) 0 0-37 | 0-14 | —0-53 | —0-77 —0-74 0 
From H-k curve] 0 0:4 0-2 —0-3 —0:8 —0:8 0 


§ 5 ConcLUSION 


The variation principle discussed in this paper throws some light on 
the method used by Howarth and Jones (1952) in their calculation of 
the energy levels of sodium. They obtain their best values of A and 4,, 
by point matching conditions on the boundaries of the unit cell, and then 
substitute in what is effectively (7) to determine E(k) by the condition 
that J shall vanish. The difference between X and H(k) is used as a 
criterion for the accuracy of the wave function. In the present paper it 
is suggested that the stationary property of J gives a better method of 
obtaining the energy and coefficients than any necessarily arbitrary 
matching conditions. With the same number of terms in the expansion 
of s, the numerical labour is greater than that required in point matching : 
it has been found, however, that a smaller number of terms is necessary 
to obtain the same degree of approximation. Because of this and of the 
ease of obtaining V,/(k), the variation procedure, besides being more 
general and rigorous, appears to be of considerable practical value. 
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XIII. CORRESPONDENCE 


The Propagation of Heat in Liquid Helium 


By J. M. Ziman 
Clarendon Laboratory, Oxford* 


| Received October 30, 1953] 


THE results of experiments to measure the velocity of second sound in 
He II below about 0-:5°K have not proved easy to interpret. It has been 
suggested (Gorter 1952, de Klerk, Hudson and Pellam 1953, Atkins 1953, 
Mayper and Herlin 1953) that true second sound does not exist at the very 
lowest temperatures in pure *He, and that the observed propagation of 
heat pulses is due to the scattering of phonons by the walls of the appar- 
atus. There does not, however, seem to have been any quantitative 
consideration of the effects to be expected if this interpretation is correct, 
for comparison with experiment. This may be because of the mathe- 
matical complexity of the problem, and the labour required to obtain a 
general exact solution. But it may not be generally realized that a 
reasonable approximate solution can be obtained very simply for the case 
of a tube which is long compared with its diameter. 

The essence of the argument is that we may neglect those processes in 
which the phonons reaching the detector have made few collisions with the 
walls of the apparatus. Since their wavelengths are small, the phonons 
will presumably be scattered at random, and the problem can be reduced 
to a complicated variant of the problem of the ‘ random walk ’, with the 
appropriate probability distribution functions. But it is known 
(Chandrasekhar 1943) that, in the limit of many collisions, this problem 
may be reduced to the solution of a diffusion equation, 

oT ord 
Tira Oe 


for the temperature 7, which is equivalent to the phonon flux, at any 
point of the surface of the tube. This equation may evidently be identi- 
fied with the equation of heat conduction, and comparison with Casimir’s 
analysis (1938) for the boundary scattering thermal resistance of a cylinder 
of radius R, shows that D=2vR/3, where v is the mean phonon velocity 
in the medium. Again, as a reasonable approximation, to avoid mathe- 
matical complications, we consider an infinite tube, ignoring the end effects 
which would presumably correspond to scattering of the phonons back 
and forth along the length of the tube, i.e. reflection of the received pulse. 


EE ee ee eee 


* Communicated by M. H. L. Pryce. 
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It is then easy to show (e.g. Sneddon 1951) that the temperature 7'(x, t) 
at time ¢ at a point distant x from the origin is given by 


2 (co) 
Te, t) = — | T(0, t—x?/4Dn?) e-” dn, 
eH) Va ere ieee) m 


where 7'(0, t) is the temperature at the origin at time ¢ (+0) and the tube 
is initially at zero temperature throughout. Fora square pulse, of height 
7, and length t), applied at the end of a section of length L, we have 

T(L, t)=T olerf {(r—r9)-1/?} —erf {r-/2}], 
where 7—8fRvt/3L? and 7,=—8Rvt,/3L?. In the limiting case of a very 
short pulse, we obtain 

T(L, t) « 7-8/2 exp (—1/7). 

It will be seen from the figure, in which the cases ry > 0 and 77=0-2, 0-5 
have been plotted, that the pulses are similar in shape (having been 
reduced to the same maximum ordinate). 


l 2 a a ya ; 


Pulse shape calculated for 75>0, t7=0°2, 79> =0°5 
— —~—- Observed pulse shape (Atkins and Osborne). 


For comparison, the shape of the pulse observed by Atkins and Osborne 
(1950) is also plotted on the same scale, and it will be seen that it agrees 
quite closely with that predicted with the present theory. In their case 
we have L—34cm, 2R=0-5cm, and 4=100 HSeC, 80 that To= 0-14 
(taking v equal to 2:37 x 10*cm sec"?). The leading edge lies very close 
to the curve which has been plotted for t)>=0-2, but there is a considerable 
discrepancy in the shape of the tail, although here the experimental curve 
is difficult to measure against the background of noise. 

The approximation is valid only when 2A/L is rather small ; otherwise 
the main contribution to the observed pulse comes from phonons which 
have made very few collisions, and their arrival times are no longer 
governed mainly by the random scattering process. In this case, also, 


102 Correspondence 


the replacement of the exact binomial distribution in the solution of the 
random walk problem by continuous functions derived from Stirling’s 
approximation tends to bring the leading edge too far forward. An 
alternative approximation is easily obtained for the other limiting case, in 
which the tube is very short compared with its diameter, so that most of 
the phonons pass directly between the source and the detector, but this 
configuration does not seem to have been used in any of the experiments. 


This note arose out of discussions with Professor H. A. Fairbank and 
Dr. J. Wilks. 
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Optical Determination of the Thickness of the Unsaturated Helium II Film 


By N. G. McCrum and K. MENDELSsoHN, F.R.S. 
Clarendon Laboratory, Oxford * 


[Received November 12, 1953] 


THE thickness of the saturated helium film has been measured by Burge 
and Jackson (1951) using an optical method based on the influence of 
surface layers on metallic reflection. In this laboratory volumetric 
experiments have recently been made on the isothermal formation on 
' glass surfaces of the unsaturated film with increasing pressure (Brewer 
and Mendelssohn 1953). These have now been supplemented by optical 
measurements with percentage saturations varying between 60 and 100%, 

The method used was essentially similar to that of Burge and Jacksont 
(loc. cit.), but required certain modifications in our case. In order to 
study the formation of the film, the mirror had to be enclosed in a 
separate adsorption chamber which in turn was immersed in a bath of 
liquid helium IT of constant temperature. Since the unsaturated film 
is considerably thinner than the saturated, great care had to be taken 
to screen it, when under measurement, from extraneous light. The 
beam was admitted to the apparatus by means of a small shutter in short 


bursts usually of less than two seconds duration. Balances were always 
¥ eer Sie i Se a et 
* Communicated by the Authors. 
} We are greatly indebted to Dr. L. C. Jackson and Dr. E. J. Burge for 
much advice. 
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finally checked by observing the instantaneous half field matching on 
opening the shutter. This procedure was repeated several times for 
each point, the experimental scatter in the analyser nicol reading being 
usually under +2’ of are. The illuminated portion of the mirror was 
moreover only a small fraction (~10~*) of the total surface carrying the 
film. Thus any evaporation of film from the observation area could 
be made good rapidly by superflow from the unilluminated part. Since 
the percentage saturation in the adsorption chamber is a rapidly varying 
function of temperature, the latter had to be kept accurately constant. 
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The difference between the analyser reading and an arbitrary fixed position 
against relative saturation. The quarter wave plate fast axis orienta- 
tion with respect to the vertical was 134° 33’ at 1-84°K and 127° 10° at 
1-87°K. 
This was particularly essential close to saturation where the film thickness 
is strongly dependent on pressure. By keeping the temperature constant 
to better than 0-0005°K during all observations leading to an experimental 
point, reproducible results were obtained. 

The results for four different temperatures expressed in nicol rotation, 
are plotted in figs. 1 and 2. As these graphs show, the thickness 
was found to be a rapidly increasing function of relative saturation 
(p/po) for 0-95<p/pyo<1-0. This is in good agreement with the volu- 
metric adsorption measurements (Long and Meyer 1949, Frederikse 
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and Gorter 1950, Brewer and Mendelssohn 1953, Strauss 1953 (cf. 
Long and Meyer 1953)), which though differing in other respects all give 
a very steep increase in coverage just below saturation. The same 
behaviour was found by Bowers (1953), using a gravimetric method. 
It is interesting to note that measurements of superflow of the unsaturated 
film (Bowers, Brewer and Mendelssohn 1951) indicate that even at 
the lowest temperatures measured (1:37°K) the ratio of the unsaturated 
to the saturated flow rate is less than 0-1 for p/pyp<~0-9. Any tempera- 
ture dependence for p/p)<1-0 for values of 77<7', was too small to be 
resolved in these experiments. A more detailed account of the work 
which is being continued will be given at a later date. 


Fig. 2 
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Slip Lines, Etchpits and Corrosion-pits 


By A. H. Appou 
Department of Metallurgy, Sir John Cass College, London, E.C.3* 


[Received October 29, 1953] 


Ir has been shown recently that by rapid etching of crystals the centre 
of each isolated spiral growth arising from an independent screw disloca- 
tion is replaced by an etchpit (Horn 1952). This technique for revealing 
dislocations in crystals depends on the differential rates of chemical 
attack of a dislocation and good crystal. 

Before this relation between a dislocation and an etch pit was discovered, 
the grouping of etchpits was used to reveal crystallite boundaries 
(Lacombe 1948) and polygonization walls (Cahn 1949). Amelinckx (1953) 
was the first to attempt the use of this technique in finding correlation 
between dislocations and slip lines in aluminium. He found that etchpits 
grouped themselves preferentially along slip lines and certain grain 
boundaries in addition to a general random distribution. 

Similar observations have been made by the writer on single crystals 
of high purity zinc (99-99°) grown from the melt. These were electro- 
lytically polished in 20°/, chromic acid and then deformed in tension. 
Some crystals were left in the laboratory atmosphere and others were 
etched in 20%, chromic acid for various periods. In the first case it was 
found that corrosion pits rapidly developed on the surface and their 
positions can be distinguished along slip lines and at random. Figure | 
(figs. 1-4, Pl. 1) is of a crystal left for 12 hours in the laboratory 
atmosphere and fig. 2 is of another crystal after 36 hours. 

In the case of etching, however, it was found that if the attack is 
rapid enough etchpits group themselves along slip lines only leaving the 
rest of the surface free of them as shown in fig. 3. Figure 4 shows the 
effect of prolonged etching ; existing etchpits on slip lines are deepened 
while new ones develop at random. 

This work is now being extended to find a correlation between the 
distribution of etchpits and the slip displacement and to investigate 
whether slip originates near the surface or in the interior of the crystal. 
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XIV. Notices of New Books and Periodicals received 


Crystal Growth and Dislocations. By Asir Ram _ VERMA. (London 
Butterworths Scientific Publications, 1953.) [Pp. 182.] Price 30s. net. 


Errors and omissions excepted, this is a valuable and timely little book. It is 
based on a Ph.D. thesis, and displays the merits and demerits of such works. 
The portions of the subject with which the author was most intimately concerned 
are generally well presented. The remainder is compiled from the best authori- 
ties, and usually condensed merely by the omission of steps, which does not 
introduce error, but may leave the reader uncertain about the degree of 
reliability of a conclusion. When the author adds or substitutes his own 
arguments, the result can be disastrous. Like others before him (Wulff 
included), he devises a short ‘ proof’ of Wulff’s theorem which is no proof at 
all. His particular fallacy is the assumption that two n-dimensional vectors 
orthogonal to a third are necessarily parallel to each other. There is fortunately 
a reference to the correct treatment by von Laue (1943). Minor errors of 
argument occur not infrequently : however the statement on page 28 that a left- 
handed screw can be changed into a right-handed screw by a rotation should 
serve as a sufficient reminder that critical reading is necessary. The main 
point. is that there has recently been a decisive advance in the understanding 
of crystal growth, and that silicon carbide, the principal material of Dr. Verma’s 
researches, exhibits a greater variety of the newly comprehended phenomena 
than any other substance, together with a few that have yet to be understood. 
This documentation of the subject should be of real value to all for whom 
crystal growth or crystal surfaces are of importance. 

The principal blemish of the book appears to be accidental. It contains a 
considerable number of photographs and drawings provided by other authors. 
With one exception their names do not appear in the legends, and they cannot 
in general be easily discovered from the text.* The preface states that a complete 
list of acknowledgments will be found on p. xii. On examining this page (which 
carries no names) it appears that it is a list of Journal references for the author’s 
own re-published illustrations, with the addition of only two references to the 
work of others. Presumably the wrong list has been printed. Publisher and 
author must both be held to blame for allowing this very regrettable lapse to 
occur, E.C.F. 


The Mathematical Theory of Non-uniform Gases. By S. CHapMaNn and T. G. 
conan 2nd edition. (Cambridge: University Press.) [Pp.431.] Price 
3. 
THE second edition of this standard text book in theoretical physics is almost 
identical with the 1939 edition, apart from 25 pages of added material to 
account for more recent developments. The added pages deal with the viscosity 
on. the basis of the 7-7, r~183 model, with various aspects of diffusion and thermo- 
diffusion, with a third order approximation of the velocity distribution function. 
They also contain some consideration of the work of Kirkwood and Born and 
Green on the theory of liquids. De 


Radioactivity and Radioactive Substances. By Sir JamES CHADWICK. [Pp. 120.] 
(Isaac Pitman.) Price 12s. 6d. 


Tae tothe ie ; 
Tuts is the fourth edition of a simple, concise and accurate statement of the 
main facts and theories of the subject. It first appeared in 1921, but has 
now been revised and supplemented by Professor J. Rotblat. 


* ~P, oher’o N. 2 Mhia por , 7 i i 
* Publisher's Note,—This error was noticed and the necessary correction slip has been 
printed for insertion in the book, 
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Scintillation Counters. By J. B. Brrxs. (Pergamon Press Ltd., London.) 

[Pp. vili+142+index ; 59 diagrams, 4 plates.] Price 21s. 

On the jacket it is claimed that this book forms a complete review, that all 
factors are discussed, that a new theory is presented in detail, that every major 
application is reviewed, that the references are exhaustive and completely up 
to date and that the treatment is wholly new. However, the reader will be 
relieved, on proceeding further, to find that the book takes the normally 
expected form, except perhaps in a limited number of respects. Thus it is 
possible that the phrase “ wholly new” refers to the author’s advocacy in 
chapters 6 and 7 of the so-called photon exchange mechanism of energy migra- 
tion in the fluorescence of organic substances. Here the author’s intention is 
to propagate vigorously the view that this mechanism (one of four recognized 
mechanisms of migration) and no other is involved. The possibility that two 
or more of the mechanisms might occur simultaneously in the rather complex 
phenomena of fluorescence apparently is not considered. It is difficult to agree 
that the author’s views are established. He dismisses rather too readily work 
which does not fit into the proposed scheme. He has given a detailed account 
of his own experiments, performed and to be performed. 

Where the author is not concerned with photon exchange the book is well and 
concisely written with numerous, but not exhaustive, references. The first 
five chapters constitute a good condensed account of the main factors involved 
in scintillation technique. There are, however, dangers to the uninitiated in 
this condensation as for example in the diagrams for the energy distributions 
of secondary ionizing particles produced in scintillators by y-rays and fast 
neutrons. These depend greatly on the size of the scintillator which may range 
in practice from a few mm? to several cubic feet. Just over three pages are 
specifically devoted to consideration of the luminescence process in inorganic 
phosphors, a difficult and extensive field meriting longer discussion. Inform- 
ation on crystal growing is scanty and there is no separate treatment of circuitry 
problems arising in the use of very fast crystals and liquids. It is doubtful if 
the somewhat controversial treatment of energy migration can be justified in 
view of such brevity. 

The book is for the most part a useful and compact survey of much of the 
field of scintillation technique and its presentation of factual material is sound 
if somewhat over-condensed. When the author departs from such presentation 
to suggest improvement in the work of others the result is less happy as on 
page 117 where it appears that he misunderstands the method of fast neutron 
spectrometry under discussion. 

“The miraculous speed of publication ” is contrasted with that of certain 
scientific journals but such a criticism of the journals will appear unwarranted 
to most readers at this time. In any case the problems are rather different in 
the two cases and comparison will be valid only if each member of the series of 

nographs appears equally speedily. 

The Rast is likely sane peal to nuclear physicists rather than to those in 
other sciences where the chief task, generally speaking, consists in using the 
technique to its best advantage in a particular set of ee a : 


> 


Methods of Mathematical Physics. Vol. I. By R. Courant and D,. HILBERT. 
(Interscience Publishers, Inc.) [Pp. 561.] Price $9.50. 


THis is a very welcome translation (and slight revision) of the well-known 
German book. It treats such topics as Linear Integral Equations, the Calculus 
of Variations, and Eigenvalue Problems in a manner which is useful to both 
the mathematician and the physicist. It is to be hoped that the translation 
of Volume II will also be available soon. R. W. 
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The Theory of Metals. By A. H. Witson. 2nd _ edition. (Cambridge : 
University Press.) [Pp. 346.] Price 45s. 


Tuts book looks altogether different from the first edition dating from 1936. 
A great number of recent developments in the theory have been incorporated, 
causing a considerable increase in size of the volume. The chapters on the 
formal theory of conduction and the mechanism of conductivity have benefited 
from Dr. Sondheimer’s contributions to the theory. Conduction in semi- 
conductors receives due attention. A chapter on ferromagnetism gives an 
account of the collective electron theory of ferromagnetism. The diamagnetism 
of free electrons is dealt with in a new way, but one has to consult the first 
edition for the expressions for the diamagnetism of electrons in a periodic 
field. A minor point for criticism in this part of the book is that no mention 
is made of the importance of the ‘ thermodynamics of irreversible processes 
in connection with certain symmetry relations and the derivation of the Kelvin 
relations in thermoelectricity. 

Since the emphasis in Wilson’s book is on the formal rather than the physical 
aspect of the theory, it is somewhat disappointing that no great improvements 
have been made in the chapters dealing with the band theory of solids. The 
systematic discussion of the zone-theory of three-dimensional lattices on the 
basis of the reduced zone picture is completely absent. Brillioun zones are 
discussed with reference to structure factors, but it is not clearly pointed out 
that this is only significant for those solids in which the free electron approxi- 
mation is a good one. Strictly speaking the discussion of the diamond structure 
is incorrect, while the suggestion implied in the information on the (non integral) 
number of electrons in the first zone of more complicated structures is very 
misleading. 

The optical properties of metals are not discussed in the second edition and 
superconductivity is also left out. This is in line with the apparent intention 
of the author, not to touch upon subjects of which the present theory is still 
highly unsatisfactory. There is no doubt that the edition will be a most useful 
textbook, especially for the detailed theory of conductivity. Ps 


Chemical Constitution. An Introduction to the Theory of the Chemical Bond. 
By J. A. A. Kerecaar. [Pp. 398.] (London: Elsevier Publishing Co.) 
Price 40s. 


Tats book is a very good translation by Dr. L. C. Jackson of the second Dutch 
edition (1952). Although the contributions of the Dutch school of physical 
chemistry to the study of chemical combination has most notably been con- 
cerned with ionic bonding, the present book is not overweighted in this direction. 
The nature of ionic bonding is, of course, well treated, and covers a variety of 
topics, including ionic complex formation, volatility, hardness and solubility. 
Quantitative illustrations are largely restricted to those instances where the 
bonds may reasonably be regarded as purely ionic. 

The chapters on covalent bonding are rather more conventional and include a 
discussion of colour, free radicals, and chemical reactivity with a brief mention 
of recent work on electron deficient structures. 

The last two sections of the book are concerned with metallic and van der 
Waals binding. The former is briefly but clearly treated with some reference 
both to alloys and to semi-conductors. The latter is discussed very well and is 
related to the subjects of molecular compound formation, melting and boiling 
points, miscibility and solubility. Hydrogen bonding is included in this section. 

In general the presentation is likely to stimulate the reader to further thought 


and, at a not excessive price, the book can be recommended to both honours 
students and research workers. G. E. C. 
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Introduction to Solid State Physics. By ©. Krrreni.. [Pp. 396.] (London : 
Chapman and Hall.) Price 50s. 


THE scope of this book can best be illustrated by a summary of its contents. 
The first 90 pages deal with crystal structure, lattice energies and thermal 
capacities ; 110 pages follow on dielectric and magnetic properties of crystals ; 
100 on electrical properties ; 50 on imperfections in solids, including firstly 
Frenkel defects. and related topics, and secondly dislocations and mechanical 
properties. Finally there are 50 pages of Appendices elaborating some of the 
more mathematical points. 

In his preface, the author gives an impressive list of seventeen topics that have 
been omitted, and explains that those subjects have been favoured which can 
be discussed in terms of simple, concrete, and well-developed models. Even 
so, the result is a highly condensed piece of writing which can be called an 
‘ Introduction * only in virtue of an excess of modesty on the part of the author. 
Lack of space has prevented the inclusion of the experimental details relating 
to the phenomena dealt with. But further than this, it is assumed that the 
reader has a fair knowledge of the results of such experiments, and the discussion 
proceeds almost immediately to a consideration of models which will account for 
these results. The emphasis here is rather on the mathematical development, 
than the underlying physical concepts. 

The book is clearly written and well documented ; it will prove invaluable to 
many readers—not least to those who have to prepare lectures on its subject 
matter. But to imply, as is done in the preface, that it constitutes a text-book 
suitable for students, either just before or just after graduation, is to suggest a 
standard of attainment in American Universities a good deal higher than in 
their British counterparts—a conclusion which could not be supported by most 
people in a position to judge. 


The Outer Layers of a Star. By R. v. bp. R. WooLtEy and D. W. N. Srisss. 

[Pp. xi+303.] (Oxford: Clarendon Press, 1953.) 40s. net. 

A CONSIDERABLE body of knowledge has grown up, especially during the past 
thirty years, about those outer layers of the stars which we can see. The subject 
shares most of the complications of meteorology, as well as the characteristic 
difficulties of astronomy, but there is by now general agreement on much 
observational material and on the usefulness of a substantial body of theory, 
although most workers in this field would readily admit that the real stellar 
atmospheres are far more complicated than our present pictures or models. 
For many years there has been a need for a general book in English on the 
subject, not too difficult or specialized, and this volume may go some way 
towards meeting the demand. 

The mathematical techniques involved in the theory of stellar atmospheres 
have already been thoroughly treated in two other of the Oxford Monographs on 
Physics, those by Chandrasekhar and by Kourganoff and Busbridge respectively. 
Although The Outer Layers of a Star pays more attention to the physical 
problem than these other books, the general emphasis in my opinion could 
with profit have been shifted still more in this direction. Important progress in 
this subject during the past twenty years has come predominantly from new 
physical methods and ideas, and the big obstacles to further advances are still 
essentially physical. The book deals with a great many of the standard 
problems, although certain irregular phenomena such as sun-spots, flares, 
prominences, etc. have been specifically excluded. There is surprisingly little 
attention paid to the question of the chemical composition of stars, although 
spectrum analysis of their outer layers is our chief source of extra-terrestrial 
information. ‘The style of writing is clear and there are sufficient references to 


guide the reader to the more important literature. The printing Pe oorie 
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“Experimental Nuclear Physics, Vol. I. Edited by E. Seeré. (Chapman and k 
Hall.) [Pp. 600.] Price 96s. + 

Tuts is the second volume of a comprehensive treatise aiming to bring the | 

experimentalist up to date in techniques, and indicate the broad lines of 

theoretical interpretation. Its bibliography is complete up to the end of 1952 

in the two fields with which this volume is concerned, namely (i) a survey of 

nuclear reactions, and (ii) the neutron. The third volume is announced to be 

in preparation. 


Flames, their Structure, Radiation and Temperature. By A. G. GAYDON and 
H. G. WotrHarp. (Chapman & Hall, 1953.) [Pp. 352.] Price 55s. 


FLAMES are regions where energy changes form as rapidly as Proteus, and they 
are equally fascinating to the scientist and the engineer. The chemist wants _ 
to know the steps in the reactions, the intermediate products and the degree to 
which equilibrium is reached, the physicist wants to understand the way energy 
is released and distributed and the effects of the flame on flow. The engineer 
on the other hand wants to be able to produce a given heat release in a certain 
space which may be large or small, and to give the resulting flame certain heat 
transfer characteristics which may be high or low emissivity coupled with high 
or low convection coefficients. Gaydon and Wolfhard look at flames primarily 
as physicists, secondarily as chemists and tertiarily as engineers. For this 
reason they are most interested in the stationary non-turbulent premix flame 
typified by the Bunsen burner, because it can be most readily explored by 
physical techniques. Of these techniques, spectroscopy has in their hands made 
an outstanding contribution to the physical understanding of the energy release 
processes, especially when applied to premix flames at low pressures which have 
a much thicker flame front than those at atmospheric temperatures. Laminar 
diffusion flames have a chapter and the effect of turbulence on both premix and 
diffusion flames are considered, but clearly the difficulties in the way of obtaining 
fundamental understanding of the physical and chemical processes are much 
greater in these cases than with the simple premix flame, which has a character- 
istic burning velocity. Propagating and explosion flames are not covered at all 
except where they are of interest as providing another method of measuring 
burning velocity. 

An illustration of the authors viewpoint is provided by the fact that the 
combustion of powdered coal which is responsible for a large and increasing 
proportion of all the electricity made in this country, is mentioned (p. 310) under 
the heading ‘‘Some unusual flames’’. Nevertheless, no engineer can afford to 
ignore the very real progress in understanding which is coming from the study of 
the simple Bunsen flame, while physicists must welcome the fact that flames are 
emerging from the natural history stage and becoming regions to which physical 
formulae can be applied as instructively as they were applied to simple gases at 
very low pressures in the first decades of this century. | 

Gaydon and Wolfhard’s book can be strongly recommended as an excellent 
review of our present state of knowledge in its chosen field, both describing in 
fair detail the experimental methods used for producing the flames, measuring 
flow temperature (two chapters), spectral radiation characteristics, and ion 
concentrations and the theories for flame shape, reaction stages, flame propa- 
gation, soot formation, and energy liberations. The experimental evidence is 
carefully studied and wherever possible conclusions are drawn as to which theory 
gives the best agreement with this evidence, although naturally, in such a 
developing subject many rival theories remain in the field. M. W. T. 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 


